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Recently a new class of quantum phases of matter: symmetry protected topological states, such 
as topological insulators, attracted much attention. In presence of interactions, group cohomology 
provides a classification of these [X. Chen et al. . larXivrH 06.4772v5 (2011)]. These phases have short- 
ranged entanglement, and no topological order in the bulk. However, when long-range entangled 
topological order is present, it is much less understood how to classify quantum phases of matter in 
presence of global symmetries. Here we present a classification of bosonic gapped quantum phases 
with or without long-range entanglement, in the presence or absence of on-site global symmetries. 
In 2-1-1 dimensions, the quantum phases in the presence of a global symmetry group SG, and 
with topological order described by a finite gauge group GG, are classified by the cohomology 
group H'^{SG X GG,U{1)). Generally in d-|-l dimensions, such quantum phases are classified by 
H'''^^{SG X GG,U{1)). Although we only partially understand to what extent our classification 
is complete, we present an exactly solvable local bosonic model, in which the topological order 
is emergent, for each given class in our classification. When the global symmetry is absent, the 
topological order in our models is described by the general Dijkgraaf-Witten discrete gauge theories. 
When the topological order is absent, our models become the exactly solvable models for symmetry 
protected topological phases [X. Chen et al., arXiv:1106.4772v5 (2011)]. When both the global 
symmetry and the topological order are present, our models describe symmetry enriched topological 
phases. Our classification includes, but goes beyond the previously discussed projective symmetry 
group classification. Measurable signatures of these symmetry enriched topological phases, and 
generalizations of our classification are discussed. 
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I. INTRODUCTION 

Recently there has been significant interest in topo- 
logical phases of matter, which are quantum phases of 
matter beyond the Ginzburg-Landau symmetry-breaking 
description^. After the discovery of fractional quan- 
tum Hall states, the notion of topological order was 
proposed^. Topologically ordered phases of matter fea- 
ture ground state degeneracies on torua^, and anyonic 
quasiparticle excitations in the bulk in 2-1-1 dimensions^. 
These features are robust against arbitrary local pertur- 
bations. Therefore, the global symmetry is NOT a re- 
quirement for topologically ordered quantum phases. 

More recently, symmetry protected topological (SPT) 
phases have been discovered. SPT phases are defined to 
have no topological order in the bulk (and thus no anyons 
in the bulk nor ground state degeneracies on torus); nev- 
ertheless their distinctions are protected by the global 
symmetry^"— . Examples of SPT phases include topo- 
logical insulators and superconductors^"—. One exper- 
imental signature of the SPT phases are the symmetry 
protected gapless boundary states, which can be obtained 
from a Chern-Simons based classification of phasesi^.. 
Another classification of interacting bosonic SPT phases 
for on-site global symmetries is provided in the origi- 
nal work by Chen et al."'^^ using the group cohomology 
method. At the superficial level, there is no relation 
between the SPT phases and the topologically ordered 
phases. For example, in l-|-ld, it can be shown that in 
the presence of interactions, there is no topological or- 
der but there are nontrivial SPT phasesi^i"— , such as 
the AKLT— 1^ integer spin chain. A recent beautiful 
work^l shows that SPT phases and topologically ordered 
phases are related via certain duality in spatial dimen- 
sions higher than one. 



How to understand/classify gapped quantum phases 
when both the topological order and the global sym- 
metry are present? This is an important question for 
both fundamental and practical purposes, and this pa- 
per is an attempt to answer it, at least partially. To 
illustrate the importance of this question, as an exam- 
ple, we can consider one famous topologically ordered 
phase: the Laughlin's = 1/3 fractional quantum Hall 
liquid^ (FQHL), which has three- fold ground state de- 
generacy on torus^ and anyonic quasiparticle excitations 
in the bulk. In the physical realization of the Laughlin 
FQHL in 2DEG, there is also a global symmetry: the 
C/(l) charge conservation for electrons. One can imag- 
ine what would happen if the C/(l) charge conservation 
was absent, for instance, if a small electronic pairing was 
introduced via proximity effect. Because the topological 
order is robust towards arbitrary perturbation, the three- 
fold ground state degeneracy and the anyonic statistics 
of quasiparticles would still be present. 

Is the C/(l) global symmetry unimportant for the 
FQHL physics then? Obviously, this is not the case. In 
fact, this C/(l) symmetry allows one to find two striking 
experimental signatures of Laughlin's state: the quan- 
tized Hall conductance iJxy = e^/3/i, and the e* = e/3 
fractional charge carried by quasiparticles. The second 
signature is very interesting: the quasiparticles of a topo- 
logically ordered phase can carry a fraction of the quan- 
tum number of the fundamental degrees of freedom (elec- 
trons here) in the quantum system. Such phenomena are 
often referred to as "symmetry fractionalization" . This 
phenomena only occur when the system has topological 
order. The e* = e/3 charge of quasiparticles is a remark- 
able demonstration of how the global symmetry can "act" 
on the topological order in a non-trivial fashion^Sri^. 

Another collection of fascinating quantum phases is 
the quantum spin liquid (QSL). Quantum spin liquids 
are often defined to be featureless Mott insulator phases, 
namely phases that respect full lattice symmetry as 
well as the SU{2) spin rotational symmetry, with a 
half-integer spin per unit cell. Based on the Hast- 
ings' generalization's of Lieb-Schultz-Mattis theorem^ 
in higher dimensions, we know that gapped quantum 
spin liquids in two and higher spatial dimensions must 
host non-trivial ground state degeneracies on torus. But 
because there is no symmetry-breaking-induced ground 
state degeneracy, this indicates that the gapped QSLs are 
topologically ordered. 

How can one classify/understand QSL phases? For in- 
stance, recent numerical simulations'^'* point out that the 
spin- 1/2 Heisenberg model on a Kagome lattice hosts a 
gapped quantum spin liquid phase. It is then an impor- 
tant issue to understand the nature of this QSL phase. 
As a matter of fact, numerical evidence for topological or- 
der described by a Z2 gauge theory has been found^i^^. 
Is this topological order enough to determine the nature 
of this QSL phase? The answer is negative. It turns 
out that there are more than one QSL phase on the 
Kagome lattice even for a given Z2 topological order—"—. 
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Their distinctions are protected by the global symme- 
tries. Roughly speaking, the way that the global symme- 
tries act on the topological order are different for differ- 
ent phases. These phenomena have been called "symme- 
try enriched topological phases" or "symmetry enriched 
topological order"—"—. When the global symmetries are 
absent, all these phases are no longer distinguishable and 
are adiabatically connected to one another. But when 
the global symmetries are present, one necessarily en- 
counters phase transitions while going from one phase to 
another. Therefore, for the Kagome lattice gapped spin 
liquid example, it remains an unresolved issue to under- 
stand which among all the symmetry enriched topological 
phases is the one found in the numerical simulations. 

The above physical examples motivate us to consider 
the following questions: How are symmetry enriched 
topological (SET) phases generally classified? Or, how 
can one classify different ways in which the global symme- 
try "acts" on the topological order? What are the exper- 
imental/numerical signatures of different SET phases? 
The last question is quite urgent for the above Kagome 
QSL example: although there are nice numerical meth- 
ods (e.g., the topological entanglement entrop y'^^'^'' ) to 
detect the Zi topological orde r^^i^^" — , due to the lack 
of theoretical understanding it is still unknown how to 
numerically distinguish different SET phases. 

This paper attempts to address these questions to a 
certain level. We consider on-site global symmetries only; 
namely the global symmetry transformation is a direct 
product of unitary transformations, and each transfor- 
mation only acts in the local Hilbert space. In addition, 
we focus on bosonic systems with finite unitary symme- 
try groups SG and topological orders that can be de- 
scribed by finite gauge groups GG. Generalizations of 
these conditions will be discussed at the end of the pa- 
per. Under these assumptions, we propose that gapped 
bosonic quantum phases with SG and GG are classified 
by group coholomogy H^{SG x GG,U{1)) in 2-1-1 di- 
mensions, and generally H^~^^{SG x GG,U{1)) in d+l 
dimensions {d > 2). Here "x" is the direct product (or 
the cross product) of two groups, and we will explain the 
notion of group cohomology shortly. 

Let's consider some special limits of our classifica- 
tion. When the system does not have topological order, 
GG = Zi, our classification becomes H'^+^{SG,U{1)). 
This in fact goes back to the group cohomology classifi- 
cation of SPT phaseaiS,. When the system does not have 
global symmetry, SG = Zi, our classification becomes 
H'^^^{GG, U{1)). In 2-1-1 dimensions, this coincides with 
the Dijkgraaf-Witten classification^' of topological quan- 
tum field theories with discrete gauge groups. 

When both the SG and the GG are non-trivial, we will 
show that the indices of the classification H'^^^{SG x 
GG, U{1)) can be expanded as: 

H'^+^{SG X GG, U{1)) = H'^+^{SG, U{1)) 

X H'^+^{GG, U{1)) X SET{SG, GG), (1) 

where SET{SG,GG) will be introduced later. 



SET{SG, GG) describes the non-trivial interplay 
between the topological order and the global symmetry, 
and classifies the symmetry enriched topological phases. 

Some detectable signatures of SET phases, for exam- 
ple, the symmetry protected degeneracy of excited states, 
are studied in this paper. We leave the general numeri- 
cal/experimental signatures of SET phases as a subject 
for future investigation. Nevertheless, we provide ex- 
actly solvable local bosonic models for every phase in 
our classification, in which the topological order is emer- 
gent. These models would be useful tools to further study 
the properties of these phases, including detectable sig- 
natures. 

The plan of this paper is as follows. In Sec|ITl we 
provide the mathematical background of our classifica- 
tion. We review a previously known partial classification 
of SET phases: the projective symmetry group (PSG), 
and comment on the general notion of "symmetry frac- 
tionalization" . In particular, we show that our classi- 
fication includes the mathematical structure underlying 
the PSG classification in 2-1-1 dimensions, and goes be- 
yond it. Namely, our classification contains phases that 
are not described by the PSG. In 3-1-1 dimensions, our 
classification becomes very different from the PSG clas- 
sification and we will explain the reason in Sec lVIl In 
Sec lIIIl we focus on 2-1-1 dimensions and present the ge- 
ometric interpretation of group cohomology, leading to a 
class of exactly solvable models. Each model corresponds 
to a phase in our classification. Generalizations to higher 
dimensions will be briefly discussed. Staying in 2-|-l di- 
mensions, in Sec. lIVI we study the elementary excitations 
of these models, namely gauge fluxes and charges, by 
introducing string-like operators. In Sec|Vl we will solve 
these models in 2-1-1 dimensions for some illuminating ex- 
amples. One particularly important case is the simplest 
example that is NOT described by the PSG classification 
nor "symmetry fractionalization" . In that example the 
global symmetry transformation interchanges the quasi- 
particle species. Detectable signatures of these examples 
will be studied. In Sec lVIl we consider generalizations of 
our study, comment on relations with previous work, and 
conclude. 



II. THE CLASSIFICATION 
A. Mathematical preparation 

1. Definition of the cotiomology group 

We begin with a brief introduction to group cohomol- 
ogy. A detailed introduction can be found in Ref.^^; in 
this paper, we will not present the most general definition 
of group cohomology. 

For a finite group G, and an abelian group M (M 
does not need to be finite or discrete), one can con- 
sider an arbitrary function that maps n elements of G 
to an element in Af; uj : G" -^Mot equivalently 
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uj{gi,g2, ■■■,gn) e M, Vgi,52, •■•ffn e G. Such a group 
function is called an n-cochain. The set of all n-cochains, 
which is denoted as C'^{G^ M), forms an abelian group in 
the usual sense: (wi-W2)(5i)ff2, •■•,5n) = ff2, 
^2(31,32, ••■,5ra), in which the identity n-cochain is a 
group function whose value is always the identity in M. 

One can define a mapping S from C'^{G,M) to 
C"+i(G',M): Vw e C"(G,M), define 6lo e C"+i(G, M) 
as 

...,.g„+i) = a;(g2, • a;(~^)"^'(gi, ...,.g„) 

n 

X ■ g^+i,g^+i, ..,g„+i). (2) 

1=1 

It is easy to show that the mapping 6 is nilpotent: 
= 1 (here 1 denotes the identity (n+2)-cochain). In 
addition, for two n-cochains aJi,aj2, obviously S satisifies 
• 0J2) = {6uJi) ■ {Su}2)- 

An n-cochain uj(gi, ...gn) is called an n-cocyle if and 
only if it satisfies the condition: Suj — 1, where 1 is the 
identity element in C'^~^^{G, M). When this condition is 
satisfied, we also say that uj(gi, ...gn) is an n-cocycle of 
group G with coefficients in M. The set of all n-cocycles, 
denoted by Z''{G, M), forms a subgroup of C"(G, M). 

Not all different cocyles are inequivalent. Below we 
define an equivalence relation in Z^{G,M). Because 6 
is nilpotent, for any (n-l)-cochain c{gi, ...,gn-i), we can 
find the n-cocyle Sc. And if an n-cocyle b can be repre- 
sented as 6 = Sc, for some c e G"~^(G, M), b is called 
an n-coboundary. The set of all n-coboundaries, denoted 
by B''{G,M), forms a subgroup of Z"(G,M). Two n- 
cocycles uji,uj2 are equivalent (denoted by lui ^102) if 
and only if they differ by an n-coboundary: coi — 102 ■ 
where 6 e B"(G,M). 

The n-th cohomology group of group G with coeffi- 
cients in M, H"{G,M), is formed by the equivalence 
classes in Z"(B, M). More precisely: H"{G,M) = 
Z''{G,M)/B''{G,M). 

In this paper we will make a lot of use of 3-cocycles uj. 
We will always choose them to be in "canonical" form, 
which means that (51, 52, 33) = 1 if any of 51,92,53 
is equal to 1 (the identity element of group G). For 
any of the inequivalent cocycles mentioned above, it is 
always possible to choose a gauge such that to becomes 
canonicali^i. Specifically, the explicit elementary cocycles 
that we will use in studying examples of our models in 
Section |V] are going to be canonical. 

So far the notions of cocycle and cohomology group 
are quite formal. But it turns out that they have clear 
geometric/topological meanings, which we will describe 
in SecHn 



2. Examples 

H^(G,U(1)) and one-dimensional representa- 
tions of groups: Let's consider the first cohomology 



group of a finite group G with coefficients in t/(l): 
H^{G, U (1)). In this case, the cocycle condition becomes: 

'^(ffi) •'^(32)/t^(5i -92) = 1 (3) 

This means the 1-cocycle u}{g) is a one-dimensional uni- 
tary representation of the group G. And clearly different 
1-cocycles are different representations. A 0-cochain is 
defined to be a constant cq S Consequently, a 

1- cocycle is a 1-coboundary if and only if it is identity: 
w(g) = co/co = 1. We conclude that the H^{G,U{1)) 
is formed by inequivalent one-dimensional unitary repre- 
sentations of G. 

For instance: 

i/i(Z„,[/(l)) = Z„, H\zlU{l))^Zl (4) 

More generally, for any finite abelian group G, due to 
a fundamental theorem, we know that G can be decom- 
posed as G = X .^^3... X Zn\- Because we know the 
one-dimensional representations of all the components, 
clearly, 

H^{G, U{1)) = G, V finite abelian G. (5) 

H^(G, Z) = Zi, V finite G: Following the above dis- 
cussion, H^{G,Z) is formed by the group of all group 
homomorphisms from G to Z. It is straightforward to 
show that the only group homomorphism between a fi- 
nite G and Z is the trivial one. 

H^(G,U(1)) and projective representations of 
groups: The condition for 2-cocycles is: 

^{91,92) •w(5i -52,53) =^^(52, 53) -^(51, 52 -53) (6) 

In fact, 2-cocycle is related to the so-called projective 
representations of groups. In usual unitary group repre- 
sentations, each group element 5 in G is represented by 
a unitary matrix D{g), which satisfies: D{gi) ■ D{g2) — 
D{gi ■ 52). But for projective representations, this re- 
lation can be modified by a phase factor 01(51,52) £ 
U{1): D{gi) ■ D{g2) = 0^(51,52)^^(51 • 52)- And the 
phase factor ^(51,52), which is a function of 51,52, is 
called a factor system. A factor system cannot be ar- 
bitrary. In order to satisfy the associativity condition: 
[D{gi)-D{g2)]-D{g:,) = D{g,) ■ [D{g2) ■ Dig^)], the factor 
system must satisfy Eq.® — the same condition as for 

2- cocycles. 

What is a 2-coboundary? A 2-coboundary 0(51,52) 
can be written as 0(51,52) = c{gi) ■ c{g2)/c{gi ■ 52) for a 
certain 1-cochain c{g). If two 2-cocyles, 01,02 , differ by 
a 2-coboundary: 

f \ f \ ^(51) • c(52) 

^^1(51,52) = ^2(51,52) ^, (7 

c(5i • 52) 

it is obvious that they correspond to equivalent projective 
representations, because one can absorb the 1-cochain 
into D{g) by redefining D{g) = c{g) ■ D{g), after which 
the two factor systems becomes the same (this is actually 
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the definition of equivalent projective representations). 
We conclude that H^{G,U{1)) classifies all inequivalent 
factor systems of projective representations. 

The calculation of H^{G, U{1)) is non-trivial. We list 
some useful results: 

H\z^, u{i)) = z,, H^zl [/(!)) - 

H^{Zn X [/(!)) = Zgcd(n,m)i (8) 

where Zi is the trivial group. 

H3(G,U(1)). The 3-cocycle condition is: 

(^{gi, 92,93) ■uj{g2, 93,94) ■uj{9i,92 -93,94) 
=w(gi -92,93,94) -uj{gi,g2,93 -94), (9) 

and a 3-cocyle a; is a 3-coboundary iff it can be repre- 
sented as 

/ X c{92,93) -c{9i,92-93) 
w 51,52,53) = -7 ^ — 7 7- 10 

These equations may look strange. But after we intro- 
duce a geometric interpretation in Sec lIIIi their meanings 
will become clear. 

We list some useful results for H^{G, U{1)): 

i?3(z„,[/(l)) = Z„, 

H^iZtUil)) = zfe+'=(fc-l)/2!+fc(fc-l)(fc-2)/3!_ (11) 

For instance, H^{Z^,U{1)) ^ Z| and C/(l)) = 

zl 



3. Some useful theorems 

First, it is known that for any finite group G, its every 
n-th cohomology group with n > is a finite Abelian 
group. Below we list a couple of theorems on group co- 
homology that will be used in the following. 

Universal coefficients theorem: This theorem re- 
lates cohomology groups with different coefficients: 



iI"(G, B) = [iJ"(G, Z)(g>B]x Tor(i/"+i(G, Z), B 



(12) 



This formula allows one to compute cohomology groups 
with coefficients in some Abelian group B by using the 
cohomology groups with coefficients in the group of inte- 
gers Z. 

Here "x" is the usual direct product of groups, and we 
need to define the two new operations: "(8>" and "Tor". 
"(g)" stands for the "symmetric tensor product" (over Z) 
between two abelian groups; while "Tor" stands for the 
"torsion product" . 

Instead of explaining the rigorous mathematical defini- 
tions of these products, we simply list some useful results. 



B always equals B ® A (up to isomorphism) , and 

Zn ® Zm = •^gcd(n,m), 
Zn ® Z — Zn, 

Zn®U{\) = Zl, 
Z(®U{1) - {/(I), 
Z'g)Z ^ z, 
(Ax B)(g,C ^ {A(g,C) X {B(g,C) (13) 

The last relation means that (x) is distributive. 

Concerning the torsion product, one also has 
Tor(A,B) = Tor(B, A). In addition, 

Toi-{Zn, Z„i) = Zgi.d(^n,m), 

Tor(Z„,t/(l)) = Z„, 

Tot{A,Z)^Zi, VA, 
Tor(A X B,C)^ Tor(A, G) x Tor(B, G). (14) 

The torsion product is also distributive (the last relation 
above) . 

Using the universal coefficients theorem, one can com- 
pute the cohomology groups with coefficients in Z from 
those with coefficients in U{1), and vice versa: 

i?"(G,C/(l)) 
-[i/"(G, Z) C/(l)] X Tor(F"+i(G, Z), (7(1)) 
=W'+\G,Z), forn>0, (15) 

where we used the fact that H"{G, Z) is a finite abelian 
group for n > so that [H°{G,Z) ® U{1)] = Zi. Note 
that the above equation is invalid if n = 0. At this mo- 
ment, let's define the 0th cohomology group. In this pa- 
per, H°{G,M) = M. Therefore [H°{G,Z) ® U{1)] = 
U{l)^Zi. 

Using Eq. ([T5l) and Eq.©, we have: 



i/2(G, Z) = G, V finite abelian G. 



(16) 



The K linnet h formula: This theorem allows one to 
compute the cohomology group of a direct product of 
groups, using the cohomology groups of its components: 



H'^iAxB, Z) = Y[[H'{A, Z) ® H'^-'iB, Z)\ 

n+1 

X [| Tor(i7*(A, Z), H"+^-'{B, Z)), (17) 



where "fl" is the usual direct product. For example, this 
formula and the following basic results: 



Zn if p is even; 
HP{Zn,Z) = { if p is odd; 

Z ifp = 0, 



(18) 



together with Eq. (|15p , allow one to obtain the results for 
HP{Z^,U{1)) listed previously. 
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B. The notion of symmetry fractionalization and 
the Projective Symmetry Group 

The main goal of this paper is to address the non-trivial 
interplay between global symmetry and topological or- 
der. What is already known about that interplay? One 
phenomenon famously connected to such an interplay is 
the so-called "symmetry fractionalization" . 

Topologically ordered phases feature anyonic quasipar- 
ticle excitations in the bulk. In fact, in some sense these 
quasiparticles are non-local, because one cannot create a 
single quasiparticle excitation in a system with periodic 
boundary conditions (PBC). One must at least create a 
pair: a quasiparticle and its anti-quasiparticle. There- 
fore, a single anyonic quasiparticle state is not in the 
excitation spectrum of a system with PBC. 

The fact that a single anyonic quasiparticle excitation 
is not a physical excitation in a system with PBC has an 
important physical consequence )2£ When the system has 
a certain global symmetry group SG, based on quantum 
mechanics, we know that all excitated states of the quan- 
tum system can be labeled by irreducible representations 
(irreps) of SG. The irreps characterize how the ground 
state and excited states transform under the global sym- 
metry. Then it is natural to imagine that each anyonic 
quasiparticle also has to transform as a certain irrep of 
SG. But this does not need to be true, exactly because a 
single quasiparticle is not a physical excitated state! 

Let's consider a famous example, the = 1/3 Laugh- 
lin state. There the physical system has a U{1) charge 
conservation symmetry, and therefore any physical state 
should be an irrep of this U{1): jV") e^"^^\ip), for 
Ve'^ G U{1). Here the integer m G Z is nothing but 
the total electric charge of the state. But we were also 
told that the anyonic quasiparticle carries 1/3, a non- 
integer, electric charge. A fractional charge is NOT an 
irrep of the global symmetry. 

Clearly, the fractional charge of an anyonic quasipar- 
ticle can be realized exactly because the single anyonic 
quasiparticle is not a physical state. Only when there 
are three (generally multiple of three) quasiparticles in 
the bulk can it be a physical state, which carries one 
more electric charge compared to the ground state. This 
phenomenon is called symmetry fractionalization. 

In this example, we can ask a further question: why 
the quasiparticles have to carry 1/3, not 1/5, or some 
other fraction of electric charge? Or, what is the guiding 
principle that dictates this fractional charge? 

One obvious guiding principle is the fusion rule. We 
know that three quasiparticles become an electron af- 
ter fusion, which must carry electric charge one. Conse- 
quently each quasiparticle must carry 1/3 charge. In fact, 
this point of view is conceptually very general. For ex- 
ample, one may even be able to consider topologically 
ordered phases with non-abelian quasiparticles. How- 
ever, the mathematical framework behind this point of 
view, for the most general topologically ordered phases, 
is technically highly nontriviaF^ , and is beyond the scope 



of this paper. 

In this paper, we choose a different point of view, which 
involves a simpler mathematical framework — projective 
representations of symmetry group. The trade-off will be 
that we can only use this point of view to understand 
symmetry fractionalization in certain subclasses of topo- 
logically ordered phasesS. However, this is enough for 
the purpose of this paper. 

The point of view that we choose is the following. Be- 
cause only multiples of three quasiparticles correspond to 
physical states, we can define a so-called Invariant Gauge 
Group (IGG): IGG = Z3 = {1, e^^^^/^^ c^^tt/S}^ ^.^n 

multiply each quasiparticle in the system by a fixed ele- 
ment in IGG. Clearly the total phase becomes unity and 
the physical wavefunction is not modified. 

This IGG tells us that when we implement the global 
C/(l) transformation on each quasiparticle, it is perfectly 
fine to have a phase ambiguity, if and only if (iff) this 
phase ambiguity is an element in IGG — because this am- 
biguity does not modify the physical state at all. There- 
fore, a single quasiparticle does not have to form an irrep 
of SG, but it can form a so-called projective representa- 
tion of SG with coefficients in IGG. Formally, this means 
that a single quasiparticle can transform under the global 
U{1) as: ipqp D{e^^)%p, Ve*^ G U{1), where L»(e'^) 
only needs to be a projective representation: 

£l(e*''0 •^(e'"') ='^(e"'Se"'^)£'(e'(''^+''^^), (19) 

where a;(e*^i, e'^^) G IGG (this is why we say that the 
projective representation has coefficients in IGG). In ad- 
dition, as we learned in Sec lII A2l the associativity con- 
dition is satisfied iff a;(e*«Se*^^) G H'^{U{l),IGG) — a 
2-cocycle of U{1) group with coefficients in IGG. 

In fact, the 1/3 electric charge (in general n/3 frac- 
tional charge with n being an integer) is exactly a pro- 
jective representation of U{1) with coefficients in IGG. 
One can check it explicitly: Let's define the transforma- 
tion law of a single quasiparticle under U{\) as D(e'^) = 
e^e/s^Ve G [0,27r). Clearly w{e^'^\e''^^) = e^^'^/^^ iff 
ei+e2> 27r, V6'i,2 G [0,27r), and uj{e^^\e^^^) = 1 other- 
wise. This w(e*^i',e'^2) is a 2-cocycle G H'^{U{l)JGG), 
because the associativity condition is obviously satisified 
by D{e'%. 

In this example, we learned that a quasiparticle of 
a topologically ordered phase can transform under the 
global symmmetry group SG as a projective representa- 
tion of SG with coefficients in a certain abelian group 
IGG. This point of view is also quite general and is 
enough to characterize symmetry fractionalization in this 
paper. Actually, we will use this point of view to classify 
symmetry fractionalization. Firstly, we comment on our 
choice of notation. 

In this paper, "symmetry fractionalization" is a phrase 
reserved to characterize how the global symmetry is im- 
plemented locally on a single quasiparticle. Here "lo- 
cally" is the key word. It basically means that when we 
claim symmetry fractionalization, we already made a ba- 
sic assumption — that the global symmetry transforma- 
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FIG. 1: Illustration of the basic assumption of symmetry frac- 
tionalization Ea.(|2Up: Under a symmetry transformation U (g) 
with Vg G SG, an excitated state is transformed by the prod- 
uct of local transformation operators Ui{g), with each opera- 
tor only acting on one quasiparticle locally. 

tions can indeed be implemented by local transformations 
of each quasiparticle. 

The Basic Assumption of Symmetry Praction- 
alization: Consider an excitated state of a topologi- 
cally ordered phase with a global symmetry group SG, 
having n-quasiparticles (which do not have to be of the 
same type) spatially located at positions ri,r2...r„, far 
apart from one another. Let's denote this state by 
\ip(ri,T2, ...r„)). For any symmetry transformation U{g) 
by a group element g e SG, clearly U{g) will gen- 
erally transform this state into another state: U{g) : 
|'0(ri,r2, ...r„)) -i- |Vi(ri, r2, ...r„)). The basic assump- 
tion of symmetry fractionalization is that there exist lo- 
cal operators Ui{g),U2{g), Un{g), such that Ui{g) is a 
local operator acting only in a finite region around the 
spatial position r^, and does not touch the other quasi- 
particles; in addition, Ui{g),U2{g), ■■■,Un{g) satisfy: 

Uiig) ■ U2ig) ■ ■ ■Unigmvi,r2, ...vn)) 
=C/(g)|^(ri,r2,...r„)) = |^(ri, r2, .., r„)) (20) 

Pictorially, this assumption is shown in FigHJ 

Wen^^ first attempted to classify symmetry fraction- 
alization while investigating the parton mean-field states 
of quantum spin liquids in presence of global symmetries 
such as lattice space group symmetries, time-reversal 
symmetry, and spin-rotation symmetry. In Refls^. Wen 
introduced the notion of Projective Symmetry Group 
(PSG), a very useful tool to classify parton states, as 
well as the low energy gauge fiuctuations, in the pres- 
ence of these global symmetries. We will not provide 
a detailed review the PSG classification here, however, 
we introduce the core mathematical structure underly- 
ing the PSG, and comment on the connection between 
PSG and the current work, in Sec. Ill CI Also, a brief 
introduction to the parton construction and PSG can be 
found in Appendix \X\ 

One way to understand PSG is via the low energy effec- 
tive theory of a state with topological order. Let's con- 



sider the following situation: there is a quantum state 
whose topological order is described by a gauge group 
GG. Therefore, there are gauge charge excitations of the 
gauge group GG in the system, which are anyonic quasi- 
particles. In order to write down an effective theory in 
terms of these gauge charge quasiparticles, it is crucial to 
understand how they transform under the global symme- 
try SG. Similarly to the usual Ginzburg-Landau theory, 
any local term that is allowed by global symmetry will 
appear in the effective theory. Naively, one may expect 
that the gauge charges must form representations of the 
symmetry group SG. 

However, Wen pointed ouiM that these gauge charges 
do not have to form representations of the SG. Instead, 
gauge charges transform under a larger group: the PSG. 
The relation between PSG, SG and GG is given by: 

PSG/GG = SG. (21) 

Mathematically, PSG is a group extension of SG by the 
group GG. From here on to the end of this section, we as- 
sume for simplicity that GG is a finite abelian group. In 
this case, PSG can be shown to be the central extension 
(i.e., GG is in the center of PSG) of SG by GG (see Ap- 
pendix . The classification of how gauge charge quasi- 
particles transform under the SG becomes the classifica- 
tion of all inequivalent central extensions of the group. 
There is a nice mathematical theorem (see, for example, 
Ref§2) stating that all inequivalent central extensions of 
the group SG by GG are classified by H^{SG, GG). 

The mathematical structure H^{SG,GG) underlying 
the PSG classification, which has been independently ob- 
served by several people^^, is somewhat mysterious at 
this moment. But in fact its physical meaning can be 
easily understood. To proceed, again for simplicity, let's 
assume GG has the form GG = ZnX Zm, and one can eas- 
ily generalize the following discussion to any finite abelian 
gauge group. In this case, in order to understand the 
symmetry transformations of the gauge charges, we only 
need to consider two fundamental gauge charge excita- 
tions: ipi and which carry gauge charge (1,0) and 
(0,1), respectively. (Note that we adopt the notation 
(a, b) to label gauge charge here, 0^a<n,0^6< m.) 
This is sufficient because one can build any gauge charge 
quasiparticle by fusing tpi and V'2- 

What are the most general possible ways in which ipi 
and '02 can transform under SG? This is a big question 
and we will attempt to provide an answer later in this 
paper. In this section, however, let's consider a smaller 
question: Under the assumption of symmetry fractional- 
ization, what are the most general possible ways in which 
"01 and "02 transform under SG? 

Under this assumption, symmetry transformations of 
quasiparticles are realized by local operators, which can- 
not change the quasiparticle's species (or more precisely, 
the superselection sector of a quasiparticle). Therefore, 
the gauge charge will be invariant under SG transforma- 
tion: ipi only transforms into ipi while ip2 only transforms 
into "02 ■ However, similarly to the situation with frac- 
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tional charge in fractional quantum Hall states discussed 
above, V'l (or '02) does not need to form a representa- 
tion of SG. This is because any excited states with PBC 
must contain a multiple of n number of Z„ gauge charges 
ipi, and a multiple of m number of Zm gauge charges 
ip2- Consequently, when we define symmetry transfor- 
mations of tpi ('02), it is perfectly fine to have a phase 
e*2i--J- ^gj27r-^^ ambiguity. The state ipi (■02) only needs 
to form a projective representation of SG with coeffi- 
cients in the Zn {Zm) subgroup of U{1), which is ex- 
actly classified by H^{SG,Zn) {H^{SG,Zm)). Finally, 
because we can pair up any two transformation laws 
of '01 E^nd V'2i the symmetry transformations of gauge 
charge quasiparticles with GG = Zn x Z,„ are classified 
hy H^SG,Zn)xH^{SG,Zm). 

Based on the universal coefficients theorem, it is 
straightforward to show that: 

H^{SG, Zn X Z„,) = H^iSG, Zn) x H''{SG, Z„,), (22) 

which is exactly the mathematical structure H^{SG, GG) 
underlying the PSG classification. 

Through this example, we learned that H^{SG, GG) is 
a classification of different ways in which anyonic quasi- 
particles transform under the global symmetry SG, under 
the assumption of symmetry fractionalization. There- 
fore, in this paper we will refer to H^{SG,GG) as the 
symmetry fractionalization classification, and the classes 
contained in H^{SG, GG) as the symmetry fractionaliza- 
tion classes. 

There is one important point that we have not men- 
tioned. We have shown that H^{SG,GG) classifies how 
gauge charges transform under SG. But we also know 
that there are other quasiparticle excitations in the sys- 
tem, such as gauge fiux excitations. For instance, in a 
GG = Z2 topologically ordered phase, there are three 
species of non-trivial quasiparticles: Z2 gauge charge e, 
Z2 gauge fiux m, and their bound state em. In a usual 
topologically ordered state described by an abelian gauge 
group GG, the gauge charges and gauge fluxes are dual 
to each other in 2-|-l dimensions. For instance, it does 
not matter if one labels e or m as the gauge charge in 
the usual Z2 gauge theory. 

In fact, the above discussion indicates that 
H'^{SG,GG) is only a classification of symmetry 
fractionalization for GG gauge charges (or gauge fluxes) 
only, but not for both gauge charges and gauge fluxes. 
That means that the full classification of symmetry frac- 
tionalization should go beyond H^{SG,GG). However, 
we will see shortly in Sec lII CI that our classiflcation 
of symmetry enriched topological phases only contains 
H'^{SG,GG). In addition, in our exactly solvable 
models, we will show that this H^{SG,GG) only 
corresponds to the symmetry fractionalization of the 
gauge fluxes. It turns out that in these exactly solvable 
models, the gauge charges always have trivial symmetry 
fractionalization.'''* We will comment on this issue in 
SeclHCl and in SeclVll 



Now let's consider some simple examples to see the 
power of H^{SG,GG). For the reason mentioned in 
the previous paragraph, in the following examples we 
describe H^{SG,GG) as the symmetry fractionalization 
classifcation of the gauge fluxes. 

1. GG = Z2, and SG = Z2. 

Let's denote the generator of SG as cr, and denote 
by Dm{<j) the transformation of the Z2 gauge flux 
m by cr. Because SG = Z2, we have: 

CT^ = 1. (23) 

The universal coefficients theorem allows us to com- 
pute: 

H^{SG, GG) = H\Z2, Z2) = Z2. (24) 

It means that there are two symmetry fractionaliza- 
tion classes of the Z2 gauge flux. They correspond 
to: 

A„(a)2 = ±1. (25) 

These two possible signs are exactly the two in- 
equivalent cocycles in H^{Z2, Z2). The positive 
sign is the trivial symmetry fractionalization class, 
while the negative sign is the non-trivial class. 

2. GG = Zl and SG = Z2 

Let's denote the generator of SG by a. Now there 
are two fundamental Z2 gauge fluxes: mi, the tt- 
flux in the flrst Z2 gauge group, and m2, the tt- 
flux in the second Z2 gauge group. Straightforward 
computation gives: 

H^{SG, GG) = H'^{Z2, Zl) = Zl. (26) 

There are 4 classes. The corresponding transfor- 
mations of the Z2 gauge fluxes mi, m2, denoted by 
Ani(a'), An2(o'), satisfy: 

DmA^f - ±1, Dm^of = ±1. (27) 

3. GG = Z2, and SG = Zl 

Let's denote the two generators of SG by cr, t. Be- 
cause SG = zl, we have: 

CT^ = 1, = 1, err = Tcr. (28) 

Straightforward computation gives: 

H^{SG, GG) = H^{Zl, Z2) = Zl. (29) 

There are 8 classes. The corresponding transfor- 
mations of the Z2 gauge flux m, denoted by Dmicr) 
and D„i{t), satisfy: 

D„,{<j)^ = ±1, A«(t)' = ±1, 
Dm{a)Dm{T) = ±Dm{T)Dm{a). (30) 
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C. The classification and connection to previous 
work 

Quite some time ago, Dijkgraaf and Witten pointed 
out that the topological orders in 2+1 dimensions, 
described by discrete gauge theories with a gauge 
group GG are classified by its third cohomology group: 
H^{GG,U{1)).^^ Different topological orders labeled by 
H^{GG, U{1)) can be viewed as different discrete versions 
of the Chern-Simons terms^"—-. For example, because 
H^{Z2,U{1)) — Z2, there are two distinct topological or- 
ders described by a Z2 gauge group. In the language of 
the if-matrix, the two topological orders are described by 

2 and K ^ ^ , respectively. The first 

one is the usual Z2 gauge theory while the second one 
is the so-called double-semion theory. The quasiparticle 
anyonic statistics in the two theories are different. 

Recently, an original work by Chen et al*i2, showed 
that bosonic SPT phases protected by a global (unitary) 
on-site symmetry group SG in 2-1-1 dimensions are also 
classified by H^{SG, U{1)). Here different phases labeled 
by H^{SG,U{1)) can be viewed as different topological 
0-terms on a discrete space-time. For instance, because 
H^{Z2,U{1)) — Z2, there are two distinct Ising param- 
agnetic (namely disordered) phases (without topological 
order) in 2-1-1 dimensions. One is the usual Ising para- 
magnet, while the other one is the non-trivial Ising SPT 
phase which features symmetry protected gapless edge 
states. 

It appears that the mathematical object H^{G,U{1)) 
shows up in these two completely different physical con- 
texts, and one may wonder if there is a certain underly- 
ing relation between them. A recent beautiful work by 
Levin and Gu^^ demonstrated such an underlying rela- 
tion explicitly. It was known that the deconfined phase 
of a usual Z2 gauge theory is dual to the usual Ising 
paramagnetic phase. What was shown in Ref.27 is that 
following the same duality, using exactly solvable models, 
the double-semion gauge theory is dual to the nontrivial 
SPT phase. And it was proposed that such dualities be- 
tween the Dijkgraaf- Witten theories and the SPT phases 
are general^^. 

The observation made by Levin and Gu is illuminat- 
ing and motivated us to consider the cases where both 
a global on-site symmetry group SG and a topological 
order described by a gauge group GG are present. Let's 
consider such a gapped quantum phase. On one hand, 
one can imagine following the route of duality transfor- 
mation to transform SG into a gauge group, and even- 
tually having a quantum phase with topological order 
described by the gauge group SG x GG. One the other 
hand, one can follow the backward duality transforma- 
tion to transform GG into a global on-site symmetry, 
which eventually gives a quantum phase with a global 
on-site symmety SG x GG. If the initial phases are dis- 
tinct, it is natural to expect that the phases after duality 
are also distinct, and vice versa. 



global symmetry SG only acts here 




topological order GG only lives here 

FIG. 2; A bilayer system in which the topological order and 
the global symmetry decouple. 



Therefore it is reasonable to expect that, in 2+1 di- 
mensions, bosonic phases with both a global on-site sym- 
metry group SG and a topological order described by a 
gauge group GG are classified by H^{SG x GG,U{1)). 
We will construct exactly solvable models for these 
phases shortly, and we will solve these models in some 
examples and discuss the measurable differences between 
different phases. 

Intuitively, a classification of phases having both SG 
and GG should at least include H^{SG,U{1)) and 
H^{GG,U{1)). This is because one can always consider 
a system where the degrees of freedom which give rise to 
the topological order and the degrees of freedom on which 
the symmetry group acts completely decouple from each 
other. For instance, we can consider a bilayer system, in 
which the global symmetry SG only acts non-trivially on 
the first layer, while the topological order described by 
GG only lives on the second layer (see Figl2]). In this case, 
the possible phases living on the first (second) layer would 
be classified by [/(!)) (H^ {GG , U (l))). Because 

one can tune these phases separately, a classification of 
phases with both SG and GG should actually at least in- 
clude the cross product: H^{SG, U{1)) x H^{GG, U{1)). 
These indices are labeling the phases with no interplay 
between the global symmetry and the topological order. 

This intuitive argument also indicates that, if a clas- 
sification contains more phases than H^{SG,U{1)) x 
H^{GG,U{1)), the extra phases must have non-trivial 
interplay between the global symmetry and the topolog- 
ical order. 

Using the Kiinneth formula and the universal coeffi- 
cients theorem, we can immediately examine whether or 
not our classification is consistent with the above physical 
intuition: 

H^{SG X GG, U{1)) ^ H^{SG x GG, Z) 

=H'^{SG,Z) X H^{GG,Z) x SET{SG,GG) 

^H^{SG,U{1)) X H^{GG,U{1)) x SET{SG,GG). 

(31) 

Note that to obtain the first two terms, we have used 
H°{G, Z) = Z and K iS> Z = K, y abelian finite group 
K. Here we define the abelian group SET{SG, GG) as 
all the other terms in the Kiinneth expansion formula. 
For reasons that will become clear shortly, we further 
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decompose SET{SG, GG) into two parts: 

SET{SG,GG) = SFC{SG,GG) x EXTRA{SG,GG), 

(32) 

where 

SFC{SG, GG) = [H^{SG, Z) H^{GG, Z)] 

X ToiiH^iSG, Z), H^{GG, Z)), (33) 

and 

EXTRA{SG, GG) = Toi{H^{SG, Z),H^{GG, Z)), 

(34) 

where we used the fact that H^{G, Z) — Zi, V finite G. 
In Ea.(|5T]). we see that indeed our classification contains 
H^{SG,Uil)) X H^{GG,U{1)), which is what one ex- 
pects. When we choose the indices in SET{SG,GG) 
to be trivial, i.e., the identity group element in 
SET{SG, GG), these terms label the phases in which the 
topological order and the global symmetry are decoupled. 
Clearly, the indices in SET{SG, GG) are characterizing 
the non-trivial interplay between the topological order 
and global symmetry; namely global symmetry and topo- 
logical order together enrich the classification. The no- 
tation SET{SG, GG) follows from "symmetry enriched 
topological order" . 

The potential physical meaning of SFG{SG, GG) be- 
comes clear if GG is an abelian group. In this case we can 
consider the symmetry fractionalization classes, which 
are given by H^{SG,GG) as discussed in Sec lIIBI Us- 
ing the universal coefficients theorem: 

H^(SG,GG) 
= [H^{SG, Z) (g) GG] X Tor{H^{SG, Z), GG) 
=SFC{SG, GG), if GG is abelian, (35) 

where we used the fact that H^{GG,Z) = 
H^{GG,U{1)) = GG, if GG is a finite abelian 
group. Indeed, in this case, SFG{SG ,GG) has exactly 
the same mathematical structure as the symmetry 
fractionalization classification, leading to the notation 
"SFC" . 

When GG is non-abelian, the projective symmetry 
group is no longer related to the central extensions of the 
SG by GG, and H'^{SG, GG) is not even well-defined. In 
this case, the mathematical structure underlying PSG, 
for symmetry fractionalization classes, was unknown. 
However, SFC{SG,GG) in Eq.® is still well-defined. 
We propose that SFG{SG,GG) is the correct counter- 
part of H^{SG, GG) when GG is non-abelian. 

At this moment, the expansion formula Eg. pip is 
completely mathematical. It appears that the above 
discussion is attaching physical meaning to the terms 
in this formula, such as symmetry fractionalization for 
SFC{SG, GG), without justification. In fact, we will not 
mathematically prove our physical interpretation of the 



formula Eq. ([3T|) generally, although we believe it. How- 
ever, because we have exactly solvable models for every 
phase in the classication H^{SG x GG, U{1)), we can at 
least justify our physical interpretation in some examples 
by solving these models. We will show in Sec|V]that, in 
all the examples that we study, our physical interpreta- 
tion is correct. 

As mentioned earlier, a full classification of symmetry 
fractionalization classes should go beyond H^{SG,GG) 
even when GG is abelian, because one should at least 
consider the symmetry fractionalization classes for both 
gauge charges and gauge fiuxes. However, in the expan- 
sion Eq.dSll), only SFG{SG,GG) appears. We wih show 
that in the exactly solvable models, this SFG{SG, GG) is 
characterizing all the symmetry fractionalization classes 
for gauge fluxes only. It turns out that gauge charges 
in these models always have trivial symmetry fraction- 
alization. Intuitively, this means that our classification 
for the symmetry fractionalization is incomplete. This 
may be due to the fact that we only consider quantum 
phases with exactly solvable model realizations, which 
puts constraints on our classification. 

The extra indices EXTRA{SG,GG) in the expan- 
sion Eq. (j3ip have a completely different mathemati- 
cal structure than symmetry fractionalization classes, 
and intuitively this term must be related to the non- 
trivial interplay between the global symmetry and the 
topological order, but should not be associated with 
symmetry fractionalization. Indeed, we will show that 
EXTRA{SG, GG) is related to the phenomena in which 
global symmetry transformations interchange the quasi- 
particle species (or more precisely, the superselection sec- 
tors). For instance, in the example mentioned in Sec lIIBi 
in which GG = Z„ X Zm, EXTRA{SG,GG) character- 
izes the phenomena where the global symmetry could 
transform a -01 gauge flux into a tp2 gauge flux under 
certain conditions. Such a non-trivial interplay between 
the global symmetry and the topological order is beyond 
symmetry fractionalization, because it violates the basic 
assumption of symmetry fractionalization: it is impos- 
sible to change quasiparticle species by operators acting 
on the quasiparticles only locally. 

Before we move to the exactly solvable models, let's 
present the examples that we will solve in Sec|Vl We 
consider three simple cases: 

1. SG — Z'2, GG ~ Z2- 

H^{SG X GG, C/(l)) = Zl (36) 

and among these: 

H^{SG,U{1)) ^ Z2, H^{GG,U{1)) ^ Z2 

SFC{SG, GG) = Z2, EXTRA{SG, GG) = Z^. (37) 

This means that among indices, one Z2 is la- 
beling the two SPT phases, one Z2 is labeling the 
two Dijkgraaf-Witten topological orders. And the 
remaining Z2 is labeling the symmetry fractional- 
ization classes, whose physical meaning is presented 



11 



in Ea. (|25|) . In this case there is no SET indices be- 
yond the symmetry fractionahzation classification. 

2 . SG ~ Z2 ,GG = Z2 • 

H^{SG X GG, U{1)) = Zl, (38) 

and among them: 

H^{SG,U{1)) = Zl H^{GG,U{1)) = Z2 

SFC{SG, GG) = Zl EXTRA{SG, GG) = Zi. (39) 

This means that among indices, one Zf is la- 
behng the 8 SPT phases, one Z2 is labehng the 
two Dijkgraaf-Witten topological orders. The re- 
maining Z| is labeling the symmetry fractionahza- 
tion classes, whose physical meaning is presented 
in Eg. (|30|) . In this case there is also no SET in- 
dices beyond the symmetry fractionahzation classi- 
fication. 

3. SG — Z2, GG = 

H^{SG X GG, U{1)) = Zl, (40) 

and among them: 

H^{SG,U{1)) = Z2, H^{GG,U{1))= Zl 

SFC{SG, GG) = Zl EXTRA{SG, GG) = Z2. (41) 

This means that among Z\ indices, one Z2 is la- 
beling the two SPT phases, one is labeling 
the 8 Dijkgraaf-Witten topological orders. One Z\ 
is labeling the symmetry fractionahzation classes, 
whose physical meaning is presented in Eq. (P7|) . Fi- 
nally, the remaining Z2 in EXTRA{SG, GG) labels 
the phases beyond the symmetry fractionahzation. 
This is the simplest example in which SET phases 
beyond symmetry fractionahzation are realized. 

III. EXACTLY SOLVABLE MODELS 

In this Section we introduce the exactly solvable mod- 
els which exhibit all the phases from the general classi- 
fication introduced above. First we recall the Dijkgraaf- 
Witten topological invariant, and then introduce the gen- 
eral form of our exactly solvable models. 

A. The geometric interpretation of group 
cohomology and the Dijkgraaf-Witten topological 
invariants 

1. The geometric interpretation of group cohomology 

In Sec lII Al we introduced group cohomology, which ap- 
pears to be a group theoretical concept. However, group 
cohomology is actually about topology. In this section. 
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FIG. 3: The 3-cocycle aj(/ii, /i2, /13) assigns a (7(1) complex 
number (i.e. a phase) to a 3-simplex (tetrahedron), (a) Left 
to Center: the "ordering" of tetrahedron's four vertices, we 
choose here 1 — s- 2 — s- 3 — s- 4. An edge is oriented from lower to 
higher vertex, so no triangle forms an oriented loop. (Alterna- 
tively, one orients all edges without forming oriented triangle 
loops, and a consistent underlying vertex ordering is guaran- 
teed.) Center to Right: "Coloring" assigns group element gtj 
to j ^ i edge, with gji — g~j^ . (Four out of six elements are 
shown explicitly.) The shown tetrahedron 1 — >■ 2 — >■ 3 — >■ 4 
is assigned the phase cj (543, 532, 521)^ = 0^(^43,^32,^21). The 
exponent e = ±1 is determined by: (b) Chirality. For tetra- 
hedron 1 — 5- 2 — s- 3 — )■ 4, looking from vertex 1, (counter- 
)clockwise loop 234 means e — —1 (-1-1), which is real- 
ized in the right (left) tetrahedron, (c) The zero-fiux rule 
applies to all tetrahedron faces, i.e. triangles. Generally, 
9ki ■ gij ■ gjk = 1, the group identity element. Recall that 
9jk ~ 9k j ■ Choosing an ordering and assigning elements to 
ordered bonds, like shown, leads to the constraintz = y ■ x. 

we introduce the geometric interpretation of group co- 
homology, which is the mathematical foundation of our 
exactly solvable models. 

An n-cocyle lu G H"{G,U{1)) of a group G al- 
lows one to construct a topological invariant for n- 
dimensional manifolds. Generally, different elements of 
i/"(G, C/(l)) correspond to different topological invari- 
ants of n-manifolds. Below we will illustrate the con- 
struction of such topological invariants. 

Let's consider a 3-dimensional manifold as an exam- 
ple. We know that tetrahedra can be viewed as building 
blocks for arbitrary 3-manifolds. To begin with, we show 
that a 3-cocycle uj £ H^{G,U{1)) allows one to assign 
a complex number to a tetrahedron following a simple 
procedure. 

The procedure contains two steps (see Fig|31). The first 
step is called ordering, in which one chooses an ordering 
of the 4 vertices of the tetrahedron. We can represent this 
ordering by assigning arrows going from lower to higher 
ordered vertices on the edges of the tetrahedron. For any 
given face (i.e., a triangle) of the tetrahedron, obviously 
the three arrows never form an oriented loop. 

The second step is called coloring, in which one as- 
signs a group element to every edge of the tetrahedron. 
The coloring must be consistent with certain rules be- 
low. Note that an edge already has an arrow, or orien- 



12 




1-4 




FIG. 4: The 1-4 move (3 dimensions) changes triangulation 
but not total product of phases Yii W{ai)'^' of 3-simplices 
aj. Left: The initial tetrahedron is assigned the phase Wo = 
'^(fl43i 532, ff2i)~^ (see Fig. Right: The vertex a is added, 
and we choose the ordering such that 1-4-2— >-3^4^a 
(obvious from chosen orientations of new (red) edges). There 
are now four smaller tetrahedrons, with phases: 1. Tetrahe- 
dron 1 ^> 2 — s> 4 a: Wi = uj^gai, gi2, g2i)', 2. Tetrahedron 
2 — s- 3 — > 4 — 5- a: 14^2 = cj((?a4, (?43, (?32); 3. Tetrahedron 
1 — >• 3 — 5- 4 — s- a: Wz = uj(ga4, g43, g3i)~^ ; 4. Tetrahedron 
1 ^> 2 — >• 3 — >■ a: W4 = cj((7a3, 532, (/2i)~^- The 3-cocycle con- 
dition, Eq. ((9l), says the total phase does not change by the 
move: Wo ~ W1W2W3W4. Note that only one independent 
new group element is introduced (we marked the gai), and 
from zero- flux rule gas = (?a4 • 343. Changing our choice of 
ordering for a relative to 1, 2, 3, 4 would lead to an equivalent 
3-cocycle condition. 
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FIG. 5: The 2-3 move (3 dimensions) changes triangulation 
but not total product of phases Ilj ^(""J^)*^^ • Left: Two 
initial tetrahedrons, 1234 and 1235, are assigned the to- 
tal phase Wo = 1^(543,932,521)1^(543,932,521)"^ (see Fig. [3)). 
Right: One edge (red) is added, and we choose the order- 
ing 4—^5. The volume is now divided into 3 smaller tetra- 
hedrons, with phases: 1. Tetrahedron 1 — >■ 2 — > 4 — >■ 5 
Wi = a;(554,542,52i)~^ 2. Tetrahedron 2 — )► 3 ^> 4 ^ 5 
W2 = 1^(554, 543, 532)"^^; 3. Tetrahedron 1 — >• 3 — 5- 4 — !> 5 
Ws = 0^(554,543,531). The 3-cocycle condition, Eq. says 
the total phase does not change by the move: Wo ~ W1W2W3. 
Note that the new group element 554 is not independent, e.g. 

554 = 552 • 542^ • 



tation, associated with it. The assigned group element 
for a given edge should then be understood in the fol- 
lowing way: if we assign the group element g E G to 
follow the direction of the arrow, then we automatically 
assign group element G G to the direction opposite 
to the arrow. Let's denote the group element assigned to 
the bond connecting vertices j and i as gij , following the 
orientation from j to i: j i. We then automatically 
assign g^i = .g.^\ 

In addition, the three assigned group elements for any 
given face must satisfy the constraint: gij ■ gjk ■ gki = 1, 
where 1 is the identity element in group G, and i,j,k 
are the three vertices of the face. We will call this con- 
straint the "zero-flux rule" throughout this paper. With 
this constraint, it is easy to show that among the six 
group elements for the 6 edges of the tetrahedron, only 
three are independent. In particular, let's denote the or- 
dered vertices by 1,2,3,4; then, 543,(732,521 completely 
determine all the other group elements. 

Given a 3-cocycle uj{x,y,z) E H^{G,U{1)), one as- 
signs the complex number ^(543, 532, 521)'^ to an or- 
dered and colored tetrahedron. (Sometimes we use the 
'-^'^(543)532,521) notation.) Here e — ±1 depending on 
the chirality of the ordered vertices. One can determine 
this chirality by the right-hand rule: imagine looking at 
the face formed by vertices 2-3-4 from the vertex- 1; if the 
vertices 2-3-4 form a counter-clockwise (clockwise) loop, 
the chirality of the ordering is positive (negative) and 

e = l (6 = -l). 

This assignment of ^(543, 532, 521)*^ to an ordered and 



colored tetrahedron allows a simple geometric interpre- 
tation of the cocycle condition Eq.([S]), see FigU) To see 
this, consider an ordered and colored tetrahedron and 
the associated complex number 0-1(543, 532, 521)*^- One can 
now add one more vertex-a inside the tetrahedron. With 
vertex-a, the original tetrahedron can be triangulated 
into 4 smaller tetrahedra. One can further continue the 
ordering and coloring procedure for the 4 smaller tetra- 
hedra. Since we already have the ordering and coloring 
for the large tetrahedron, we only need to assign an or- 
der to vertex-a, as well as to color the four newly created 
edges la, 2a, 3a and 4a. Actually, according to the zero- 
flux rule, it is easy to show that only one of the four new 
edges is independent. After we complete ordering and 
coloring the 4 small tetrahedra, we will have 4 new com- 
plex numbers, each of which is associated with a small 
tetrahedron. It is straightforward to show that, no mat- 
ter how one performs the complete ordering and coloring 
procedure, the cocycle condition Eq. dictates that the 
product of the 4 new complex numbers exactly equals the 
orginal complex number 0^(543, 532, 521)"^. 

Such a procedure of completing the triangulation, or- 
dering and coloring of tetrahedra after adding a vertex 
is called a 1-4 move. A specific example of a 1-4 move is 
shown in FigH) 

Similarly, there is a 2-3 move, Fig. [5j Namely, one can 
consider two face-sharing tetrahedra, both of which have 
been ordered and colored. There are then two complex 
numbers, each of which is associated with a tetrahedron. 
One can now connect the two vertices that are on op- 
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posite sides of the shared face, and the volume enclosed 
by the original two tetrahedra can be triangulated into 
three tetrahedra. One can continue the ordering and col- 
oring procedure for the three tetrahedra and obtain three 
new complex numbers. It is also easy to show that, no 
matter how one performs the further ordering and color- 
ing procedure, the cocycle condition Eq.(^ dictates that 
the product of the three new complex numbers equals the 
product of the two orginal complex numbers. A specific 
example of such a 2-3 move is illustrated in FiglSl 

In this paper we will use "canonical" 3-cocycles oj, 
meaning that a; (31,32,^3) = 1 if any of 91,92,93 is equal 
to 1 (the identity element of group G) . It is always pos- 
sible to choose a gauge for ui such that it becomes canon- 
icali^. Specifically, the explicit elementary cocycles that 
we will use in studying examples of our models in Sec- 
tion |V] are going to be canonical. 



2. The Dijkgraaf- Witten topological invariants 

The above examples of 1-4 move and 2-3 move sug- 
gest that the products of the assigned complex numbers 
^{9ik,gkj,gjiY for a given volume may be related to a 
certain invariant that is independent of the triangulation, 
ordering and coloring procedure. This is indeed true, as 
stated by two mathematical theorems presented in the 
following. 

Let us first consider a closed 3-manifold M without 
a boundary. One can triangulate M by a finite num- 
ber of 3-simplices (i.e., tetrahedra), and then order the 
vertices of this triangulation. Next, one can have a col- 
oring of all the edges in the triangulation obeying the 
zero-flux rule. Note that under a fixed triangulation and 
ordering of vertices, there can be many different color- 
ings. Let's denote a 3-simplex of the triangulation, to- 
gether with the ordering of its vertices, by ct/, where 
I = 1, 2, S* labels 3-simplices and S is the total num- 
ber of 3-simplices. For a given coloring (f, let's also de- 
note the assigned complex number uj{gik, gkj, gjiY for the 
simplex-cT/ as W{ai, ipy^'^'\ and we can further compute 
the product of all these complex numbers for 3-simplices: 
Yii^i^i'^iiVY ■ For each given coloring ip, we will 
have one such product. 

Theorem 1: The sum of such products for all possible 
colorings, with an appropriate normalization factor, is a 
topological invariant of the closed manifold M 

^^'-\i\v E f[W{ai,^)<'^'l (42) 

' ' ip£ all 1=1 
possible 
colorings 

Here IG] is the number of elements in group G, and V 
is the number of vertices in the triangulation. Note that 
without Theorem-1, one would naively expect that Zm 
depends on both the triangulation and the ordering of 
vertices (while different colorings are already summed 



over). But with Theorem-1, we know that Zm does not 
depend on either of them — it only depends on the topol- 
ogy of the manifold M and the 3-cocyle uj e H^{G, U{1)). 
One can further show that equivalent 3-cocycles (i.e., 3- 
cocycles differing by a 3-coboundary) give exactly the 
same topological invariant Zm^ namely, Zm only de- 
pends on inequivalent elements in H^{G, U{1)). 

The topological invariant Zm is exactly the partition 
function of the Dijkgraaf- Witten (DW) topological quan- 
tum field theory (TQFT) for discrete gauge group G in 
2-t-l dimensions^i^. In order to have a well-defined 
TQFT, it turns out that one not only needs to define 
partition functions for closed space-time manifolds, but 
one also needs to define quantum transition amplitudes 
for space-time manifolds with boundaries. This is given 
by the second theorem. 

Consider a 3-manifold M with boundary dM . dM is 
formed by a collection of closed 2-manifolds. One can 
triangulate dM by a finite number of 2-simplices (i.e., 
triangles), order the vertices of the 2-simplices, and then 
color their edges again following the zero-flux rule (i.e., 
9ij9jkgki = 1 for all 2-simplices). Let's denote the trian- 
gulation, ordering, and coloring of the boundary dM by 

T. 

Next, we fix the coloring t and extend it into the bulk 
of M. This means that we consider a triangulation of 
M, an ordering of its vertices, and a coloring if such that 
they become exactly the same as r when limited to the 
boundary dAI. In this case, we also say that the bulk 
triangulation, ordering, and coloring in M are compat- 
ible with T on dM. For instance, the triangular faces 
of a tetrahedron can be viewed as the boundary of a 3- 
dimensional ball. Then a 1-4 move can be viewed as a 
specific extension of the boundary t into the bulk of the 
ball. 

Now let's fix the bulk triangulation and ordering of ver- 
tices in M that is compatible with t. There are still many 
possible colorings tp in M that are compatible with r, 
and they form a set which we denote as Col{M, r). As in 
Theorem-1, with a fixed (p £ Col{M, t) one can compute 
the product of complex numbers n/=i W[ai, ipY'-"'^ as- 
signed to all the 3-simplices in the bulk of M. It turns 
out the sum of all such products satisfies the following 
theorem: 

Theorem 2: The complex number Zm (r) does not de- 
pend on the triangulation of M nor ordering of its ver- 
tices, whenever the topology of M and r on dM are 
fixed^^: 

ZM{r) = E n m<^iM''^''- (43) 

\G\ 2 ipGCol(ALT} 1=1 

Here V is the total number of vertices inside M (i.e., not 
including dM), while VgM is the number of vertices in 
dM. Obviously Zm{t) becomes Zm in Ea. (l42]) when Af 
does not have a boundary. 

To see the physical meaning of Zm{t), let's consider a 
special case: M = B x [0,1] where B is a certain closed 



14 




FIG. 6: The Dijkgraaf-Witten topological invariant Zm and 
TQFT for manifold M with boundary, (a) Two copies of 
2-manifold B: Bi,B2, with colorings ti,T2, form the bound- 
ary of 3- manifold M = B x [0, 1]. The triangulation and 
ordering in Bi and B2 are chosen identical, as sketched. The 
triangulation and ordering in the bulk of M (example red ver- 
tices and edges) does not influence the value of DW invariant 
Zm{t = Ti U r2). A coloring Ti of B represents a quantum 
state In), and then Zm(t) = (T2I P In) is a quantum ampli- 
tude of operator P in the DW TQFT associated with B. The 
image of P is the ground state manifold of the TQFT. (b) The 
operator P is a projector: The quantum amplitude {t2\ P^ |ti) 
is sketched, representing a product of amplitudes of P from 
Ti to T3, and ra to T2, with a sum over colorings T3. This 
becomes a sum over internal colorings in space M = _B x [0, 2] 
with boundary Bi U B2 ■ The latter is equal to the amplitude 
(T2I P |ti) from panel (a). 

orientable 2-manifold. dM is formed by two disconnected 
but identical closed 2-manifolds: Bi = B and B2 = B 
corresponding to S [0, 1] and 1 £ [0, 1] respectively (see 
Fig 16]). We can then triangulate and order the vertices 
on Bi and B2 in the same fashion. 

For each edge-ij of the triangulation of Bi, we con- 
struct a |G|-dimensional local Hilbert space "Hf^^ = 
{gij £ G}; namely each group element labels a quan- 
tum state in Ti-ij- Then we consider the tensor product 
of all such local Hilbert spaces 'H'^^ = edges^^?^' 
Now we can associate each possible coloring on Bi with 
a quantum state in T-L^^ . Because a coloring must sat- 
isfy the zero-flux rule, clearly all possible colorings of Bi 
span a sub-space C . All possible colorings 

of B2 then also form the exact same Hilbert sub-space 

■flDW 

Let's choose a coloring |ti) G on Bi, and an- 

other coloring |t2) S on B2. n and T2 completely 

specify the triangulation, ordering and coloring r on dM. 
Theorem-2 means that there is a well-defined quantum 
transition amplitude from the state |ri) to the state |r2): 

{r2\P\n) = Zm{t). (44) 

Because all possible |ti) (|r2)) form a basis of Ti^^ , this 
equation defines a quantum operator P on . 

Theorem-2 immediately dictates that P is a projec- 
tor: — P. This is because after we insert an iden- 
tity operator Era ^3)^3! = ^ '^''^ ^ ^2|-P'|ri} = 
Et3(^2|P|''3)(t3|P|ti} has a simple geometric interpre- 
tation (see Figin]): we can consider two copies of the 
manifold M, Mi = B x [0,1] and M2 = B x [1,2], so 
that (t3|P|ti) ((r2|P|T3)) is the quantum amplitude due 



to an internal triangulation and ordering of Mi (M2). We 
can then glue the top boundary of Mi with the bottom 
boundary of Af2- After the gluing, the vertices on the 
glued boundary B x {1} become internal vertices. Then 
Et, ('''2|^l'''3)('''3|^l'''i) can be simply interpreted as the 
quantum amplitude due to an internal triangulation and 
ordering of Mi U M2 = M, which must be the same as 
(T2|P|ri} according to Theorem-2. 

Because P is a projector, the image of P forms a sub- 
space in the Hilbert space Img{P) C in which P 
acts as identity. Img{P) turns out to be the ground 
state sector associated with the Dijkgraaf-Witten TQFT 
for the closed 2-manifold B. One can also prove^^^ that 
the dimension of Img{P) (i.e. the ground state degen- 
eracy of the TQFT) and the partition function of the 
closed space-time 3-manifold Af = B x are identical: 
dim[Jmg(P)] = ^m=bxS1 ■ 

At this point, it is useful to introduce an example. 
Consider the simplest group G — Z2. According to 
Eqs. ^ and (HH), H'^{Z2, U{1)) = Z2. This means that 
there are two inequivalent 3-cocycles and let's choose the 
trivial one: uj(x,y,z) = 1, \/x,y,z e Z2, which gives 
rise to a Dijkgraaf-Witten TQFT. This particular TQFT 
turns out to be a familiar one: the Z2 gauge theory of the 
toric code model^^. We can then use Theorem- 1 to com- 
pute the partition function Zm for a closed 3-manifold 
M, and use Theorem-2 to compute the ground state de- 
generacy via the projector P. For instance, for a 3-sphere 
and a 3-torus, Zg3 = 1/2 and Z^^ = 4, respectively. The 
latter result implies that the ground state degeneracy 
on a torus is 4, since — x S^. More generally, 
the ground state degeneracy on a closed orientable 2- 
manifold B is dim[Img{P)] = 4^, where g is the genus of 
B. 



3. The generalization to other dimensions 

The above discussion has been limited to 2+1 dimen- 
sions. In fact, the geometric interpretation of an n- 
cocycle can be easily generalized to any n ^ 2 space-time 
dimensions. Some aspects of this generalization have 
been discussed in ReflTol. Here for the purpose of the 
current paper, we briefly discuss the n = 2 and the n = 4 
cases. 

Geometric interpretation of a 2-cocycle uj{x,y) £ 
i/2(G,C/(l)). (See FigH) Let's choose a 2-cocycle 
uj{x,y) e H'^{G,U{1)). Consider a 2-simplex (i.e., tri- 
angle). Again one needs to perform the ordering and 
coloring precedure. Let's choose an ordering of the ver- 
tices 1—5-2-^3. We then color the edges by group ele- 
ments 331,532,321 under the zero-flux rule: 531 = 532521- 
Therefore, we can choose 532,521 to be the only inde- 
pendent elements. Next, we assign the complex number 
a; ((732, 521)*^ to this 2-simplex. Here e = ±1 depending on 
the chirality of the ordering of vertices: if 1 — > 2 — ^ 3 is 
clockwise (counter-clockwise), e = 1 (e = — 1). 

The geometric interpretation of the 2-cocycle condi- 
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FIG. 7: The 2-cocycle Lj(hi, /12) assigns a ?7(1) complex num- 
ber (i.e. a phase) to a 2-simplex (triangle), (a) Left to Cen- 
ter: the "ordering" of triangle's three vertices, we choose here 
1 — !> 2 — > 3. An edge is oriented from lower to higher ver- 
tex, so no triangle forms an oriented loop. (Alternatively, one 
orients all edges without forming oriented triangle loops, and 
a consistent underlying vertex ordering is guaranteed.) Cen- 
ter to Right: "Coloring" assigns group element Qij to j ^ i 
edge, with gjt = g"^. The shown triangle 1 — > 2 — >■ 3 is 
assigned the phase W = a;(g32,ff2i)^ = 1^(543, 532, 521)- The 
exponent e = ±1 is determined by: (b) Chirality. For triangle 
1 ^ 2 ^ 3, (counter-)clockwise loop 123 means e = -1-1 (— 1), 
which is realized in the left (right) triangle, (c) The zero-flux 
rule applied to a triangle. Generally, Qki ■ gij ■ gjk = 1, the 
group identity element. Recall that gjk — g^^. Choosing an 
ordering and assigning elements to ordered bonds as shown 
leads to the constraint z = y ■ x. 



tion Eq.([B]) can now be understood as the invariance of 
the product of these assigned complex numbers in a 1-3 
move or a 2-2 move (see FiglS]). For instance, in a 1-3 
move, we consider an ordered and colored triangle, to- 
gether with the assigned complex number oj(g32T g2iY ■ 
Then we add one new vertex inside the triangle. After 
connecting the new vertex with the original 3 vertices, 
3 new edges are created and the original triangle is thus 
further triangulated into 3 smaller triangles. We now 
can choose any ordering and coloring of the new vertex 
and new edges under the zero-flux rule, which then as- 
signs 3 new complex numbers for the 3 smaller triangles. 
The 2-cocycle condition Eq.([n]) dictates that the product 
of the 3 new complex numbers equals the original one 
^(332,521)"- 

Theorem- 1 and Theorem-2 can also be generalized to 
2-manifolds and 2-cocycles. For example, let's consider 
Theorem-2. For a 2-manifold M with boundary dM , 
one firstly chooses a triangulation (using 1-simplices, i.e., 
line segments), an ordering of vertices, and a coloring 
on 9M, which we denote by t. Note that now there 
are no zero-flux rule constraints for r, because there is 
no triangle in a 1-simplex. Then one can extend the 
triangulation, ordering and coloring into the bulk of M 
(where the zero- flux rule holds). We denote the assigned 
complex number for a 2-simplex cr/ in M as W{ai, ipY'^'^'^ 
as before, where Lp denotes the bulk coloring. With a 
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FIG. 8: The 1-3 and 2-2 moves (in 2 dimensions) change 
triangulation but not total product of phases Y\i WiaiY' of 
2-simplices cr/. (a) The 1-3 move. Left: The initial triangle 
is assigned the phase Wo = Lu{g32,g2i) (see Fig. [T]!. Right: 
The vertex a is added, and we choose the ordering such that 

1 — >■ 2 — 3 — >■ a (obvious from chosen orientations of new 
(red) edges). There are now three smaller triangles, with 
phases: 1. Triangle 1 — >■ 2 — > a: Wi = uj{ga2, 921)', 2. Triangle 

2 ^ 3 a: W2 = uj{ga3, 932); 3. Triangle 1 ^ 3 ^ a: W3 = 
cj((7a3, <73i)~^. The 2-cocycle condition, Eq. ([6]), says the total 
phase does not change by the move: Wo = W1W2W3. Note 
that only one independent new group element is introduced 
(we marked the gas), and from zero-flux rule ga2 = ffas • <?32. 
Changing our choice of ordering for a relative to 1, 2, 3 would 
lead to an equivalent 2-cocycle condition, (b) The 2-2 move. 
Left: Two initial triangles, 123 and 124, are assigned the 
total phase Wo = [j((?32, (?2i)i^(g42, ff2i)~^. Right: The area is 
divided into two triangles by the other possible edge (red), and 
we choose the ordering 3 — 4. The phases of new triangles: 
1. Triangle 1 ^ 3 ^ 4: Wi = a;(g43, ffai)"^ 2. Triangle 2 

3 — >■ 4: W2 = uj{g4,3, 932). The 2-cocycle condition gives Wo = 
W\W2- Note that the group element 343 is not independent, 
e.g. 9i-3 = ff42 • 932- 



fixed bulk triangulation and ordering, there will be many 
possible colorings that are compatible with r. Theorem-2 
for 2-manifolds and 2-cocyles states that Eq. ((43)) defines 
a complex number Zm{t) which is independent of the 
choice of bulk triangulation and ordering of vertices, as 
long as T is fixed. 

Following the discussion from the previous sec- 
tion, Theorem-1 and Theorem-2 suggest that a cocyle 
uj{x, y) G H'^{G, U{1)) may define a 2D TQFT. This is in- 
deed true and was discussed in a mathematical context^^. 
In 2-|-ld, we know that different topological orders can 
be characterized by different TQFTs. One may ask: does 
this mean that there are non-trivial topological orders in 
1+ld? However, we also know from previous research 
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that non-trivial topological orders do not exist in l+ld^. 
It turns out that the 2D TQFTs induced by 2-cocycles, 
do not give rise to physically non-trivial topological or- 
der. This is because the ground state degeneracy is not 
robust, as one can lift it by a local perturbationj22 

Geometric interpretation of a ^-cocycle uj(x,y, z,u) S 
H'^{G,U{1)) . Similarly to the n = 3 case, for a given 
4-siniplex one can choose an ordering of its vertices 1 — 
2 ^> 3 ^> 4 -> 5, color the edges following the zero- flux 
rule, and assign the complex number uj^g^i, 343, 5321 521)*^ 
to it. Here again e is determined by the chirality of the 
ordering of vertices. The 4-cocycle condition duj = 1 in 
Eq.(l2]) when n = 4 can be understood as the invariance 
of the product of the assigned complex numbers to 1-5, 

2- 4, and 3-3 moves. Theorem-1 and Theorem-2 also hold 
for 4-manifolds. For any fixed 4-cocycle uj{x,y, z,u) e 
H'^{G,U{1)), these theorems give rise to a 4D TQFT. 
Equivalent 4-cocycles induce the same TQFT. These 4D 
TQFTs characterize different topological orders in 3-|-ld. 

B. Exactly solvable models 

In this section we define our exactly solvable models. 
Although we discuss the generalization to other dimen- 
sions in Sec. I VI Al from now on we constrain ourselves to 
the 2-1-1 dimensional case. It will become clear that our 
models exhibit both a global symmetry forming a group 
SG, as well as topological order described by a discrete 
gauge group GG. We will explain our models' relation to 
both the SPT models of Refill and the Dijkgraaf-Witten 
gauge theories of ReffsTI. Through these connections it 
will also become clear that each inequivalent choice of 

3- cocycle in our models leads to a model with specific 
topological and symmetry properties, as labeled by the 
classification in Sec. HIl 

We consider a triangular (two dimensional) lattice with 
oriented edges (bonds), where these orientations are com- 
patible with an ordering of lattice sites, i.e. each edge 
oriented from a lower to higher ordered site and no tri- 
angle edges form an oriented loop, just as we discussed 
in Section UlI Al and Figure [T] For concreteness, in Fig.[S] 
we show our choice of edge orientations on the triangular 
lattice. We next introduce the "coloring" (p by assigning 
an element hij £ GG to each oriented edge j ^ i, again 
as discussed in Section FlII Al however, we now also assign 
a group element ui G SG to each lattice site i. 

Actually, we further introduce the group element 

9ij = hij ■ Ui ■ uj^ e G (45) 

as our variable on edge ij. This definition might seem 
somewhat redundant, since the Ui elements appear both 
on sites, and on edges through gij. It however has im- 
portant meaning. As discussed in Section III C[ it is al- 
ready known that some cohomology based classifications 
of phases with symmetry and phases with topological or- 
der can be explicitly connected by duality. Due to the 
direct product structure of the group G = SG x GG we 



consider here, it is simple to dualize either SG or GG (en- 
tire groups or their subgroups) without having to change 
the formalism. We will in fact use dualization of SG ex- 
plicitly when considering symmetry protected degeneracy 
in examples. Section IV Gl 

Let us then briefly consider how SG is dualized. The 
group variables Ui defined on lattice sites i are replaced 
by group variables utj living on oriented edges j ^ i 
according to the rule2^ 

Uij =Ui- uj^. (46) 

Due to importance of duality, we want to ensure that 
all gauge degrees of freedom present in a theory are 
treated equally. This is our motivation for using the edge 
variables gij defined in Eq. (|45p as degrees of freedom on 
equal footing with hij . 

An arbitrary quantum state in the Hilbert space Ti of 
our model is therefore labeled by |i) ~ \{ui} , {gij}) , or 
by |i) = |{uj,{/iy}}. 

The elementary building block for the theory is the 
operator labeled by a group element s G G = SG x 
GG, and a plaquette p containing six triangles sharing 
a lattice site i at the center. The plaquette operator 
therefore acts on seven group elements, one at the site i 
and six on the edges that share this site. To define its 
action, we need an additional edge oriented vertically up 
into the third dimension at site i, to which we assign the 
element s G G which can always be uniquely factored as 

s = hs-s, (47) 

with hs e GG, s e SG (Fig.©. The operator trans- 
forms the seven values of elements in the plaquette by s, 
preserving the orientation of edges, and these new values 
are represented on edges lifted above the original ones, 
see Fig. IHl With p centered on site i, we have: 

s • u, (48) 
hij y hg • hij 
hki hki ■ h~^, 

leading to 

gij s ■ gij (49) 

gki gki ■ s~^- 

Further, non-zero matrix elements of B^, 

B;^{iis)\B;\i), (50) 

are assigned the following quantum amplitude 

6 

B;^Y[W{aj,^r^'^'\ (51) 

7=1 

where the six 3-simplices aj are built using the six tri- 
angles of the plaquette p (with the initial group element 




FIG. 9: Action of the operator Bp, s = hs ■ s, with hs G GG, 
s G SG, on plaquette p centered on site i: For six el- 
ements hii £ GG on edges il it resembles a local gauge 
transformation preserving zero-flux through all triangles, e.g. 
h[j = hs ■ hij and h'j^i — hki ■ hj^ , while on-site element 
transforms as u'^ — sui, leading to e.g. g'ij = s ■ Qij and 
g'ki ~ Qki ■ s~^\ additionally, there is an overall phase factor 
which is the product of complex numbers assigned by the 3- 
cocycle ijj £ Z^{SG x GG, (7(1)) to the six 3-simplices (tetra- 
hedrons) forming the 'pyramid' (see Fig. and Eq. (jSip ). 
Note that fixing the initial and final state (i.e. the values 
Ui,u'i, etc.) leads to a unique value of s for which the opera- 
tor matrix element does not vanish. In that sense, this picture 
also describes the action of the plaquette operator _Bp. 

values assigned), the vertical edge (assigned the group 
element s), and the six lifted edges (assigned the final 
element values according to Eq. (|49|). This action is 
shown in Fig. [HI The orientation of new (lifted) edges is 
chosen to match the orientation of original edges upon 
downward projection. 

It is important to note that the zero-flux rule (dis- 
cussed in Sec. IIII Ap is by construction satisfied on all 
faces (triangles) of the six tetrahedra, if it is satisfied 
in the 6 triangles of the plaquette p. The zero-fiux rule 
must hold on all faces of the tetrahedra for which we are 
calculating the phase W. The operator is therefore 
defined only in a Hilbert subspace /Cp which consists of 
states having the zero-flux rule satisfled in all six trian- 
gles of the plaquette p. Finally, note that choosing a final 
state |f) in a non-zero matrix element fixes a unique value 
of s. 

We finally define the plaquette operators Bp as having 
matrix elements 

To explicitly illustrate the plaquette operator (as well 
as the full Hamiltonian defined below) through examples, 
in Appendisin] we will consider the models for two Z2 
topologically ordered phases, i.e. the well-known "toric 
code"— and "double semion" theory^. 

Returning to the most general case, our plaquette op- 
erators Bp turn out to be projectors. Namely, using the 
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properties of 3-cocycles, one finds that applying a 
operator twice at the same plaquette leads to a group 
multiplication in the amplitude, 

(f|B^B;'|i)=i?-^'. (53) 

Using the normalization in Eq. ((5^ it follows that 

(f| BpBp |i) = Bp, (54) 

and also that Bp is a projector. 

Crucially, we will show further below that the plaque- 
tte operators commute: 

[Bp,Bp,]=Q, \/p,p'. (55) 

Let us next introduce the operator Qt, which projects 
flux in a triangle t to zero, i.e. it enforces the zero-fiux 
rule discussed in Section ITlIAl In other words, Qt is non- 
zero (and equal to 1) only when acting on a triangle t 
made out of lattice sites i,j,k such that 

hij ■ hjk ■ hki = Igg (56) 

where 1 is the group identity in G, and the second line 
follows directly from the definition Eq. (|45|) . 
We can at last define the Hamiltonian as 

^=-EQ*-E^pnQ*' (57) 

where the label t ^ p enumerates the six triangles mak- 
ing up the plaquette p. As mentioned above, the fac- 
tor HtGp '5* is actually crucial to ensure that H is well- 
defined: it ensures that Bp acts within the subspace }Zp 
on which it is defined (see discussion after Eq. (ISTI) '). 

Further, it is easy to see that plaquette operator term 
Bp Htep actually commutes with the projectors Qf. 
Namely, the transformation rule by s in operator B'L, as 
introduced above, preserves the product rule Eq. (1561) on 
all triangle faces of simplices in Fig. [9l if it is satisfied in 
either the upper or lower triangles, i.e. either in the |i) 
or |f) state. Obviously then the zero-flux rule enforced 
by action of Qt' commutes with the action of Bp Htep Qt 
even when t' belongs to the plaquette p. 

Our model has the global symmetry group SG, follow- 
ing from the fact that the Hamiltonian commutes with 
the global symmetry operations 

Ui ^ Ui ■ s~^ , V«, (58) 

and s G SG. The symmetry operation obviously does 
not influence the zero-flux rule in Eq. (j56p . and therefore 
commutes with every Qt- Considering a plaquette oper- 
ator, the symmetry operation leaves the edge elements 
gij invariant, and also the final value of site elements Ui 
is the same no matter the order in which apply Bp and 
s, due to the group property. 
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Next, our model H is exactly solvable: All terms in the 
Hamiltonian H commute with each other (we still have 
to prove the commutation of Bp, Bpi), so the model is 
exactly solvable. Let us now consider the ground state 
manifold of our model. Since all the terms in H are also 
projectors, the ground state manifold is the image of the 
projector P = Y\p^pY\t£pQt- Actually, it is also easy 
to see that BpBp = Bp due to Eq. ([55t and the group 
property. This means that for a ground state it also holds 
that = l,Vp,s. 

First, let's consider the special cases in which SG — Z\ 
is trivial. In this case, the projector P = J^^ Bp Htep Qt 
is exactly the projector in the Dijkgraaf-Witten theory, 
Eq. Namely, applying all Bp operators in the plane 

creates a lifted copy of the plane, leaving the volume 
between them triangulated by tetrahedrons; the transi- 
tion amplitude for this operation is equal to the prod- 
uct of H/^i ^{pit'vY^"^^ phases contributed by aU the 
J|pi?p. The initial and final states fix the coloring r 
on the two planes, so the transition amplitude exactly 
equals the Dijkgraaf-Witten topological invariant Zm{t) 
(Eq. evaluated on the manifold having the two 

planes as boundaries. (Note that there are no vertices in- 
side the volume, and the number of plaquettes p is equal 
to VaA//2 since there are two planes in dM, leading to 
correct prefactor from Eq. (|43p .) We can therefore con- 
clude that when SG = Z\, the ground state sector of 
our model, to which P projects with eigenvalue 1, is also 
the ground state sector of the Dijkgraaf-Witten TQFT-Si 
defined on the triangular plane. We will generally study 
the topological order of our models in Appendix |B1 

On the other hand, we can consider the opposite situ- 
ation where the gauge group is trivial GG ^ Z\, so that 
G = SG. In that case, our models become equal to the 
exactly solvable models for symmetry protected topolog- 
ical phases constructed by Chen et al. in Section IIP of 
Ref,19,. Namely, since the only degrees of freedom on the 
edges are from SG, see Eq. (gS]) , the zero- flux rule is au- 
tomatic so Qt = 1. The Hamiltonian is just a sum of the 
Bp plaquette operators, and these are obviously identical 
to the plaquette operators forming the Hamiltonian in 
Reffiol. We can conclude, as claimed in the introduction, 
that our exactly solvable models in the case of trivial 
gauge group GG become the models for symmetry pro- 
tected topological phases classified by H^{SG, U{1)). 

Our exactly solvable model is in-between these two ex- 
treme situations (the DW theory and the SPT model), 
and it can be understood as a partially dualized version 
of either of them. 

In general, the ground states of our models do not 
break the physical symmetry SG so that they describe 
symmetric quantum phases. The simplest way to con- 
vince oneself of this is by noting that the Bp operators 
in our models (see Eq. ([50l) ) create/annihilate small do- 
main walls in a quantum state when s G SG. Since in a 
ground state it holds that 13^ = 1, Vp, the ground state 
is a domain wall condensate — i.e., the symmetric phase. 



Let us now prove Eq. ((55)) by using the Dijkgraaf- 
Witten topological invariant from Eq. P5)) . Consider the 
picture of action of two overlapping plaquette operators, 
described by matrix elements 

BBi = (f| BpBp, |i) , (59) 

and 

BB2 = {i\Bp,Bp\i), (60) 

where p and p', the plaquettes centered on sites i,j re- 
spectively, necessarily share two triangles, while only the 
ij edge is acted ugon by both operators. Fig. [TOl The 
operator product BpI3pi is obviously defined within the 
Hilbert subspace ICpp' = ICpD Kpi . 

First note that because the final state is the same in 
both BBi 2 cases fFig. [TOk.b). the final values on the ij 
bond, g" and g" respectively, have to be equal (the initial 
value is g = gij). Applying the rule Eq. (j56p on trian- 
gular faces created by s,s',g edges shows that indeed 
g" = g" = s ■ g ■ s'~^, with conventions as in Figs. I9I10I 
Next, note that fixing the initial and final states amounts 
to choosing a coloring on the surface of the tent-shaped 
body formed on top of the p,p' plaquettes. The only un- 
constrained internal edges are s,s'. By construction of 
the model and the properties of the ground state man- 
ifold, the edge orientations and the constraints on ele- 
ments are consistent with a triangulation and coloring 
of the tent-shaped manifold as required in the definition 
Eq. of Zm{t). The surface coloring r is fixed by the 
choice of initial and final state, while the sum over s, s' 
in the expressions for BBi , BB2 is the sum over internal 
colorings in the definition of Zm{t), Eq. ([43]). The BBi 
and BB2 are therefore equal to the DW invariant Zm{t) 
of the tent-like object in Fig. [TUl and they differ from 
each other only in the choice of triangulation, i.e. the 
position of one internal edge (notice that the value of el- 
ement on this edge is also different in the two cases, but 
in both consistent with general coloring demands from 
subsection nil A\ . According to the properties of the DW 
invariant expressed by Theorem-2 (Eq. (|43)) ). this differ- 
ence in triangulation does not change its value, meaning 
that BBi = BB2. 

IV. ELEMENTARY EXCITATIONS 

In this section we introduce the low energy excitations 
in our models, and study their general properties. We 
define ribbon operators which describe excitations at the 
ends of open strings in Sec. IIVBI having first introduced 
the motivation for the definition in Sec. IIV Al In Sec- 
tion IIV CI we will use the algebra of ribbon operators 
(extended by some local operators) to study the general 
structure of these excitations. The 3-cocycles present in 
our models introduce a "twist" into this extended ribbon 
algebra and therefore play a key role. 

Further, we will explicitly show in Sec. |V] on exam- 
ples that excitations in our models can have distinctive 
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FIG. 10: Overlapping plaquete operators commute. As in Fig. (5] the action and matrix elements of BBi — (f| BpBpi |i) (right) 
and BB2 ~ {f| BpiBp |i) (left) are shown. Note that fixing |f) and therefore g'/j in both cases is consistent, giving g'/j = s-gij-s' ^ . 
Concerning the overall phase factor due to the 3-cocycle factors in Eq. (|5ip . the two images differ only in the choice of the single 
internal edge (dark blue on left, pink on right) of the triangulation of the tent-shaped object. The summation over elements s, s' 
inherent in the plaquette operators, Eq. (|52p . amounts to the sum over internal colorings, while |i) , |f) fix the surface coloring 
of the 'tent'. The phase factor in both quantities -6-61,2 becomes just the Dijkgraaf-Witten topological invariant (Eq. (143 [) ) of 
the 'tent' calculated with different triangulation choices, therefore having the same value. 



symmetry protected properties. Appendix IB] presents in 
detail the general braiding and fusion of quasiparticles 
based on the twisted extended algebra. 

Up to now, the SG and the GG groups in our mod- 
els were either Abelian or non-Abelian. From now on, 
for simplicity we assume both the SG and the GG to he 
Abelian. 



A. Towards ribbons: Loop operators 

We study closed-string (loop) operators in this subsec- 
tion, which will motivate the subsequent expression for 
open string (ribbon) operators. The loop operators we 
will describe commute with the Hamiltonian in Eq. ([57)1 . 
The open-string operators will inherit this property lo- 
cally along their string, except at the string ends, where 
the excitations are located. 

To define a loop operator, let us consider a contiguous 
area A of the lattice. This area is bounded by a sequence 
of connected edges on the triangular lattice forming the 
lattice loop P. Next, if a lattice site i is inside the area 
A, or is lying on its boundary P, we define the plaquette 
p centered on i to be "inside A" , i.e., p £ A. Now, the 
loop operator is just a product of plaquette operators Bp 
inside the area: 

Lp^U (61) 

where the ordering of the product is defined below, al- 
though it is physically irrelevant since the plaquette op- 
erators Bp operators commute for p ^ p' . Obviously the 
Hilbert subspace on which the loop operator is defined 
has the zero-flux rule obeyed in all triangles belonging to 
all plaquettes p £ A. This space is given by ICp^Ci. . .n/Cp^ , 
with pi, . . . ,pa the plaquettes in A. For the purpose of 
this subsection, we can for simplicity consider only states 



which satisfy the zero-flux rule in all triangles of the lat- 
tice. 

The loop operator will, for s = hsS = kg ■ TLsGi G 
GG, have an action only on the boundary P of the area, 
and therefore we label Lp by P only. 

To prove this basic property of the loop operator, start 
by considering the bulk of A, meaning the sites, edges 
and triangles within A including its boundary P. We 
now need to flx the choice of ordering the operators Bp 
in the product Eq. (|5T|) according to their plaquettes p. 
A natural choice is according to the order of lattice sites 
i = p on the lattice, putting highest p rightmost in the 
product. (As before, i = p means the site i on which 
the plaquette p is centered.) This choice turns out to be 
the simplest and most convenient for calculations. The 
action of Lp in the bulk of A is then given by the ex- 
pressions presented in Fig. [TTk . The SG elements on the 
sites are unchanged due to s = 1, while the GG elements 
on the edges get conjugated by s (contrast to Fig. [S]). 
Since we focus on Abelian groups, the loop operator acts 
trivially on edges lying inside A, including its boundary 
P. 

The non-trivial action of the loop operator is therefore 
limited to the edges which lie outside A but still share 
a site with the loop P, and this action is due to the B^ 
operators having p centered on a site on the loop P, see 
Fig.inh. 

Returning to the contribution from the bulk of A, it 
reduces to the phase factors of 3-simplices (tetrahedrons) 
lying on top of the bulk of A, according to the action of 
operators Bp. To evaluate these, we need to consider 
further the chosen ordering of p's in the product. Con- 
sider an oriented edge j ^ i inside A (including P), so 
that its arrow points towards i. By our choice Pp=j is 

applied before Bp^j. The action of Pp=j assigns the ori- 
entation j ^ i to the lifted edge as in Fig. [SI so the lifted 
edge must have its arrow pointing "upwards" (towards 
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FIG. 11: Creating a loop operator Lp = HpgA multiplying plaquette operators centered on sites within area A bounded 

by a lattice loop P. This leads to two types of phase factors in the bulk of A, coming from: (a) Left: down triangles; and 
(a) Right: up triangles. The Bp's in the product are chosen to be ordered according to order of vertices i on which they are 
centered (p = i), with highest p being rightmost (i.e. applied first). For an Abelian group and the choice s = hs ■ IsG the 
elements within A and P do not change, i.e. u'l — U4, = 332, etc., and phase factors from all tetrahedrons cancel up to a 
total phase due to a 2-cocycle evaluated on edges along P, Eq. (|7Up . (b) The loop operator only acts on triangle edges (thicker 
black) lying outside A but sharing one site with P (red edges). 



the lifted i vertex). This property apphes to all edges in 
the plane, so every lifted edge, connecting a vertex to a 
lifted vertex, will have its arrow pointing upwards, i.e. 
towards the lifted plane. 

Having determined this fact, it is simple to determine 
the tetrahedrons formed by action of Lp in bulk of A, 
Fig. [TTk . It becomes obvious that the bulk of the area is 
spanned by only two different types of triangles. For the 
up (A) and down (v) triangles on the lattice, as labeled 
in Fig. [TTK . the resulting phases are 



(62) 



where using the 3-cocycle uj we have introduced the func- 
tion: 



Cs(5i,52) 



^ ^(5, gi, 52)^^(51,52,5) 
~ ^^(51,5,52) 



(63) 



assigned to the 2-simplices (i.e. triangles) on our ordered 
and colored lattice, according to the general considera- 
tions from Section IIII Al More precisely, just as in that 
section, a 2-simplex a is defined by a triangle with or- 
dered vertices and group elements assigned to its edges. 
The 2-simplex is assigned the complex number 



e^(a,^) = iy(a,^)^(-). 



(64) 



where W{<T,ip) = Cs{gkj,gji) for a triangle with ordered 
vertices i j — t- k. The sign e = ±1 is given by the 
chirality, see Fig. [7] and Fig. [T^ 

The bulk of A is formed by the 2-simplices aj , so the 
total phase contributed by the bulk of A is: 



(65) 



with 51,52 G G. The function 0^(51,52) is most generally 
parametrized by an arbitrary group element s £ G, even 
though here s € GG, and it can be directly shown that 
this function satisfies the 2-cocycle condition introduced 
in Eq. (gl). 

On the other hand, the 2-cocycle (cs) value appearing 
in the phases Eq. (p^ exactly corresponds to the phase 



which one can calculate by changing the triangulation of 
this area as in Fig. [T^ (Recall the this re-triangulation 
will not change the total phase due to the allowed 
"moves" from Fig. |51) Namely, we make all the trian- 
gles share the vertex i = N G P. Labeling the vertices i 
on the lattice loop in CW order P = {i \ i = 0,1, . . . , N}, 
see Fig. [T^t, the phase $^ = 9p becomes 



J 



6p — Cs{gN,N-l, gN-l,N-2) Cs{gN,N-l ■ 5Af-l,Af-2,5Af~2,Af-3) ' ' ' Cs{gN,N-l ' ' ' 521, 5lo), 



(66) 



r 



where we took into account the zero-flux rule as well as 
9ij — dji^- E'l- (ESI) shows explicitly that the contribution 
of bulk of area A depends only on the edges along its 
boundary loop P. 



This expression will be useful for us later on, but let us 
now consider in more detail the case when the 2-cocycle 
Co is trivial, which means that it can be rewritten in the 
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FIG. 12: (a) The loop operator Lp, with P the loop of lattice edges bounding area A, has total contribution from triangles 
within A which equals the product of 2-cocycle Ca (Eq. (|67p l evaluated on the triangles as oriented 2-simplices (compare to 
Figs. [Tiandllip. (b) An example loop P (red edges), (c) The phase can be calculated by re-triangulating the area A inside the 
loop P. The resulting expression for the phase Op, Eq. (|66p . is given in terms of elements along edges on P with clockwise 
ordered vertices i = 0, . . . , N . The expression is used for the definition of an open string, Section [IV Bl 



form 



5(51,52) 



gs(gi)£s(ff2) 
£3(51-52) 



(67) 



Because in expression Eq. (|67p the three elements 
5i)52,5i ■ 52 belong to the three sides of the triangle, 
it is easy to check that contributions to 9* from an edge 
shared by an up and down triangle cancel. The total 
phase then becomes a simple expression obtained by go- 
ing around the loop: 

Qp = £s{gN,N-i) £s{gN-i,N-2) ■ ■ ■ £5(520) £s(5lo) = 

= []£,(<7,+i,,). (68) 
cw 

Let us now return to the general case of arbitrary 2- 
cocycle Cg. Having dealt with the bulk contribution of A, 
the action of the loop operator is presented in Fig. [TTb . 
where a segment of the loop is marked by red edges, as in 
Fig-US Note again that along the loop i, i' e P, elements 
Mi, gui are unchanged between the initial and final states, 
just as they are not changed inside the loop. 



1. Summary for loop operators 

In summary, the loop operator Lp only acts on two 
edges per triangle that has a vertex or edge on the loop P. 
This action follows from the action of Bp (see Eq. (|48)) ') 
with p centered on the loop P. (Such edges are marked 
thick black in Fig. [TTb.') Additionally, the loop opera- 
tor has a phase factor contribution from the 3-cocycle of 
tetrahedrons at on top of each such triangle (see shaded 
triangles in Fig. [TTb ) , which we do not explicate as they 
cannot be simplified further; this phase factor we label 
succinctly by 



0-* 



(69) 



Finally, there is the overall phase Op, giving the total 
amplitude for the loop operator: 



n n 93+i..3,9i,i~i 
n£.(5.+i,on^('^*)'^''^ 



(70) 



cw 



teP 



teP 



where the last line holds only for trivial 2-cocycles Cg 
and the CW loop P contains vertices P = {i \ i = 
0,1,. ..,iV}. 

One can show, using a direct method as in Fig. |TU1 that 
the plaquette operator Bp overlapping with the loop P, 
e.g. p Cz P, commutes with the loop operator. 

The final expression Eq. (|70p motivates a definition of 
an open string operator, to which we now turn. 



B. Ribbon operator 

In this subsection we introduce the open-string (rib- 
bon) operators, which describe the low energy excitations 
in our models. 

Let us start by defining a geometric object: the open 
ribbon F. The ribbon has two ends, end-A and end-B, 
which we need to define first. Choosing two neighboring 
vertices on the triangular lattice, vertex-i^ and vertex-i^, 
there is a unique 2-simplex (triangle) formed by vertex-z^i 
and vertex-z^ and another vertex that is not contained 
in F. We denote this 2-simplcx by tA- We then use the 
label A = (i^, tA), the collection of the vertex-i^ and the 
triangle tA, to formally define the end-A of the ribbon F. 
Similarly, we will use B = {iB,tB) to define the end-B 
of F. The 2-simplices tA,tB are not inside F. We also 
define the edges of ribbon F at the two ends: for end-A, 
there is a single 1-simplex that is shared between F and 
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Ia, which we denote as edge-^. Similarly we can define lattice edge ij inside F, such that it connects the inner 



edge-B. 

Having defined its ends, the open ribbon F is finally 
specified by its two "ribbon edges". Namely, F has an 
"inner-edge" , which is a sequence of connected edges on 
the triangular lattice, going from the vertex-M to the 
vertex-is. Further, F also has an "outer-edge", running 
from the vertex-i^ to the vertex-i^, which is displaced 
from the inner-edge of F by one lattice spacing. (See 
Fig im for an illustration of these definitions.) As a ge- 
ometric object, the ribbon F contains all the vertices on 
the inner and outer edges, together with all 1-simplices 
and 2-simplices connecting these vertices (for these sim- 
plices, we will say that they are inside F, or we will write 

We now proceed to define the operator ^('''^^(F) for a 
given open ribbon F. 

Fet us define a Hilbert sub-space /C(F) C TL, formed by 
those states which satisfy the zero-fiux rule everywhere 
inside F. Namely, V IV^) G /C(F) and Vt e F, Qt \tp) = 
(Note that however, may violate the zero-fiux rule at 
tA^ts, for instance.) For the purpose of this paper, we 
only define the operator (F) in the Hilbert sub-space 

/C(F). 

Based on the understanding of the loop operator 
from previous subsection, we define the ribbon opera- 
tor F^'^'S^r) such that it modifies the gauge degrees of 
freedom living on the lattice edges connecting the inner 
and outer ribbon edges of F, while leaving all degrees of 
freedom living elsewhere unchanged. In particular, for a 



and outer ribbon edges, the group element h^j living on 



it is changed into h'^j — n ■ ihj [n-ij — ihj 
lattice edge is oriented to point towards the inner [outer] 
ribbon edge. The operator F^''^'^'>iT) therefore has non- 
zero matrix element only between states {/i^j}) and 
\{ui},{hij}). 

Finally, we need to define the non-zero matrix element 
of F'^'''^)(F). This matrix element has two factors: one 
chosen in accordance with the closed-loop operator, and 
the other dependent on the degrees of freedom at the two 
ends of the open ribbon F. We define: 

({"J, {M,}|F('''^)(F)|{uJ, {K,}) =fA-fB- fAB ■ wlig), 

(71) 

where: 1) wj^{g) is motivated by the closed- loop opera- 
tor and will be presented shortly; and 2) fA,fB,fAB are 
rather complicated phase factors depending only on the 
degrees of freedom living on ends of F, and we present 
them in Appendix |B] To motivate the expression for 
wj^{g), let us start from the expression for loop opera- 
tor Eq. ([70)1 . Fet the F ribbon's inner edge go along the 
sequence of lattice sites {i \ i = 0,. . . ,N}, where now 
the sites i = and i — N are not nearest neighbors, 
but rather the vertex-is and vertex-i^, respectively. For 
convenience (see Fig. [T51 for the pictorial definition), we 
define the group elements a„ G G, n — l,...,N hy 
Oi = gi.i~i- We then define the phase 



h ■ h^j [h[. 



if the 



N 



(72) 



vl{g) = ^^e^^\^^g,,,.^,gj ee 

= (^Yl ^i'^tY'^'^'^^ Ch{t,aN)ch{aN,aN^i) Ch{aN ■ aAr_i,ajv-2) • • •c/i(ajv • ■■ 02,01) S{ai ■■■aN,g) = 

= WicTtY'^'^*^^ Ch{g ■ aj"\ai) Ch{g ■ aj"^ • a^\a2) • ■ ■ Ch{g ■ aj"^ • ■ ■ a~^\,aN-i) Ch{g ■ aj"^ • • • a~\ ca?) (5(ai ■■■on, 



with 5{g,g') the Kronecker delta function, which we used 
in the second line to obtain a simple pictorial definition. 
Fig. [131 

A way to understand the meaning of phase in Eq. (|72[) 
is to note that "cutting open" the loop to get an open 
string introduced the parameter g € G which is related to 
the charge carried by the excitations at the string ends. 
We will consider it in more detail below. (Compared to 
the loop, there is one extra factor c;i(l,aAr), which is 
inconsequential for canonical cocycles, see Eq. (1870^ .) 

In accordance with the definition of loop operator, the 
element h S GG C G, see after Eq. ([CT|) . Physically, 



the loop operator Lp can be seen as a closed, P loop- 
shaped, domain wall inside which we acted by element h 
(through action of 13^ ). For a gauge element h G GG 
there is actually no transformation inside the domain 
wall. It is further possible to create an 'open domain 
wall' (i.e. open string F) in the gauge transformation h, 
and it defines gauge excitations at the ends of F. (On the 
other hand, one could define Lp for a global symmetry 
operation h G SG, which would create a closed domain 
wall inside which the transformation h G SG is applied. 
However, it is physically unsound to try to define an open 
domain wall of such a transformation h G SG.) 
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FIG. 13: The phase contribution wl{g) = *f 6f 5 (^11™=! 9 J matrix element of ribbon operator F'^''-^\T), with s £ GG, 
on a length N = 5 example. (Left) The operator acts non-trivially only on elements adjacent to path F (thick black edges) 
Note that elements m = gi,i-i = hij-iUiU~}-^ axe defined to be directed along the path, so for n = 4, 5 they are opposite to the 
standard definition, which is along the edge orientation. Elements at ends A, B of the ribbon F, i.e. vertex-iA and vertex-i_B, 
are us € SG and ito € SG, respectively. The phase factor 'I'f is the product of 3-simplex phases W{aiY'''^'^ shown on top of 
shaded triangles. (Right) The phase Gf is the product of 2-simplex phases W^(crt)^'''^'' shown, where 1 is the identity element, 
and the parameter g = a\ ■ ■ ■ aN £ G. 




FIG. 14: Inner segment of ribbon operator _F'°'^^(F), with s G GG, commutes with plaquette operator Bp having p centered 
on inner edge of F, and therefore also with the Hamiltonian. (a) Only elements on two edges (thick black) are changed, but in 
the same way for F13 and BF due to Abelian G. (b) Quantum phase for FB (top row) and BF (bottom row) differ only due 
to shown 3-simplices which are on top of three blue shaded triangles in (a), (c) The phase contribution due to string phase O 
is shown, (d) The 3-simplices contribution to phase factor ratio WB/BW (from panel b) equals precisely the phase W {atY^"*'' 
due to two blue 2-simplices (see Eq. ([78])). The WB (top) and BW (bottom) case are therefore graphically seen to be the 
same, as interior points of the polygon can be removed in a 1-3 move (Fig. (H)). 



It is important to realize that the delta function in 
the ribbon operator, 6{ai ■ ■ ■ aj\j , g), with = gi,i-i — 
hi^i-i-Ui-u~}i, has different effect on the global symmetry 
and gauge parts of G = GG x SG. Namely, given the 
factorization g = hg ■ g, the delta function separates into 

S{ai ■ ■ ■ aM,g) ^ 5{hio ■ ■ • hM.N-i, hg) S{uMU^^,g). 

(73) 

The delta function therefore constrains the product of 
gauge degrees of freedom along the inner edge of ribbon 
to the value hg, while leaving only one constraint on the 
undualized elements of SG at the two ends of the ribbon: 
UN = 9 ■ uq. (Recall that actually un is the element on 
site vertex-iyi, while uq is on vertex-is.) 

A special case that offers insight occurs when the 2- 
cocycle Ch is trivial. In that case, using the property 



Eq. (p7|) of trivial cocycles, we get 

<{9) = \ 5{a,---aN,g)el\g)\{eH{a,), 
Vter / i=i 

(74) 

where Eh is a family of U{1) valued functions on G (i.e. 
a 1-cochain) parametrized by the element h. 

We now show that away from its endpoints the ribbon 
operator F''^'^'>{T) commutes with , and therefore also 
with the Hamiltonian in Eq. (|57p . Obviously, the non- 
trivial situation occurs when the plaquette p is centered 
on a site z € F which is positioned on the inner edge of 
r, see Fig. [Tlk. Note that by definition of the operators, 
elements s € GG and s' G G. The product of the two 
operators is defined in the Hilbert subspace A^(r) fl ICp. 
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Physically, we do not allow states with a fliix-carrying 
excitation positioned on the ribbon F or inside the pla- 
quette p, as we consider this commutator. We also choose 
p that does not overlap with tA or Ib, i.e. we consider 
the commutation away from ribbon ends. 

It is clear that the final states of the system are the 
same no matter the order of applying the two operators, 
since G is Abelian. The two resulting quantum phases 
(f| F^''3^r)B^' |i) and (f| i?^'F("'S)(r) |i) we label by wB 
and Bw, respectively, and we will show that their ratio 
is 1. 

The wB and B w differ in their phase factor 8f (which 
was defined in Fig. [T3| , and these are shown on the right 
of Fig.fTlTb). RecaUing that Ui = gi.i-i = ft-i.i-i ■Ui-uY_}i 

and i — 1, . . . ,N, the difference is due to operator B^ 
sending 



a, s ■ Qi 



0.1+1 fli+i • s' 



1-1 



(75) 



see the basic definition, Eq. p8)) . and we remind that p 
is centered on the site i. Due to the product structure 
in 0p, there are only two factors affected by this change, 
and these are the shaded 2-simplices in Fig. [TUb) (right). 
The ratio of quantum amplitudes is 

QwB _ Cs{A ■ s'~i • a^^s' • at)cs{A ■ a^"^ ■ ajr^\,ai+i • s'' 

©Bu. 



Cs{A ■ a^ ^,ai)cs{A • a ■ ^ • a^_^\,ai+i) 



(76) 



Csjs' ^ ■ at+i,s') 
Cs(,s',aO 



where A is the product of ai'<i elements, and we used 
the 2-cocycle condition (see Eq. ©). 

Further, it is obvious that the phase difference coming 
from 3-simplices is only due to ones stacked on top of 
three triangles that are shared by the ribbon and the 
plaquette. Fig. [TiT a'). and these simplices are presented 
in detail in Fig. iMT b). Grouping the phases according to 
the triangles from left to right, and using the definition 
of the 2-cocycle Eq. (|63|. we get 



f— ^ Cs{ai+i,h) {s' ■ ai+i,s' ■ h)LO {s',h,s) 



Bw 



^2 



^1 



wB 



(77) 

: w(s', s, h) u)^^{s, s', h) U!{s, s', • h) uj{s' , s, ai ■ h) 



Cs(s', flj) a;(s', s, ■ h)uj ^{s, s , ■ h), 



where h is the element on the leftmost edge of left triangle 
in Fig.fTHa). The total phase due to 3-simplices becomes 



Bw 



Cs{s' 1 • ai+i, s') 



(78) 



This phase ratio exactly cancels the contribution from 
the string of 2-cocycles in Eq. ((75)) . completing the proof 



that the inner piece of ribbon operator commutes with 
plaquette operators. 

We point out the pictorial interpretation of the re- 
sult, that the 3-simplex contributions to the phase ratio, 
Eq. ([75)) . are exactly equal to the phase of the two blue- 
shaded triangles in Fig. [MlcV The total phase ratio is 
then equal to the ratio of the upper and lower polygons 
in Fig. [T3^c), which obviously equals 1 due to the rule 
that allows removal of internal points in polygons (rule 
from Fig. HI). 



C. Local symmetric operators and the twisted 
extended ribbon algebra 

In this subsection we will introduce local symmetric 
operators which have a non-trivial algebra with the rib- 
bon operators. This will allow us to understand the 
general structure of excited states. The focus will be 
on states with a single pair of excitations, or two pairs 
when discussing braiding. The obtained results will be 
directly relevant for studying examples in Sec. [V] The 
fully general case of many quasiparticle pairs, including 
their braiding and fusion properties, will be studied in 
detail in Appendix |B] by using the extended algebra. 
^ Before giving the formal definitions, let us remark that 
given the positions of excitations, the "extended algebra" 
contains the set of ribbon operators F (with their strings 
connecting pairs of excitations), as well as a set of local 
operators D acting at the positions of excitations. How- 
ever, the presence of the 3-cocycle u adds a "twist" in 
the algebra and crucially determines the resulting prop- 
erties of excitations. It has been shown^^ that certain 
broken gauge theories with Chern-Simons terms lead to 
discrete (group H) gauge theories having such twisted 
algebra describing their quasiparticles^S.. In that situa- 
tion, the cocycle is generated by the Chern-Simons term, 
and the resulting discrete gauge theory can be classified 
using H^{H ,U{1)),^* making connection with discrete 
DW TQFTs. Importantly, a non-trivial cocycle twist- 
ing of the algebra can render the resulting theory non- 
Abelian even though its gauge group H is Abelian^. Our 
models inherit this interesting property, and we will dis- 
cuss this briefly concerning properties of multiple pairs 
of excitations in Section IV Gl In Appendix |B] we will 
also show explicitly how the considered operators of our 
models form a Hopf algebra (more precisely, a quasi- 
triangular quasi- Hopf algebra due to the cocycle twist, 
see Refs lssllGlh . and how they describe the braiding and 
fusion of excitations^. 

It is important to emphasize that in the present context 
the local operators D are also crucial for determining the 
interplay of symmetry and topological order: After we 
explicitly construct the D operators, we will show that 
they are symmetric, i.e. commute with transformations 
from SG, and further we expect them to span the alge- 
bra of all local symmetric operators. Thereby the set D 
will provide us all symmetry allowed local perturbations 
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FIG. 15: (a) Ribbon operator with both ends (red dots), shown schematically in the bottom. It only acts on the thick black 
edges, consistent with zero-flux rule on all 2-simplices above and within the ribbon. The transition amplitude of the ribbon is 
given in Fig. [O] without the phase factor /a/s/as local at the string ends; / is involved (Appendix|B|, but includes the edges 
of two triangles tA,tB at the ends (shaded light orange). Dashed violet lines indicate edge elements that are not acted on by 
the operator, (b) The local operator at position B = (i,t), with i the lattice site and triangle t (dark blue) having vertices k 
(shaded dark blue), has the transition amplitude (f| g)(B) |i) — 5{gij ■ gjk ■ Qki, h) Wh Wa. It acts as the plaquette operator 
I3p (Fig. |9]) with plaquette p centered on i (changing only thick black edge elements), except that it also projects flux in t to 
value h, and also has an additional phase factor Wh which is the phase of the 2-simplex (shaded light blue). (For operator 
D{A), the triangle iki' is used instead.) The phase factor We depends on six drawn 3-simplices (tetrahedrons) having altered 
elements o (red edges) such that it is well-defined even if the zero-flux rule is violated in the plane (as occurs at ribbon string 
ends), see after Eq. (17^ . 



and therefore the possibility of calculating the symme- 
try protected degeneracy and other properties of excited 
states in Section IVl 

Let us briefly recall some relevant details about the rib- 
bon operator F^^'S) (F) from Sec. HVBl see Fig. [HI By 
definition, the element h G GG while g € G. The oper- 
ator definition demands that the zero-flux rule is satis- 
fled for all triangles in ribbon F, i.e. it acts within the 
Hilbert subspace /C(F). The two ends of ribbon F we 
label by A, B. For this subsection it is important to re- 
call the structure of the ribbon ends (Fig. [15]), which are 
completely determined by a site and a lattice triangle 
A — {iA,tA), B — (j_B,t_B) {tA,tB are not considered to 
be within F). 

We next deflne the local operator Df^^^g) [B) positioned 
a.t B = {iB,tB), Fig. [ini (Operator D{A) will be very 
similar, below.) Let the triangle ts have ordered ver- 
tices fc — > j — > i as in Fig. [TS] The D(^h,g){B) operator 
acts on the elements in the plaquette centered on is in 
the same way as the plaquette operator B^^^^ ; however, 
D(h.g){B) additionally projects the 2-siniplex at Ib to 
having flux h G GG, and also has an additional phase 
factor. Actually, in contrast to B^=ig , we define the op- 
erator D(^iig){B) in the entire Hilbert space H, as will 
become clear soon. Let us first state the quantum am- 
plitude of the local operator: 

(f| i^(„,3)(B) |i) = S{g,, ■ g,k ■ 9H, h) Wh W^ii), (79) 

where Wh = Wh{fTij,(pY'''^^^^ is the phase of 2-simplex cr.y 
(light blue shaded in Fig.llSp formed by edge ij and verti- 



cal edge g, so it equals Ch{g,gij) (see Fig. [15]). Of course, 
by the definition of group elements on edges, Eq. (j45l) . 
Hdij ■ 9jk ■ gki,h) = S{hij ■ hjk ■ hki,h) depends only on 
the elements of gauge group GG. 

The phase We{i) is due to 3-simplices (tetrahedrons) 
on top of the plaquette; although analogous to the I3p 
operator (Fig. [9]), in this case operator D has to be well- 
defined even if the zero-flux rule is violated in the plane, 
as can occur at the tA,tB of the ribbon F. On the other 
hand, the 3-cocycle can assign a phase to a tetrahedron 
only if the zero-flux rule is satisfied on all its faces. Be- 
cause of this, just for the purpose of calculating the phase 
due to the six tetrahedrons Il/=i ^{'^i^ 

rede- 
fine the values oi , . . . , og of elements on six outer edges 
of the plaquette (red in Fig. [TSb) such that the zero-flux 
rule is satisfied in all six triangles of the plaquette in the 
plane. The six internal edge elements of the plaquette in 
the plane are considered unchanged from their value in 
i), and they suffice to fix the redefined values o[, . . . ,o'q 
on the outer edges of the plaquette according to the zero- 
flux rule in all six triangles of the plaquette. Formally, 
this phase contribution is: 

Weil) = (o'l, . . . ,Og; /i, ...JelBp \o[,. . .,o'Q;ii, . . . ,^5) , 

(80) 

with ii, . . . , ig the initial values of elements on six internal 
edges of plaquette, and /i, . . . , /e their flnal values. (Note 
that the redeflned values on outer edges o[, . . . ,o'q have 
to stay the same in the initial and final state.) 

It is now clear that operator D(^h,g){B) is well defined 
in H. Actually, it is also well-defined within the subspace 
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/C(r) because it does not influence the flux in triangles 
belonging to ribbon F (recall that by definition Ib is not 
inside F). 

The definition and properties of the D{A) operator, 
relevant due to its action at the A-enA of a string, is 
almost identical to the D{B) case just described. The 
only difference is in the phase Wh^ which in this case is 
not the phase of the 2-simplex iji' (light blue triangle in 
Fig. [TSb). but the 2-simplex iki' instead. 

Obviously, a non-trivial algebra between D{A), D{B) 
and F is due to their overlap at the triangles tAjtg. For 
concreteness, we presented ribbon ends of the form in 
Fig. [151 i-e. tA,tB being bottom up-triangles, omitting 
versions rotated by multiples of 60°. 

The local operators D(iig){C), with C — A,B, are 



symmetric. As noted before (see discussion of plaquette 
operator. Sec. IIII B|) . a global symmetry transformation 
s G SG does not alter the elements on the edges, leaving 
the phases in Eq. (1751) intact. Further, even though the 
local operator acts on the site ic by the element g g SG, 
where g — hg ■ g, hg G GG, this action automatically 
commutes with the action of s since we have restricted 
G to be Abelian in the present analysis of elementary 
excitations in our models. 

We emphasize again that for local operators -D(ft,^g)(C), 
with C = A,B, by definition h E GG while g £ G. 

Having the definitions of: ribbon F having ends C ~ 
A,B; the ribbon operator -F'(F); and the local operators 
D{A), D{B) at hand, one can tediously but straightfor- 
wardly derive the following algebra: 



^(/i2,S2)(C') £'(;ii,gi)(C') - 5h^^h2 Chi (52,51) P^{h^,g^g2){C) 

F(^^^s.a.)iY) i?(.,,.i,,,)(A) =c,,(/.2,/ii)c,i(52,.9i)i?(.,,3,)(A) F^^-'^^\T) 
Dik2hu92)iB) F^'^-^^^3-\T)^Cg,{h^M)chA9i.92)F'^'^-^^-\T) D^„^^g.^~,{B), 



(81a) 
(81b) 
(81c) 
(81d) 



where Eq. (|81ap holds in the entire Hilbert space H (and 
also within /C(F)), while the three other relations hold 
within the subspace /C(F), where the ribbon operators 
are well-defined. 

To understand the implications of the operator algebra 
in Eq. (|8T|) . we need to consider the precise form of the 
excited state. Starting from a two-particle excited state, 
having excitations at the two ends A, B of the string F, 
one would expect it to be given by the simple action of 
the ribbon operator on the ground state: 



^('1,9)^ = F(''.9)(r) \gs) 



(82) 



These states can be shown to be orthogonal. However, 
the space C{A, B) spanned by these states needs to be 
specified further. Namely, as Eq. (|73p shows, the ribbon 
operator puts only one constraint on the values of two 



elements at the lattice sites lAjis, i-e. Ui 



= 5 e 



SG, with the factorization g = hg ■ g. We already know 
that the local operators transform these elements, e.g. 
under action of Di^h.gi){A), the element 51 • Ui^, 

where 51 = ft-i -51, /ii G GG, 51 G SG. We therefore need 
to specify the value of one element (either of Ui^ or ) 
in the excited state. By using the projectors 



Pu{A)\{u,},{g,j})^ 5ui.,u \{ui}, {gij}) 



(83) 



we can consider the subspace {A, B) of the Hilbert 
space % spanned by projected states: 



ih..g) 



(84) 



with a fixed element ua G SG. The value of Ui^ is 
then automatically fixed by the action of ribbon oper- 
ator = g~^ 



Completely 



UA. 

analogously 



Cub B) spanned by 



define the subspace 
A.,(B)|7/-(''^5))^byus- 

ing the projector P^^ {B) at the vertex-i^. Note that the 
projectors P{A), P{B) commute with the ribbon opera- 
tor F{T). It is also easy to check that Ds^h^g) {G)Puc (G) = 
Pg-uc{C)B>{h,g){G) , where G is either A or B, and the 
usual group element factorization is g ^ hg ■ g, hg G GG, 
gESG. 

The subspaces Cue B) were introduced using the 
action of ribbon operators on the ground state. The 
end-^ and end-B of ribbon F in principle carry excita- 
tions, and Cue (^i B) does not depend on the particular 
shape of F (as long as F does not change topological 
class). It is however not trivial to prove that the sub- 
space C{A,B) = (BucGScCuci^i B), with G either A or 
B, actually exhausts all possible excited states with two 
quasiparticles positioned at A and B. In Appendix |E] 
we prove that the space C{A, B) indeed contains all such 
excited states. (We do not however have a proof that all 
multi-particle states having more than two particles, are 
also given by the action of appropriate ribbon operators 
on the ground state.) 

Now that we established the appropriate Hilbert space, 
it is easy to show that the local operators form a unitary 
projective representation of the group G within this two- 
particle Hilbert space. We give the precise definition and 
explicit proof of this fact in Appendix [El 
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Let us next show that the local operators, as well as 
the ribbon operators, form Hopf algebras. (The succinct 
notation we introduce here will be useful in Appendix IeI) 
More precisely, let us denote the algebra formed by local 
operators D(j^g'f{A), with h G GG, g € G, hy the sym- 
bol D{A). Analogously, the operators D{^}ig-^{B) form 
the algebra Since n{A) and are formally al- 

gebraically the same, we use D to denote this abstract 
algebra, keeping in mind that and D(_B) act at 

physically different positions. The algebra of ribbon op- 
erators fC'-s), h G GG, g &G, having ribbon T, wiU be 
denoted by F. We now consider the subspace /C(r) of 
Hilbert space, so that all these operators are simultane- 
ously well-defined, and Eq. (|8T|) holds. 

To start with, the operator relations in Eqs. (|8ip can 
be rewritten succinctly using double index notation: the 
Latin indices i,j,k,...,r are shorthand for group element 
pairs i = {hi,gi), j = etc., with hi, hj, . . . e GG 

and gi, gj, . . . £ G. We can then write: 



F™(r) F"(r) = A^'i^^r) 

F'"(r) D,{A) = Af nil Dj{A) f\t) 
D,{B) F^{T)^Kl^nTuF\T)D,{B) 

where we defined the tensors 



(85a) 
(85b) 
(85c) 

(85d) 
(85e) 



(86a) 
(86b) 



and summation over repeated double indices is under- 
stood as, e.g., Y.^ = Eh,eGG,g,eG- 

The tensors in Eq. ([86)) contain the twist due to the 
cocycle uj present in our models. Similar realizations of 
twisted algebra, but describing broken gauge theories, are 
analyzed in e.g. Ref.62. We also note that only when the 
cocycle is trivial, uj{gi,g2,gs) = 1 ^ Cg{g',g") = 1, these 
tensors reduce to the form occurring in the generalization 
of toric code to arbitrary finite groups, as given in Ref . Issl. 

To derive the properties of this algebra, we will have 
to use several properties of a 2-cocycle c derived from an 
arbitrary 3-cocycle uj (Eq. (pH])). listed here for conve- 
nience: 

Ca{g, h) Caigh, s) = Ca{h, s) Caig, hs) (87a) 

Caig, h) Cb{g, h) Cgia, b) Ch{a, b) = Cgh{a, b) Cab{g, h) 

(87b) 

Caig, h) ~ 1, if any of a,g, h equals 1 (87c) 
ca{g,g-^)^caig-\g). (87d) 

(We note that all these properties actually hold generally, 
i.e. for arbitrary elements a,b,g,h,s G G.) The identity 
(|87ap is the general condition for a 2-cocycle, (j87bl) can 
be derived directly using the definition Eq. (|63p . identity 



(|87cp follows from our choice to use only canonical cocy- 
cles a; (U is the identity group element), and finally the 
useful relation (j87dp follows from (|87ap and ()87cp . 

Using these identities we can prove that the ribbon 
and local operators form Hopf algebras. A basic axiom 
of Hopf algebra is the associativity of multiplication with 
identity, which holds for F and D algebras, as shown by 
the following two relations, respectively: 



A^"A^' 



A Ij A mn 
^^k ^^'J ' 



e A™^ 



^t, (88b) 

where for double indices the Kronecker delta function 

5) = Sh.MjS, 



and the functions e,- = S 



= 5, 



define the unit and counit of the algebras, If = ^kF*' , 
Id = e'^Dk, e(F'=) = e'= and t{Dk) = e^. Both Eq. (Pa|) 
and Eq. (pb)) hold due to Eq. ((57a|) and Eq. (|87cll . 

The comultiplication in the Hopf algebra is physically 
related to fusion, and instead of presenting a formal def- 
inition here, we show in Section IIVDI and Appendix [B] 
that the braiding and fusion properties of excitations are 
consistent, and contain the twist by the 3-cocycle charac- 
teristic of a quasi-Hopf algebra first introduced in Refl6ll. 

Having established the Hopf algebra relations in 
Eqs. (|88ap . (|88bp . we are in a position to prove that 
all two-particle excited states in our models are indeed 
within the Hilbert subspace C{A, B) we defined in this 
subsection. The proof is presented in Appendix [El and is 
based on a derivation given in Ref. |58|. This Hilbert space 
for an excitation pair is especially important here since 
it will be studied explicitly for several example groups in 
Section El 

In summary, the local operators are the set of non- 
trivial operators acting on excitations, and they are also 
symmetric and form a projective representation of the 
group. We will use these properties to study the symme- 
try protected properties of our models in Section |Vl 

We close this subsection by briefly introducing the 
Hilbert space of many-particle excitations. (It is stud- 
ied in more detail in Appendix |Bl) Let us consider a 
system having n quasiparticles at positions Ai, . . . ,An, 
and one quasiparticle at position B, and no other exci- 
tations. Such a system is described by the Hilbert space 
£(Ai, . . . , An,B). To study these states, we consider a 
space C{Ai, . . . , An, B) spanned by the action of ribbon 
operators, as described in the following. Let us connect 
each position Ai through a ribbon Ti having end-A^ to 
the common isolated position B. Therefore all ribbons' 
cnd-i3i coincide, and all are equal to B. (For further dis- 
cussion and lattice realization, see Appendix iBl) Focus- 
ing on the subspace lC{Ti , . . . , r„) in which all ribbon op- 
erators having ribbons Ti are simultaneously well-defined 
(i.e. zero-flux rule obeyed inside all ribbons), we can 
now define its subspaces Cub (^ii ■ ■ ■ : ^) spanned by 
states of the form: 

l^^y ■■''") = Pus {B)F^- (Fi) • • • F'^" (F„) \gs) . (89) 

The states in Eq. ([551) form the subspaces 
with ub fixed. We actually expect that the 
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with hj , hg^ G GG. This constraint to GG elements is ex- 
actly encountered when describing braiding in the next 
subsection. 

Because the Dj, operators are non-local with respect 
to the excitations at Ai while commuting with ]3){Ai) they 
are called "topological operators"—. 

Topological operators and braiding in many-particle 
states are analyzed in detail in Appendix [BJ 



FIG. 16: a) Ribbon operators for two-particle excited state, 
and braiding two excitations at Ai, A2 (see also Fig. llSf) . The 
common endpoint B contains no excitation, but "topological') 
operators of i-th ribbon D^'^^{B), i = 1,2, local at B, can 
be used to represent braiding, b) Applying different ribbons 
(blue) than in original state shown in (a) leads to a braided 
state (see Section IIVD|l . c) The resulting counter-clockwise 
braid of particles 1 and 2. 
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B) 



coincides 



with the 
...,A„,B) 



space 
(We 



have a proof of this only in the two-particle case, see 
Appendix |e1 but as mentioned this is hard to prove in 
general.) 

At the same time, the extended algebra of each given 
ribbon F{Ti) contains the local operators D^*-* = Dk{B) 
which we define to affect only that (i.e., the i-th) ribbon. 
These operators, being local at B, commute with all local 
operators Djn{Aj), j = 1 . . . at the excitations Aj. The 

-D^*' operators therefore act, in view of Eq. (jSSdp . as 



D 



(1). 



(n) 



/ fci,...,fc„\ _ Qfei r. 



(90) 

In this definition we constrain the elements in D{B) 
strictly to GG, e.g. the double index j — {hj, hg. ■ Isg) 



D. Braiding matrix 

In this subsection we calculate the braiding matrix of 
two quasiparticles, restricting ourselves to a system that 
has quasiparticles at most at three positions, Ai,A2,B 
(see previous subsection). An alternative and more ex- 
plicit way to obtain the braiding properties of quasipar- 
ticles is by considering ribbons of strings that cross, but 
this approach is applicable to Abelian quasiparticles only; 
we present it in Appendix [Cl and show there that the 
results of the two approaches agree for Abelian quasi- 
particles. The braiding as well as fusion in states with 
many quasiparticles introduce additional subtleties, and 
this situation is presented in detail in Appendix |B] Here 
we will focus only on braiding of two particles, and show 
that the braiding properties are entirely determined by 
the topological order, i.e. the gauge group GG. The ef- 
fects of interplay of topological order and symmetry are 
revealed in the concrete examples in Section |Vl 

Recall that we define the ribbon operator matrix ele- 
ment (see Eq. (|71|)) such that the operator algebra in Sec- 
tion |IVC] is obtained. (The matrix element is fully pre- 
sented in Appendix [BJ) That definition also leads to the 
following algebra for two ribbon operators having strings 
ri,r2 which share their _B-end: 



u[hi,h2, tiB) 



(91) 



r 



where /is is the value of flux in triangle Ib- The Eq. ([9T|) 
is well-defined only in the subspace /C(ri,r2) of the 
Hilbert space, i.e. when zero-flux is obeyed in every trian- 
gle inside Fi or r2 (note that is is outside both ribbons). 
Fig. [TOk sketches this situation, which will enable us to 
determine the braiding of the two quasiparticles at the 
end-Ai and end-A2- (Appendix IB] considers in detail the 
case when the end-A is shared by the ribbons.) 

The 3-cocycle factor on the right hand side of Eq. (pij) 
has to be considered with care in a many-particle situ- 
ation, as is done in Appendix |B] in this subsection we 
however focus on a system with only the two ribbons 
ri,r2. Because of this, a product of two ribbon oper- 



ators, such as appearing on both sides of Eq. (|?1]) . can 
only act on the ground state, i.e. there are no other ex- 
citations in the system. Because of this, the value of flux 
in triangle is on the right-hand side of Eq. (|9T|) is nec- 
essarily zero, which means that the element /is = 1 (see 
also Appendix [B]) . Due to our choice of canonical cocy- 
cles, the two w factors on the right of Eq. (pij) therefore 
disappear. 

Let us now formally introduce the braiding matrix rele- 
vant for ribbons sharing their end-B by rewriting Eq. (j9ip 
in a compact form: 



F\V2)F'{Vi) = R'i'-^l 



(ri) F"(F2 



(92) 
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where according to Eq. (191 
equals i?*'^ = Sh,_n,.S 



the i?-matrix of our model 



We now calculate the operator 72.cc describing the 
counter-clockwise braiding operation by 180° of the two 
excitations positioned at Ai and A2. This operator is 
defined by setting |V'°'')^^„,-^ = T^cc We consider 

the original state, Fig.llob. and the one where excitations 
are exchanged (Fig. ITBb ) by braiding particle 2 counter- 
clockwise around particle 1: 



\r.l),r^,, = PuAB)ncc\r'') 

^P,AB)F^{T[)F'>{n)\9s), 



(93) 
(94) 
(95) 



where we used the double index notation, i.e. a = 
{ha,ga),b = (hb^gb). Note that we choose to project 
the states by fixing the value of Ui^ equal io ub — the 
ribbon operators then determine the values of , 
(see the definition in Eq. ((84)) ). 

According to the definition of states in Eq. (|89|) . the or- 
der of applying ribbon operators reflects the order of par- 
ticles, while the new strings r'2^,r2 compared to original 
ones Fi, F2 represent the braiding movement. Figs. 116b . b 
illustrate this. We see that the new string of particle 1 
coincides with the old string of particle 2, F'j^ = while 
the new string F2 is topologically equivalent (and the 
operator therefore the same) to Fi only if there is no 
particle at original position of 2 at the time these two 
strings are compared^ see Fig. 116b . Therefore the new 
ribbon operator on Fj is the same as a ribbon on Fi if 
it is applied before the new T'l. This is actually true in 
Eq. ([M)) . and we therefore have 



(96) 



We now only need to commute the ribbon operators to 
make a comparison to Eq. (j93l) . The commutation rela- 
tion in Eq. ([M)) directly gives 



PuAB)ncc\r'') 



R'^nl^K^ |VC"> = (97) 

1^^^). (98) 



We note that the R"^^ matrix in this expression constrains 
the group elements in D^^'^'> to be strictly in GG, e.g. 
r = {hr, h'^-gr), hr,h'^ G GG, is constrained to 5^ = Isc! 
this implies that the site element is not changed from 
its value ub by the action of D^^^ , D^^^ . The explicit form 
of the braiding operation is therefore 



(99) 



which one can of course obtain directly using the com- 
mutation in Eq. (I^TI) . without first defining the R^^ ma- 
trix. The result non-trivially involves the cocycle uj of 
our model. 



Applying 72.cc twice, we obtain for the 360° braiding 
T^lc\K) = 



= ChAgaK ^^hb)chb{gbha^,ha) 



(100) 



(ha,hb) (hbjla) 



where we remind that D^g' = Dg{B) is acting on the i-th 



ribbon operator in the product state Eq. 

The case when the braided quasiparticles are posi- 
tioned at i?-ends of their ribbons which share the end-A, 
is similar and discussed in detail in Appendix |B] Here we 
just quote the result for the 27r braiding of particles at 
end-i?i and end-i?2: 



C^), (101) 



where we use the bar over symbols to signify that the 
state and braiding concern particles at end-i?'s of strings. 

Eqs. (jlOOp and (jlOip explicitly show that the braiding 
of quasiparticles is described by the action of "topological 
operators" (see Eq. (|90|) ). which act only on the gauge 
degrees of freedom, obviously since they are labeled only 
by elements ha, hi, G GG. This fact means that the braid- 
ing properties follow directly from the topological order 
in our models, which is described by the gauge group 
GG. In Appendix [B] we show explicitly that the topo- 
logical operators form an algebra called "quasi-quantum 
double"—!^, which is mathematically a realization of a 
quasi-Hopf algebra. The prefix "quasi" denotes the pres- 
ence of the twist by cocycle uj, which is here restricted to 
elements of GG. Appendix [B] also clarifies how represen- 
tations of the quasi-quantum double label the quasiparti- 
cle species, while the multiplication and comultiplication 
operations in the algebra directly determine the braiding 
and fusion of quasiparticles. 

The quasi-quantum double construction, i.e. the pres- 
ence of a cocycle "twist" in the algebra of braiding and fu- 
sion operators, appeared in the description of excitations 
in gauge theories broken to a discrete subgroup^. This is 
to be expected since both those theories and our models 
describe topological order classified by DW TQFT g — ' — . 

In Appendix |B] we also discuss in detail the situation 
with an arbitrary number of quasiparticles (either end-A 
or end-i? ones), describing their braiding and fusion as 
well as the quasi-quantum double mathematical struc- 
ture. 



V. EXAMPLES 

As will become clear through this section, the gauge 
charge and gauge flux quasiparticles in our theories do 
not behave in the same way regarding the phenomena 
of symmetry fractionalization. We will therefore focus 
on the gauge fluxes, which have non-trivial properties. 
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while leaving the case of gauge charges, which behave 
trivially, to subsection IV Fl 

At the end of this section, we will also discuss the de- 
generacy of states with multiple quasiparticle pairs by 
considering the dualization of the global symmetry SG. 



A. Hilbert space for a pair of gauge fluxes 

We will now analyze the properties of excitations in 
our exactly solvable models, as announced in Sec. [Til We 
will later study in detail examples where the groups SG 
and GG are products of ^2- In general, we study a pair 
of gauge fluxes positioned at the ends A and B (which 
will also label the lattice sites there) of a string F, created 
by the action of a ribbon operator 



\hy,g]UA) 



= 4,(A)F(''"^9)(F)|gs), (102) 



where we introduced more succinct notation. These 
projected states span C{A,B,hy). According to the 
general definition of ribbon operators, hy G GG and 
g ^ G = SG X GG, and we label the unique factors of the 
latter element as g = hg ■ g, with hg e GG, g € SG. The 
value oiuB is automatically fixed by the action of ribbon 
operator ub = g^^ -ua, with ua,ub S SG. We will focus 
on a pair of gauge fluxes, so that hy is chosen as a fixed 
non-trivial element of GG. (In Z2 gauge theory, such a 
flux quasiparticle is called "vison".) 

Recalling the discussion in Section IIV CI 
£i(ft„,g)(^)P„^(A) = Pg.uM)D(h^,g){A) because 
the local operators transform these elements, e.g. 
UA g ■ UA with g — h ■ g, h £ GG, g G SG. The action 
of local operators on the Hilbert space of the gauge flux 
pair having flux hy is therefore given by: 

D[gi]{A) \hy,g;uA) = Cg^^{hy,h;;^)chAgi,g) \hy,gi ■ g;gi- 

B[gi](B) \hy,g;uA) = c/ijgf^ ■ 5,51) \hy,g^'^ ■ g;uA) , 

(103) 

where we remind that gi — hi-gi, hi £ GG, gi £ SG, and 
since the flux hy is flxed we introduced the shorthand no- 
tation D[g]{A) = i?(;,-i,^)(^) and D^g^iB) = D^h^^g){B). 
Our analysis of the examples will use Eqs. (|103p to ex- 
plicitly construct the matrices of local operators in this 
basis. 



B. Construction of local symmetry operators in 
symmetry-fractionalized models 

Let us now, from a general viewpoint, consider the 
possibility of symmetry fractionalization in our models. 
According to Eq. (|20)) . we seek to factorize the global 
symmetry transformation U{g £ SG), when acting in 
the flux-quasiparticle Hilbert space C{A, B, hy), into two 
local factors: 



First of all, we observe that the global symmetry trans- 
formation g acts by transforming elements on each lattice 
site i by sending Ui — > Ui -g^^, while the operators D per- 
form a similar operation on a single site, i.e. locally, for 
example D(^h^g-){A) : Uia ^ g ■ UiA, where g = hg ■ g, 
hg £ GG, g £ SG. (Recall that we deal with Abelian 
groups G here.) This is reasonable, since for any local 
operator in C{A, B, hy), such as the tentative Ug{G), we 
expect that it is representable in terms of the local oper- 
ators D{G). 

We therefore see that whenever symmetry fractional- 
ization occurs, it should be possible to flnd a local phase ip 
such that U~g{C) = e''P^^^--'3-^^)D^,,^g-i){C), where g £ SG 
and C ^ A,B. (Note that although the phase ^p{G) can 
depend on local variables at C, e.g. on uc, the func- 
tional form cannot depend on the position C , since the 
local symmetry operation cannot depend on the spatial 
position at which it is applied.) We will now show under 
which conditions (i.e. for what kind of 3-cocycle w) it is 
possible to flnd such a phase ip. 

We can effectively use the demand on the global sym- 
metry transformation U{a) which says that U{'j) com- 
mutes with the ribbon operator creating the quasipar- 
ticles; more precisely, only ua,ub of the projected ba- 
sis vectors (Eq. (I102p ) in C{A,B,hy) are transformed 
by Uc (T^^ ■ uq. Tentatively writing Ua{A) = 
UaiB) = D(^^,<,-i)(B) (note that the 

gauge fluxes at end- A and end-B are hy^ and hy respec- 
tively), with a £ SG, we get 

Ua{A)Ucr{B) \hy,g;uA,UB) = 
<^h^,a-\gC^-i{hy,h;;^) \hy,g;a^'^ ■UA,cr^^ ■ ub) , (105) 

where we have introduced the 1-cocycle 



UA] 



(-x,y,z = -^^T^^ x,y,z£G. (106) 
cx{y,z) 



Obviously, with this choice of U{G), the quasiparticle 
state is properly transformed only up to a phase, but 
we can proceed to absorb the resulting phase by a non- 
trivial <p(C). We will need some useful properties of the 
introduced 1-cocycle: 



^x,y,z-w — ^x,y,z ^x.y.w, 

^x,y,z = ^y,z,x = ^z.x.y — ^y.x.z ~ ^x,z,y ~ ^z,y,x^ 
-1 



(107) 



x.y.z 



x,y,z- 



-1 ■ 



Using the fact that g — hg ■ uau^^ and the listed prop- 
erties of e, we can define a valid phase by: 



UM) = V^^MKh^^h,a-^,UADih,a-^){A) (108) 



U^B) = Vca-i{h,h-^)e^^^- 



)iB). 



U{g) = Ug{A) -UgiB), g£ SG. 



(104) 



Note that the h in UaiA){Ua{B)) is defined to be the 
gauge fiux at the end-A(end-i?), which can be measured 
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locally. This leads to 

Ua{A)Ua{B) \hy,g;uA,UB) 

=(^h^,a~\h,U{a) \hv,g;uA,UB) ■ (109) 

Symmetry fractionalization f/(cr) = Ua{A)U„{B) can 
therefore occur if 

eh^,a,h, - 1, yK,hg eGG,ae SG. (110) 





-1 


x,y,z 






(1,*),(1,*),(1,*) 






(*,1),(*,1),(*,1) 




(1,*),(*,1),(*,1) 



TABLE I: Inequivalent basic 3-cocycles tj*^' for G = GG x 
SG ^ Z2 X Z2 which label classes in H^{G,U{1)) = Z|. A 
cocycle uj{x,y, z) has value 1, except when the group elements 
x,y, z £ G take the special values shown in the table, in which 
case uj{x,y,z) is equal to —1. The element notation x = 
(J2) signifies the direct product factorization, i.e. gi G GG 
and g2 £ SG, and the elements of Z2 are 0, 1 (additive group 
action). The star symbol * stands for "any value of element". 



It is not physical to absorb this phase factor into the 
local operators Ua{C). The reason is that the phase de- 
pends on the element value hg carried by the quasipar- 
ticles (recall the definition of ribbon operators Eq. ([73|: 
hg = hifi ■ • ■ hjsf^isf-i £ GG for a ribbon on lattice sites 
0, . . . , iV), while on the other hand the symmetry trans- 
formation should not depend on such specific proper- 
ties of the quasiparticles. The presented argument at 
least gives an indication when symmetry fractionalization 
should be impossible. We will see that for GG — Z2X Z2 
and SG = Z2 there is one 3-cocycle, w^^^^) (see TableHIII . 
for which Eq. (|110p is violated, in accordance with our 
claim about a single Z2 index beyond SFC (see example 
"3" at the end of section Hi C|) . For that example, we will 
additionally show that global symmetry exchanges the 
two species of excitations, thus giving a strong physical 
argument against the possibility of representing symme- 
try by local operators. 

In Appendix IfI we show that this construction of frac- 
tionalized symmetry operators is generally correct for 
multiple quasiparticles in the examples below. 

Returning to the present single ribbon case, let us for 
a moment assume that symmetry fractionalization is ex- 
plicitly possible in a given model, i.e. Eq. ()110p is sat- 
isfied. We can then give general formulas for the trans- 
formation of a single quasiparticle at position G under 
fractionalized symmetry transformations U^{G), where 
again G = A, B is either end of ribbon. The fractional- 
ized symmetry transformations physically only need to be 
projective, as discussed in section Hi Bl In the examples of 
the following subsections, these general expressions will 
be used to check whether the symmetry group generators 
obey the SFC relations shown at the end of section Hi Bl 
(The results will be corroborated by explicit construc- 
tions of Ua{C)). Because our examples will be based on 
the group Z2 and on direct products of the group Z2, 
we can make some simplifications. Namely, we can set 
= hy. Further, considering arbitrary symmetry op- 
erations go, 1,52 G SG, we have that 5o = 5i = 52 = 
with 1 the group identity element. Finally, for such 
groups the inequivalent cocycles lj can be chosen to only 
take values ±1, see Tables IIIIIII Now, using Eq. (fTOS]) 
and (|108|) . as well as properties of the 2-cocyclc (see 
Eqs. (|87a|) - (l87dp ). it is easy to show that in the Hilbert 



space of a hy flux pair: 

UUC)UgAG) = euM^UgACWUC) (HI) 
UgaiAf = Ch^{go,go)cgo(hy,hy) (112) 
D[g^]{Af ^ChS9o,~9o}, (113) 

where the fluxes are positioned at ribbon ends G — A, B. 
In Eq. (|113p we show the result for the standard local op- 
erator D, to contrast it with the fractionalized symmetry 
operator in Eq. (|112p (in both cases the operator on the 
right-hand side is just the identity operator). 

Equipped with this formalism, we can proceed to solve 
the instructive examples introduced in sections Hi Bl III CI 



C. The simplest example with symmetry 
fractionalization (PSG): Vison pair in GG = Z2 and 

SG = Z2 

Let us from now on use the additive notation for our 
groups, i.e. any element g G Z2 takes values g = 0,1, 
with the identity element, and the group product be- 
comes addition gi • (72 = ffi + 52 (mod 2). 

The different classes of phases described by our mod- 
els are classified by the inequivalent 3-cocycles lo. An 
explicit expression for all representatives of inequivalent 
3-cocycles in case where the group is a product of Z„ 
factors is given in Ref. [s^, and the result for the case 
G — Z2 X Z2 is shown in Table HI This means that there 
are three indices pi,p2,P3 € Z2, labeling the elements 
of H^{SG X GG,U{1)) = Z|. The trivial phase has 
Pi — P2 = P3 = 0, and can be described by the triv- 
ial cocycle oj{x, y, z) — 1, Va;, y,z ^ G. When an index pi 
is 1, the representative 3-cocycle w G H^(SGxGG,U{l)) 
of that phase is chosen to satisfy the property described 
in the definition of w*^-^); when the index pi is 0, this in- 
dex does not affect the value of cocycle. In other words, 
the elementary cocycle properties in Table |T] gener- 
ate all representative cocycles that are elements of H^: 
The three basic 3-cocycles lo^^^ generate a total of 2^ 
3-cocycles in H^{SG x GG,U{1)) = Zj, each of them 
leading to a physically different model, as stated at the 
end of section Hi CI 
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-^[(0,1)] 
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-^[{1,1)1 
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^Py Px 




iPy 


Tx 



TABLE II: The values of local operators and fractionalized symmetries for GG = Z2, SG = Z2 for (the only possible) vison 
pair hy = (1, 0) = 1 G GG. D and D' are shorthand for D{A) and D{B), respectively. The cocycle ^(12) is the only one leading 
to non-trivial projective realization of local symmetry, Ua{C)^ — —1, and is therefore identified as the cocycle in H'^{G, (7(1)) 
which generates the non-trivial symmetry fractionalization classes, i.e. the index pi2 labels the SFC{Z2 x Z2,U{1)) — Z2- 
(The action of the matrices in the vison pair Hilbert space is given in the text.) 



Let us now consider a pair of visons ft-t, = 1 e GG = 
(1,0) e G, which are the only flux quasiparticles for this 
group G. 

Since the value of is fixed, we will omit it from 
the label of local operators, i.e. we will use the labels 
^\g={h,(j)\ with h £ GG^g G SG (this notation was intro- 
duced after Eq. (|103p ). Let us introduce Pauli matrices 
Mi T Pii'''i acting in the Hilbert space of the visons on the 
basis \hv = \,g]UA) = \{hg,g)]UA) from Eq. (|102p . with 
the usual g e G = {hg e GG,g £ SG). The action of 
the matrices is naturally defined through setting /i^ as 
(-l)s, as and finally as (-1)"-^. The local 

operators, Eq. (|103l) . together with the global and frac- 
tionalized symmetry operators, Eq. (I108|) . are presented 
in Table [n] for all basic 3-cocycles uj'--'\ 

We can see that the global symmetry U{<t = 1 E SG), 
which sends uc uc -t- 1 is in this basis equal to , and 
as expected commutes with all local operators, for any 
choice of cocycle. 

The phase in the definition of fractionalized symmetry, 
Eq. (|108p . is for all cocycles trivial (note that the cocycle 
^x,y.z contains one factor of uj for each permutation of 
x,y,z). The fractionalized symmetry is then just equal 
to the local operator Ua{G) — D[(o,i)](C), and clearly 
for each cocycle it is true that [/(cr) = Uo-{A)Ua{B) = 
-D[(o,i)](^)-D[(o,i)](5) = 

There is a non-trivial symmetry-fractionalized phase 
in case of uj^^'^\ since Ua{C)^ — —1. (One can also con- 
firm this from the general expression Eq. (|112p .) That 
means precisely that the index pi2 labels a non-trivial 
realization of symmetry fractionalization, where the lo- 
cal symmetry operation on a single vison at position G 
is projectively reahzed Ua{G)^ = {—1)'p^^1, exactly as 
claimed in Eq. (j25p ). The index pi2 is obviously the only 
index classifying the SFC, SFC{Z2 x Z2, (7(1)) = Z2, as 
presented in Eq. (1571) . 

We can now consider the influence of interplay between 
global symmetry and topological order on the physical 
system with excitations. Let us assume that the excita- 
tions are a vison pair at A,B. The topological nature of 
the state must be robust against arbitrary local pertur- 
bations, i.e. we have to consider adding arbitrary local 
perturbation terms to the Hamiltonian (Eq. ([57|). The 
local perturbations however also have to be symmetric, 
i.e. commute with the global symmetries. The set of all 
local perturbation terms is formed by the local operators, 



i.e. D = {D{A),D{B)}, so that 
Hp 

erturb = H + J2aaD", (114) 

a 

where the index a labels all operators in the set D, and Oq, 
are arbitrary coefficients. The operators D are by con- 
struction symmetric, which is easily explicitly checked 
by showing that their matrices commute with the global 
symmetry operation . For any fixed elementary cocycle 
in Table |ITJ one can see that there are no other matrices 
that commute with the perturbations (except the global 
symmetries of course). Further, all the local operators, 
for any fixed cocycle, commute with each other. There- 
fore, the algebra of conserved observables is trivial, and 
there are no degeneracies protected by symmetry. (This 
will change in the next examples.) 

Physically it is important to note that the eigenvalues 
of the local operator £'[(i.o)ggg] are actually the values of 
gauge charge of the vison excitations (note that they act 
on the hg element in the ribbon operators, where from 
Eq. ([73|: hg — hi^ ■ ■ ■ h^^N^i £ GG for a ribbon with 
inner edge on lattice sites 0, ...,iV). These operators 
are indeed the same on both visons, o)6GG] (^) = 
-D[(i^o)eGG](^)- The gauge charge value is an important 
quantity that will be non-trivial when the GG is enlarged 
in our third example. 

D. Example with symmetry protected degeneracy: 
Vison pair in GG = Z2 and SG = Z2 x Z2 

In this example, we again have only one type of vison, 
hy = 1 £ GG, but two Z2 global symmetry generators, 
a = (1, 0) £ SG = (0, 1, 0) e G and t = (0, 1) £ SG = 
(0,0,1) £ G. 

We use the elementary 3-cocycles for the group G = 
Z2 X Z2 X Z2 according to the results of Ref. [13, 
and present them in Table IIIIl This means that there 
are seven indices pi = 0, 1, labeling the elements of 
H^{G, Ull)) — Z\. The seven elementary 3-cocycles w^^^ 
generate a total of 2'' cocycles in H^{G, (7(1)) = Z|, each 
of them leading to a physically different model, as stated 
at the end of section Hi CI However, the physical proper- 
ties of the models will be very different depending on the 
particular choice of SG and GG, even though G is the 
same in this and the next example. 
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TABLE III: Inequivalent elementary 3-cocycles o;'^' for G = 
GG X SG = X ^2 X which label classes in H\G, [/(I)) = 
Zj. A cocycle cj(a;,y,z) has value 1, except when the group 
elements x,y,z £ G take the special values shown in the table, 
in which case oj{x, y,z) is equal to —1. The element notation 
^ = (<?i)S2,33) signifies the direct product factorization, and 
so in the example GG = Z2, SG = Z2 x Z2, gi £ GG, while in 
the example GG = ^2X^2, SG ^ Z2, ga £ SG. The elements 
of Z2 are 0,1 (additive group action). The star symbol * 
stands for "any value of element". 



C/(ct)2 = C/(t)2 = 1 and C/(a)[/(r) = U{T)U{cr), the 
fractionalized symmetry can be realized projectively: 

• For cocycle uj^-^^^ (which "mixes" the GG and the 
first Z2 of SG), one finds the only non-trivial case 
C/,(A)2 = (7,(B)2 = -1. 

• For cocycle w*-^'^-' (which "mixes" the GG and the 
second Z2 of SG), one finds the only non-trivial 

\2 _ TT { n^2 



TABLE IV: Relevant local operators and fractionalized sym- 
metries for GG — Z2, SG — Z2 X Z2, for (the only possible) 
vison pair fti, = (1, 0, 0) = 1 G GG located at A, B. The sym- 
metry generators are a — (1,0) G SG, t = (0, 1) G SG, and 
they act as (7(a) = and Uij) = p^- There is symmetry 
fractionalization (7(s) = 17s (j4))7s(i3), non-trivial (projective) 
for cocycles aj'^^', and aj'^^s)^ ^j^jgij "j^j^" qq ^j^-j^ ^.q 

f Table HII)) . The fractionalized symmetries differ (by a local 
operation) from local operators -D[o-], -Dfr] only for the o;^^^'^-' 
cocycle (see Eq. p08|) ). The definitions of all matrices in the 
vison pair Hilbert space are given in the text. 



After fixing it is straightforward to construct the 
matrices for all local operators by using Eqs. (jl03l) . 
as well as to construct the fractionalized symmetry 
operators using their definition, Eq. (|108p . for each 
cocycle in Table IIIIl We use the tensor product of 
Pauli matrices which are defined through setting as 
(-!)''«, as (^1)91, a, as (-l)9^ as (-1)"^% 
and pz as (— 1)"-*^ in our standard basis \hg,g;uA) = 
l^s X gi X g2;uA\ X with hg G GG, (31,32) G SG, 

[UA\,UA2) G SG. 

Direct inspection of the obtained matrices, which are 
presented in Table HVl reveals the results quoted in sec- 
tions III Bl III CI Namely, for all cocycles one indeed finds 
that both global symmetries, C/(cr) — and C/(r) = 
Px, SLie fractionalized: U{s) = Us{A)Us{B). Although 
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• For cocycle w^^^^^ (which "mixes" all three sub- 
groups), one finds the only non-trivial case 
U„{C)Ur{C) = --Ur{G)U„{C), for both C = A,B. 

One should note that only in the case w^^^'^) the local 
operators had to be modified to obtain the fractionalized 
symmetry operators, as clearly seen in the last column 
of Table IIVI We can directly confirm that these are local 
redefinitions by noticing that the operators (— l)"-*i — v^, 
(— 1)"-*^ — a,re local at A, while operators (— 1)"^^ = 
f^Pz, (—1)"^^ = [izVz are local at B. 

We can now consider the influence of interplay between 
global symmetry and topological order on the physical 
" system with excitations by repeating the analysis of local 
symmetric perturbations as in the previous example. (In 
'this example it is also easy to explicitly check that the 
local, symmetric perturbation terms D — \I)[^A),D(B)] 
indeed commute with the global symmetry operations , 

.Px-) 

We can now ask: What operators can be used to la- 
bel this state? Let us for concreteness consider the state 
having the cocycle w^^^'^^. The matrices that commute 
with the entire D algebra (which includes £'[geG] opera- 
tors not shown explicitly in Table ITV)) . and therefore with 
the arbitrarily locally perturbed Hamiltonian (Eq. (|114p , 
are: 

• Px, i^x, which arc just the global symmetries. 

• Tx, which measures the gauge charge, as will be 
explained further below. 

• a^Vz, PxTzPz, cr^Px^z, Pxtyz^xPz, which are just the 
fractionalized symmetries Ur{B), Ua{B), Ur{A), 
and Ua{A), respectively. (Note that these oper- 
ators are not symmetric themselves, as a conse- 
quence of U„{C)Ur{C) = ~Ur{C)U„{C).) 

In the algebra formed from these conserved operators, 
there are exactly two disjunct pairs, Ur{B),Ua{B) and 
Ur{A),UaiA), respectively, which anticommute, as we 
already learned through the symmetry fractionalization 
for w^^^-^) . Each anticommuting pair, acting on one of the 
visons, forces a two-fold degeneracy on the state. (We 
checked this degeneracy numerically by considering an 
arbitrary perturbation term from Eq. pi4p .) We arrive 
at the physical signature of this phase: 

The symmetry protects a 2-fold degeneracy per 
vison of the pair in the P123 = 1 phase of GG = Z2, 
SG — Z2 X Z2 model. The degeneracy due to a vison 
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at C is labeled by the fractionalized symmetry operators 

U,{C),Ur{C). 

Let us now discuss the gauge charge observables t^. 
These arc just the local operators £'[iggg]j listed in 
Table HVl Physically, the operators D[i^qq-^{B) appear 
in the braiding matrix Ti^^j, Eq. (llOOp . describing the 
braiding of two visons at A-ends of two ribbons which 
share a B-end, see Fig. [1^1 Section llV DI and Appendix IB] 

Finally, let us emphasize that, at least in the case when 
G is a product of simple Z2 factors, there is no need to 
explicitly construct all the matrices as was done here. 
One can directly resort to Eqs. pTT|) - ([Tn| to test 

whether a cocycle allows SF, and then whether this SF is 
non-trivial. For this purpose, we note that the cocycle e is 
easy to use since it contains 6 factors of uj{x,y,z), such 
that each permutation of x,y,z appears exactly once; 
importantly, when G has only Z2_factors, the generating 
3-cocycles in H^{G,U{1)) (Ref54) can only take values 
±1, as is the case in Tables HI IIIII 



E. The simplest example beyond symmetry 
fractionalization: Vison pairs in GG = Z2 x Z2 and 

SG = Z2 

In this example there are two fundamental vison types. 
We will consider in parallel a pair oi ~ (1,0) G GG = 
(1,0,0) S G visons and a pair of hy — (0,1) G GG = 
(0, 1,0) e G visons. 

We again construct the matrices of all local operators 
by using Eqs. (|103l) . Before proceeding, one should no- 
tice that the condition stated in Eq. (|110|) is violated for 
the cocycle w*^^^^^: only in that case does eh^,cr,h — ~1 
become possible since hy and hg can differ while both be- 
ing non-trivial elements of GG. According to Eq. (|109p . 
we then expect that in that case the symmetry fraction- 
alization will fail. 

Let us nevertheless continue with the explicit construc- 
tion of fractionalized symmetry operators. Using their 
definition, Eq. (|108|) . for each cocycle in Table |lll] includ- 
ing a;'-^^'^\ each fractionalized symmetry operator is just 
equal to its corresponding local operator D. We therefore 
avoid repetition in Table |Vl and one should keep in mind 
that Ua{C) = D[o-](G), where the only global symmetry 
generator is cr = 1 G SG. 

The Pauli matrices are here defined through setting 
as (-1)''»S as (-1)'*"", ct^ as (-1)^, and as (-1)'' 
in our standard basis 



ig,g;uA) = \hgi X hg2 X g;uA}, 
with {hgi,hg2) G GG, g G SG, ua G SG. 

Direct inspection reveals that the global symmetry, 
represented by U (tr) = px, is fractionalized for both vi- 
sons, i.e. U{(j) = U^^ [A)!!^" {B) = p^, except for the 
cocycle lu^^'^^\ In the case of this cocycle, we get 



u!^^=(^'")iA)u':^=(''"HB)^^,,Px 

U'jM)(^A)U'j^=('''\B)^r,Px 



(cocycle a;(^2=^)), 
(115) 
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TABLE V: Relevant local operators for GG = Z2 x Z2, SG ^ 
Z2: (top half) for a vison pair = (1,0,0) = (1,0) G GG 
located at A, B; (bottom half) for a vison pair hv — (0, 1, 0) = 
(0, 1) G GG located at A, B. The only global symmetry gen- 
erator is (T = (0,0,1) = 1 G SG, acting as U{o) = px. 
Except for there are fractionalized symmetry oper- 

ators UaiC) = Diai(C), with C = A,B, (see Eq. (fTUS)) ). 
such that U{o) — Ua{A)Uc{B). The fractionalized symme- 
tries are non-trivial (projective) for cocycles oj'-^''' and cj'^^-*, 
which "mix" GG with SG (Table |ITI1). Phase having 



,{123) 



beyond SEC, since for it 0'^^.^-''°^°'' {A)D^Jf^"'"> (B) = T^Px 

-^IctI (^)^[ct] '^'"'^ (-^) ~ f^zPx, and no local redefinition 
of _D[o.] can remove the Tz,Pz factors. The definitions of all 
matrices are given in the text. 



-)h^ = (010). 



for the vison pairs — (1,0) and — (0,1), respec- 
tively. We will further discuss this below. 

Considering all the cocycles [except w'^^'^-'), we find 
phases with non-trivial SF, i.e. where the fractionalized 
symmetry is realized projectively: 

• For cocycle cj^^'^'' (which "mixes" the first Z2 in 
GG with SG), one finds the only non-trivial case 
f/>=(i^o)(A)2 = C/>=(i'°)(S)2 = -1. 

• For cocycle a;^^'^' (which "mixes" the second Z2 in 
GG with SG), one finds the only non-trivial case 



1. 



=(0'i)(A)2 = C/.""=(°'^)(B)2 

These results match the claims in section III Bl There 
are obviously two symmetry-fractionalization Z2 indices, 
P13 and P23, in accordance with the general results: 
H^{SG, GG) = Z| and, for this group, SFC{SG, GG) = 
H^{SG,GG) (see section HTC]). 

We can try to use our explicit matrix expressions to 
"force" symmetry fractionalization in the case of cocycle 
£j(i23)^ This attempt will fail, as expected from general 
arguments above, and therefore this state is beyond SEC, 
with the index P123 labehng the EXTRA(S'G, GG) class. 
What we would need to achieve is the removal of fi^ and 
T2 factors in Eq. pi5p . Inspecting the operators 



(-1)"^ = 



Pz 



(116) 
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the first being local at A, and the second local at B, we 
see that there is no hope in manipulating the /i and r 
subspaces. Physically, /Zj and Tz measure the GG de- 
grees of freedom, and these are connected to the vison 
gauge charge which is non-local. The local operators 
^[hv] (which we will use to measure the gauge charge 
shortly) by definition change this group element (i.e., r^) 
and cannot provide us with Tz, Hz- 

We next ask: which operators can be used to physically 
label this state? Let us consider the especially interesting 
cocycle uj^^'^^\ As explained in the previous example, we 
need to find matrices that commute with the entire D 
algebra (which includes all D^g^Q] operators). These are: 

• Pa;, which is just the global symmetry. 

• ei, 62, which measure the two Z2 gauge charges in 
GG, as we explain further below. For the = 
(1,0) vison pair, ei = r^, 62 = lJ-x<^zPz, while for 
hy = (0, 1) vison pair, ei = r^azPz, £2 = ^-x- 

In the algebra of these conserved operators, exactly one 
pair anticommutes: the global symmetry anticommutes 
with the gauge charge of Z2 complementary to the visons' 
flux hy. For example, considering a vison pair with flux 
set hy hy — (1,0) G GG, the second Z2 gauge charge 
of GG (62) anticommutes with the global symmetry. A 
signature of this phase is then: 

The symmetry protects a 2-fold degeneracy of 
a single vison pair in the P123 = 1 phase of GG = 
Z2 X Z2, 5G = Z2 model. 

This result for the cocycle uj^^'^'^^ reveals its physical in- 
terpretation. To understand it better, we need to explic- 
itly consider the consequences of anticommutation be- 
tween symmetry and gauge charge. Let us proceed with 
that goal, and leave for later the detailed explanation of 
why ei, 62 can indeed be understood as gauge charge 
operators. 

The quasiparticles of this theory are labeled by the 
values of flux and charge, and we consider the elemen- 
tary ones: 1) Flux particle mi ~ 1[to2 = 1] has flux 
in the first [second] Z2 of GG, which is set by the val- 
ues of hy = (*,0)[/i« = (0,*)], and no gauge charge; 
2) Charge particle 61 = 1[62 = 1] has charge in the 
first [second] Z2 of GG, which is set by the values of 
hg = (*, Q)[hg = (0, *)], and has no flux. A general quasi- 
particle will have arbitrary values (0 or 1) for these four 
numbers, so we label it (61,62,7711,7712). For instance, 
(1,1,0,0) is particle bound state of charge and flux in 
the first Z2 of GG. When symmetry anticommutes with 
6, that means it switches the two eigenvalues of e, and 
therefore switches the particles' e number between triv- 
ial (0) and non-trivial (1). As we established, symmetry 
anticommutes with ei when the particle pair in question 
is m2, while it anticommutes with 62 when the particle 
pair is 7711. This leads to the following transformation of 
elementary quasiparticles under the action of symmetry: 

(61,7711,62,7712) ^ (61 +7772,7771,62 -|- 777l , 7772) , (117) 



where (mod 2) algebra is understood. As expected, the 
symmetry changes the charges of particles, but we can 
further assign a deeper meaning to this transformation. 
Namely, in Z2 theories the charge and flux are physi- 
cally equivalent: they are dual to each other, and we 
can rename them at will without changing the physics 
(e.g. the braiding properties) of particles. We there- 
fore assign new fluxes 777 and charges e according to the 
rules: (61,777.1,62,7772) = (61+7712,7771,777,2,62 + 7711). Note 
that the quasiparticle statistics is indeed conserved, i.e. 
bosonic except between ei and 777^ when it is fermionic. 
This leads to the elegant transformation rule: 

(ei, 7711, 62, 7772) A (62, 7772, 61, 7771 ), (118) 

so that the global symmetry transformation only ex- 
changes the particle types 1 and 2! (We note that this 
physical interpretation might carry on to higher dimen- 
sions. For instance in 3d, charge excitations are point- 
like while fluxes are string-like, so they cannot be inter- 
changed among each other; however it is still well deflned 
to change their type, e.g. 1 <-> 2.) 

We have therefore found the physical nature of this 
state: The global symmetry operation exchanges 
the two quasiparticle types in P123 = 1 phase of 
GG = Z2 X Z2, SG = Z2 model. This shov^rs that 
symmetry performs a non-local transformation on 
the excitation pair, and the state is beyond sym- 
metry fractionalization classification. 

Let us now discuss in detail the gauge charge ob- 
servables ei, 62, as promised above. We focus on a 
particular vison pair hy. As in the previous example, 
one naively expects the gauge charge observables to just 
equal the appropriate local operators: -P[(i,o)gGG=(i,o,o)] 
and i?[(o,i)6GG=(o,i,o)]j which are listed in Table Ivl 
Focusing on the cocycle w^^^'^-', we however see that 
these local operators differ at the two quasiparticles, e.g. 
-C[(i,o.o)](^) 7^ D[{ifi,o)\{B), which would mean that mea- 
suring the gauge charge of a vison and anti- vison in the 
pair would give differing answers. This is physically 
wrong, but it is easy to resolve the problem. Namely, 
the physical operators 61,62 differ from the correspond- 
ing local operators Z? by a simple local operation. Using 
the local transformations in Eq. (jll6p . one immediately 
obtains that 

• For vison pair hy = (1, 0): 

ei = ■D[(i,o,o)](^) = ^[(1,0,0)] (-S) = Tx 

62 - (-l)"-^7^[(0,l,0)](^) = (-l)"-i?[(0.1.0)](S) = Px^zPz. 

(119) 

• For vison pair hy — (0, 1): 

61 ^ (-l)"^i^[(i,o.o)](^) = (-l)""C[(i.o,o)](B) = TxazPz 

(120) 

62 = -D[(o,i,o)](^) = ^[(0,1,0)] (-8) = Px- 
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These 61,62 operators exactly match the ones we found 
to commute with the entire local D algebra. 

We already know that the local operators D^^^] [B) ap- 
pear in the braiding matrix TZ^c -^1- POOP - describ- 
ing the braiding of two visons at A-ends of two ribbons 
which share a i?-end, see Fig. [161 Physically, we expect 
these local operators to equal the gauge charge operators 
and thereby transparently provide the braiding rules of 
anyons in Z2 topologically ordered theory. A question 
that arises then is: If we redefine the local operators 
so that they become exactly equal to the gauge charge 
operators (following Eqs. (|119l) . (|120p '). do the braiding 
properties change? The answer is negative. Namely, if 
we choose to braid to visons of same type h\ = = hy, 
the braiding operator (Eq. (|100|) ) will apply the same 
D(^h^,h^){B) on both visons, and these local operators 
are not even redefined. When we braid two differing vi- 
son types, e.g., hi = (1,0), = (0, 1), the braiding will 
apply Z?(;ji ;j2) and -D(/i2 ^ji), both of which are redefined 
by (-1)"^" according to Eqs. (ITTO| . p^ . As pointed 
out earlier, our analysis of vison Hilbert space revealed 
that such braiding operations should be considered after 
choosing, and keeping fixed, some value of ub- Since the 
two ribbons share the end at B, the (— 1)"^ operation ap- 
plied on both braided quasiparticles cancels. Therefore, 
we can freely redefine the considered local operators to 
be equal to the physical gauge charge operators ei, 62. 



F. Gauge charges 

In the previous subsections we considered in detail 
pairs of gauge fiuxes, i.e. visons. The other fundamental 
type of excitations is a pair of gauge charge excitations, 
which is created by the ribbon operator: 

\g)^F('^-^-^\r)\gs}, (121) 

i.e., the special case of the flux h e GG being trivial. As 
always, g E G = SG x GG uniquely factorizes as g = 
hg ■ g, with hg G GG being related to the gauge charge of 
the excitations, and g G SG. (Recall again the definition 
of ribbon operators Eq. ([73)) : hg = hi^o ■ ■ ■ h^M-i £ GG 
for a ribbon on lattice sites 0, . . . , N .) 

Recall that in this paper we use "canonical" 3-cocycles 
w, meaning that uj{gi, g2, gs) = 1 if any of 31,32,^3 is 
equal to 1. Specifically, the elementary cocycles for a 
group G which is a product of Z„ factors are written in 
Ref|54, and are all canonical. 

Since = 1 for gauge charges, it then follows imme- 
diately from Eqs. ([TTT|) . ([TT^ . ([TU5| that for these exci- 
tations the symmetry fractionalization is always trivial. 
Further, Eq. (j99p shows that the braiding matrix between 
various gauge charges is also trivial in our models. 

As a particular demonstration, let us consider the ex- 
amples from the previous subsections. Eq. (jl03p reveals 
that the action of local operators D in our basis of the 
excitation pair Hilbert space C{A, B) does not involve 



any phase factors. It is as if the chosen cocycle is always 
the trivial one, uj — 1. To avoid repetition in writing ex- 
plicitly the local operators in the present case of a gauge 
charge pair, we refer the reader to just look at the case 
of: foj. ^-^^ GG = Z2, SG = Z2 model; a; (2) for the 
GG = Z2, SG = Z2 y. Z2 model; and finally w^^) for the 
GG ^ Z2-X Z2, SG = Z2 model. The lack of non-trivial 
phase factors in the action of local operators leads to a 
trivial algebra: Repeating the analysis from above, one 
finds that in all three examples all the D operators com- 
mute, and the only matrices commuting with them are 
the global symmetry operations, which of course com- 
mute amongst each other in these direct product groups. 
Therefore, there are no symmetry protected degeneracies 
for the gauge charges in our present examples, in sharp 
contrast to the case of fluxes. 



G. Multiple vison pairs and dualization 

Let us briefly consider a system with multiple vison 
pairs. For simplicity, we focus on the case where all these 
quasiparticles are of the same type. Since the calculation 
even for this case becomes very complicated, we will try 
to give general arguments and speculate about the sym- 
metry protected degeneracy of a general state with N 
pairs of visons. 

One way to analyze this situation is to consider the du- 
alization of SG. By this we mean the standard replace- 
ment of lattice site degrees of freedom ut by edge degrees 
of freedom Uij = Ui ■ uj^ , as we also mentioned in sec- 
tion [TTlB] The total group becomes a pure gauge group, 
G = GG X SG = G, where SG denotes the gauge group 
SG. The visons in the G model are mapped to gauge flux 
particles in the G model. Let us denote the protected de- 
generacy of N — particle-pair state by VDEGg(A^) and 
VDEGg(A^) for the models with SG and the dualized 

SG, respectively. 

One should note that this is a many-to-one mapping, 
in the sense that \SG\ states obtained by multiplying all 
Ui by an arbitrary element s G SG (in the G theory) 
get mapped to a single state (in the dualized G theory). 
In general therefore, we expect that the state with vi- 
sons can only have the same or smaller degeneracy upon 
dualization to G: 

VDEGG(Af) > VDEGq(7V). (122) 

Particularly, the mentioned many-to-one nature of the 
mapping indicates that the change in degeneracy (from 
G to G) might involve the factor \SG\. 

In fact, we can recall that the Hilbert space C{A, B) 
of an excitation pair in the SG description (section [V^ 
had to be specified by keeping track of value of the ele- 
ment ua, or Ub, which belong to SG. In the most gen- 
eral description of the 2A^ visons, one connects all their 
ribbons to a common point on the lattice xq , which con- 
tains no excitation (see discussion in section II VP I) . We 
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then have to keep track of Uxo, which takes \SG\ differ- 
ent values. Since this degree of freedom becomes obsolete 
upon dualization to G, one is again lead to the assump- 
tion that the degeneracy of the 2A^-vison state in the 
dualized theory (G) is reduced by a factor of \SG\. 

Let us test these assumptions on the results of Sec- 
tion |Vl In our first example SG = Z2, GG = Z2, a 
single vison pair had no symmetry protected degeneracy. 
Upon dualization, this model becomes the gauge theory 
G = Z2 X Z2. This is an Abelian theory, which implies 
that its excitations are anyons with quantum dimension 
equal to 1. (The quantum dimension is defined as the 
degeneracy per excitation in the limit of infinite num- 
ber of excitations.) In this case, there is no degeneracy. 
Therefore, the inequality (|122p is saturated, and there is 
no factor \SG\ = 2. 

The Abelian case, having non-degenerate anyon states, 
might be a too special case for our general considera- 
tions. Let us instead focus on the examples with G — Z2, 
and the especially interesting cocycle u;^^^^\ In those 
cases we found that a single vison pair has degeneracy 
4 and 2, in the SG = Z2 x Z2 and SG = Z2 exam- 
ples, respectively. The topological order described upon 
dualization of either of these example groups becomes 
very interesting due to the presence of the cocycle ui'-^^^h 
Namely, the topological order realized in that way can 
be described by the non- Abelian theory having the gauge 
group GnA — D4 and no twist by a cocycle^. It is known 
that the D4 theory has no degeneracy for a single exci- 
tation pair, i.e. VDEGd^ (1)=1. 

For a single excitation pair and the case where the 
topological order in the dual theory is described by 
GnA = -D4, we therefore established the anticipated 

VDEGg(A^) = |S'G|-VDEGg ^{N), (for G„a = D4„N 

(123) 

Let us proceed to the case of multiple pairs of parti- 
cles, still all being of the same type. The excitations of 
the D4 gauge theory have quantum dimension 2, there- 
fore for large N one has the scaling VDEG_d4(A^ :s> 
1) ^ 4^. Since for = 1 we already discussed that 
VDEGd4(1) ~ 4°, we expect that an interpolation for- 
mula VDEGr,, (A) - 4(^-1) should hold for aU A^. Com- 
bining this expectation with the result in Eq. (I123P leads 
us to conjecture: 

VDEGg(A^) ^ |S'G|-VDEGg ^(A^) (non-Abelian G„a, 

In this conjecture we also anticipate that the key property 
which ensured Eq. (|123p is that D4 is non-Abelian, since 
in the previous paragraphs we established that Eq. (|123p 
does not hold (even though A^ = 1) when G described an 
Abelian topological order (i.e. Z2 x Z2). 

We leave further study of these questions for future 
research. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper, we present a classification of topologi- 
cally ordered phases in the presence of an on-site global 
symmetry, equipped with local bosonic exactly solvable 
models. The discussion in Sec lII CI and the solutions 
of our models in some examples allow us to reveal the 
physical meaning of this classification, which is given 
in Eq. ipTj) and Eg. (15^ . Basically, the classification 
H^{SG X GG, [/(I)) is a direct product of many finite 
abelian groups, and each group has a clear physical mean- 
ing. Among them H^{SG,Uil)) is the index for SPT 
phases, H^{GG,U{1)) is the index for Dijkgraaf-Witten 
topological orders, SET{SG,GG) in Eq.§^ is the in- 
dex labeling the possible interplays between the global 
symmetry and the topological order (symmetry enriched 
topological phases). 

SET{SG, GG) can be further understood. It is a di- 
rect product of two abelian groups, SFG{SG,GG) and 
EXTRA{SG,GG). SFG{SG,GG) is labeling different 
possible ways to fractionalize the global symmetry by 
the topological order, while EXTRA{SG, GG) is label- 
ing the phenomena that the global symmetry transfor- 
mation could interchange quasiparticle species, which is 
beyond the symmetry fractionalization scheme. 

Some measurable consequences of the symmetry en- 
riched topological phases are discussed. In particular, 
we show that under certain conditions, symmetry en- 
riched topological phases have non-trivial symmetry pro- 
tected degeneracy for excited states. For instance, in one 
striking example beyond the symmetry fractionalization 
scheme, we show that symmetry protects a 2-fold de- 
generacy for a single pair of gauge flux quasiparticles, 
which cannot be locally associated with either quasipar- 
ticle. This degeneracy of excited states can be used to 
detect the symmetry enriched topological phases in nu- 
-'^Inerics/experiments. However, we leave the most gen- 
eral diagnosis of symmetry enriched topological phases 
for future investigation. In any case, the exactly solv- 
able models constructed in the current work can be very 
useful tools for this purpose. 

The SPT phases are known to host gapless edge states, 
topologically protected by the global symmetry. We have 
not studied the possible gapless edge states due to SET 
topological indices in this paper. We also leave this issue 
as a subject of future investigations. 

Below we discuss the generalization and limitations of 
our classification. 

>!)■ 

A. Generalization to higher dimensions 

Although we have been focusing on 2+1 dimensions, 
where the classification is given by H^{SG x GG, [/(!)), 
this can be easily generalized to H'^^^{SG x GG,U{1)) 
in general d-|-l dimensions {d > 2). One way to un- 
derstand this generalization is to dualize GG to be part 
of global symmetry, after which the system has an on- 
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site global symmetry SG = SG x GG. According to 
Reffigl. the SPT phases in d+1 dimensions are classi- 
fied by H'^'^^{SG^ U{\)), and it is natural to expect that 
these phases may still different before the duality. (In 
fact, this point of view is not completely correct, which 
we will discuss shortly in Sec lVI Cl ) 

The higher dimension classification is also equipped 
with local bosonic exactly solvable models. For instance, 
in 3+ld, one can also construct the exactly solvable mod- 
els on a 3-dimensional lattice which triangulates the 3- 
dimensional space. The GG degrees of freedom live on 
the edges and the SG degrees of freedom live on the 
vertices. Similar to the 2-|-ld case, the Hamiltonian of 
the model is a sum of local projectors. Theorem-1 and 
Theorem-2 (see Sec lIII A"2|) in 3+ld dictate that these 
local projectors mutually commute and the model is ex- 
actly solvable. We leave the complete solution of such 
3-t-ld models (i.e. the excited states) as a subject of fu- 
ture investigations. 

However, at this moment it is still possible to explore 
the physical meaning of the 3-|-ld classification using 
H'^{SG X GG, f/(l)). A Kiinneth expansion immediately 
gives: 

H'^iSG X GG, [/(I)) =H^{SG, C/(l)) x ij4(GG, C/(l)) 
xSET^{SG,GG) (125) 

Here clearly, the index H'^{SG,U{1)) is labeling the 
SPT phases, while H^{GG,U{1)) is labehng the differ- 
ent 3+ld topological orders described by the gauge group 
GG. Note that H^{GG, U{1)) is labeling the direct gen- 
eralization of the Dijkgraaf-Witten discrete gauge theo- 
ries in 3+1 dimension. Further, SET^{SG,GG), as ex- 
plicated below, is labeling the different possible interplays 
between the global symmetry and the topological order 
in 3+ld. 

Let us use SET^{SG,GG) to denote the symme- 
try enriched indices SET{SG,GG) in Eq.(l32]), empha- 
sizing that it is for 2+1 dimensions. SET'^{SG,GG) 
contains a rather different mathematical structure than 
SET^{SG,GG): 

SET^iSG, GG) = [H^{SG, Z) (g) H^{GG, Z)] x 

[H'^iSG, Z) (g) H'^iGG, Z)] X Tor[7?'*(S'G, Z), H^(GG, Z)] x 

Toi[H^{SG,Z),H'^{GG,Z)] x Toi[H^{SG, Z), H^{GG, Z)] 

(126) 

One can anticipate the possible non-trivial interplay 
between the global symmetry and the topological order in 
3+ld. Firstly, let's consider the generalization to 3+ld of 
EXTRA{SG,GG), a part of SET^{SG,GG). In 2+ld, 
EXTRA{SG, GG) is describing the interchange of quasi- 
particle species by SG action. We expect that in 3+ld, 
part of SET'^{SG, GG) also describes the interchange 
of excitation species. Note that in 3+ld, the topolog- 
ical excitations can be either point-like gauge charges 
or loop-like gauge fiuxes. More precisely, this part of 
SET'^{SG, GG) should describe the species-interchange 



of gauge fluxes and of gauge charges by the global symme- 
try (but not interchange of gauge flux and gauge charge) . 

Next, we consider the generalization to 3+ld of the 
other part of SET^{SG, GG): SFGiSG, GG). One may 
wonder, does part of SET*{SG, GG) label the symme- 
try fractionalization phenomena in 3+ld? Naively, if we 
consider the symmetry fractionalization classes of point- 
like gauge charges, we should have a similar mathemat- 
ical structure as in 2+ld (see Eg. ([33)) ). In particular, 
when GG is abelian, we should have H^{SG,GG) = 
[H'^iSG, Z)®H^{GG, Z)] X Tot{H^{SG, Z),H^{GG, Z)) 
indices labeling the different projective representations of 
the symmetry group. However, this mathematical struc- 
ture is missing in SET'^'^SG, GG) above. 

In fact, we know that symmetry fractionalization 
classes of gauge charges are actually missing from 
SET^{SG,GG) in 2+ld, as discussed in Sec. EH (see 
also HrC]) . It is not surprising that SET^{SG, GG) is also 
missing those indices. However, we know that the sym- 
metry fractionalization classes of gauge fluxes are com- 
pletely contained in SET^{SG,GG). We expect that 
SET'^{SG, GG) also contains the generalized "symme- 
try fractionalization classes" of gauge fluxes — the topo- 
logical loop excitations in 3+ld. Therefore, a part of 
SET'^{SG, GG) should describe the non-trivial action of 
global symmetry on gauge flux loops, without changing 
their species. But because gauge flux loops are extended 
objects, the action of global symmetry on them can no 
longer be implemented by local operators. So, to be pre- 
cise, we should not call this phenomenon symmetry frac- 
tionalization, as it is discussed in Sec lII Bl We will call it 
the extended symmetry fractionalization. 

Although the full understanding of the extended sym- 
metry fractionalization is beyond the scope of this paper, 
we can intuitively guess the underlying mathematical 
structure. When GG is abelian, a direct generalization 
of H^{SG, GG) to one higher dimension is H^{SG, GG). 
By the universal coefficients theorem, we have: 

H^{SG, GG) = [H^{SG, Z) « H^{GG, Z)] 

xTor[H'^{SG,Z),H^{GG,Z)], (127) 

where we used H^{GG, Z) = GG for finite abelian group 
GG. These two terms indeed appear in Ea. (|126p (the 1st 
and the 3rd term). We propose that H^{SG,GG) is at 
least part of the mathematical structure describing the 
extended symmetry fractionalization classes when GG is 
abelian. 



B. Generalization to continuous groups, and/or 
anti-unitary symmetry groups 

Our discussion has been limited to the case in which 
both SG and GG are finite groups, and SG is assumed to 
be unitary (i.e., does not contain time- reversal). These 
constraints are introduced here for simplicity rather than 
due to difficulty of principle. 
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First, it is quite straightforward to consider an on-site 
symmetry group SG containing the anti-unitary time- 
reversal transformation T. In the work by Xie et al. 
on SPT phasesiS, when SG contains T, the classifica- 
tion is given by H'^+^{SG,Ut{1)). Here C/t(1) means 
that T acts non-trivially on the U{1) group; in particu- 
lar, T sends the phase e*^ G U{\) to its complex conju- 
gate e~*^. (For the detailed definition and discussion of 
H'^+\SG, C/t(1)), see RefUl). We expect that the clas- 
sification of gapped bosonic quantum phases with topo- 
logical order described by a GG gauge group, and in the 
presence of an on-site global symmetry group SG con- 
taining r, is given by H'^+^{SG x GG, Ut{1)), in which 
only the SG part of the cross product acts nontrivially 
on C/t(1). 

Further, Reffig! considered the classification of SPT 
phases with a continuous on-site symmetry group 5G, in 
which case the Borel group cohomology was used. It ap- 
pears to us that the inclusion of a continuous SG results 
in an appropriate (but may be subtle) mathematical gen- 
eralization of i?'^+^(G, U{1)) for a finite group G. We in 
principle do not expect that this generalization is hin- 
dered by difficulties. 



C. About the completeness of the classification 

Finally we comment on the issue of completeness of 
the classification by H'^+^{SG x GG, C/(I)). Is the clas- 
sification complete, incomplete, or overcomplete? 

First, we want to comment on the following question: 
Do distinct elements in H'^^^{SG x GG,U{1)) necessar- 
ily correspond to distinct quantum phases? We believe 
that the answer is negative, and the classification is gen- 
erally overcomplete in this sense. One can understand 
this claim by considering a simple example SG = Z\ 
and GG = Z2 x Z2 in 2+1 dimensions. According to 
H^{GG,U{1)) = Z|, it appears that there are 8 differ- 
ent topological orders. But this is overcomplete. Among 
these three Z2 indices, the first(second) Z2 is labeling 
the toric-code/double-semion topological order in the 
first(second) Z2 gauge group, and the third Z2 is label- 
ing a certain extra twist of the topological order involving 
both Z2 gauge groups. Let's consider two phases labeled 
by (1,0,0) and (0,1,0), where we group the three Z2 
indices into a vector. Clearly (1,0,0) and (0, 1,0) physi- 
cally correspond to the same phase — they just differ by 
an ordering of the gauge group. This is analogous to the 
K-matrix classification of abelian quantum Hall states. 
For instance K = ( J 2 ) ^^-^ K = (01) a-i'e labeling the 
same physical phase. 

The SET{SG, GG) classification of SET phases is also 
overcomplete in this sense. For instance, the symme- 
try fractionalization classes when SG = Z2 and GG = 
Z2 X Z2 are given by SFG{SG,GG) = H'^{SG,GG) = 
H^{SG, Z2) X H'^iSG, Z2) = Z2X Z2, using the universal 
coefficients theorem. If we use (a, 6), with a, 6 = 0, 1 £ 
^2, to represent this index, then (1,0) [(0,1)] simply 



means that there is non-trivial symmetry fractionaliza- 
tion in the first [second] Z2 gauge sector, and these are 
physically the same situation. 

Second, our classification is certainly not a full classi- 
fication of all possible gapped bosonic quantum phases 
with both global symmetry and topological order. There 
are certainly topological orders that cannot be described 
by discrete gauge theories, for instance, the chiral frac- 
tional quantum Hall states. Even for non-chiral topo- 
logical orders, there are phases realized by the string-net 
models*^*^ in which quasiparticle quantum dimensions are 
not integers, which again cannot be described by discrete 
gauge theories, where quasiparticle quantum dimensions 
must be integers. 

So let us ask, under the constraint that the topolog- 
ical order is indeed described by a discrete gauge the- 
ory, is the classification complete? However, we must 
firstly specify the condition: what exactly do we mean 
by "topological order described by a discrete gauge the- 
ory"? We actually mean nothing but phases character- 
ized by H'^+^{GG,U{1)). Unfortunately this sounds like 
a circular argument, but we do not know how to trans- 
late it into a more physical statement. For example, 
H^{Z2,U{1)) = Z\, meaning that there is only one Z2 
gauge theory in 3-|-ld in our context. However the string- 
net models in 3-|-ld allow one to construct two topological 
phases, both of which look like Z2 gauge theories. The 
difference is that in one phase the Z2 gauge charge is a 
boson, while in the other phase it is a fermion. In our 
classification, the second phase is not included. 

Even under this condition, namely the topological or- 
der given by H'^+^iGG, t/(l)), it seems that the classifi- 
cation may still be incomplete. For instance, we already 
mentioned that the gauge charges always have trivial 
symmetry fractionalization in our classification. Maybe 
it is better to ask: Is the classification complete under ad- 
ditional physical conditions? And if yes, what are these 
additional physical conditions? These are difficult ques- 
tions and we currently do not know the answers. Never- 
theless, we can make some comments. Below we address 
two aspects of these issues. 

Is the classification of symmetry fractionalization 
classes complete? We want to further comment on 
the missing symmetry fractionalization classes for gauge 
charges. Let us focus on 2-|-ld, where the gauge charges 
and gauge fiuxes are dual to each other, at least for 
abelian GG. In this case, we can re-interpret the sym- 
metry fractionalization classes for gauge fiuxes as those 
for gauge charges, after performing the duality. So the 
real missing part should be those phases with non-trivial 
symmetry fractionalization for both gauge charges and 
gauge fluxes. 

Naively one may think that it is possible to construct 
such a phase by coupling two phases together. For in- 
stance, consider GG = Z2- Our models can be used to 
construct two phases: phase-a (phase-b) in which gauge 
fiuxes rua (charges Cf,) have non-trivial symmetry frac- 
tionalization (i.e. transform under SG as a non-trivial 
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projective representation), while Ca (jUb) transform triv- 
ially under SG. (Phase-b can be constructed by per- 
forming e — >■ m duality from phase-a.) We can couple 
these two phases together. When the coupling is weak 
the topological order is Z2 x Z2. After a phase tran- 
sition of condensing marrib bound states (or the eaCb), 
the topological order will reduce from a Z2 x Z2 gauge 
theory to a Z2 gauge theory. But the condensed uiamb 
(or CaGb) quasiparticles actually transform non-trivially 
under SGI And consequently the new Z2 topologically 
ordered phase breaks SG. 

However, it seems possible to construct an effective 
lattice gauge theory for the phases with non-trivial sym- 
metry fractionalization (for on-site SG) for both gauge 
charges and gauge fluxes. So it is feasible to expect that 
these phases do exist. But the above discussion signals 
that these phases may not be adjacent to the phases clas- 
sified by H'^'^^{SGxGG, U{1)) in some sense. The reason 
that we miss those phases in H'^~^^{SG x GG, U{1)) may 
be because they cannot be described by exactly solvable 
models. 

Is the classification of the interplay between SG and 
GG beyond symmetry fractionalization complete ? Let us 
consider the classes indexed by non-trivial elements in 
EXTRA{SG, GG) in 2+ld. We show that in the ex- 
ample of GG — Z2 X Z2 and SG = Z2, the minimal 
model for a non-trivial EXTRA{SG, GG), such a phase 
means that SG can interchange the species of quasiparti- 
cles. However, we know that in 2-|-ld, even for the usual 
toric code topological order with GG = Z2, it is fine to 
imagine that SG — Z2 could interchange e and m quasi- 
particles, leaving the fusion and braiding algebra invari- 
ant. In fact, the translational symmetry (not an on-site 
symmetry) along the 45 degree axis of a square-lattice 
toric code model indeed interchanges e and m. Such an 
e and m interchange induced by SG is also missing from 
H^{GG,U{1)). It is however actually reasonable that 
such phenomena are missing in our classification. That 
is because our classification can be generalized to arbi- 
trary higher dimensions, while the e-m interchange can 
only occur in 2+ld. For instance, in 3-|-ld e is point-like 
and m is loop-like, so they can never be interchanged. 

As we prepared this manuscript, we noticed the recent 
work by Ling-Yan Hung and Xiao-Gang WenS, which 
discusses the general duality between SPT phases and 
the Dijkgraaf-Witten TQFTs. Another recent work, by 
Andrew M. Essin and Michael Hermele^, classifies the 
general symmetry fractionalization of gapped Z2 quan- 
tum spin liquids. Also, Yuting Hu, Yidun Wan, and 
Yong-Shi Wu^^. have very recently analyzed in detail the 
topological phases described by models similar to ours in 
absence of global symmetries. 

YR thanks helpful discussions with Fa Wang and es- 
pecially Michael Hermele. This work is supported by the 
Alfred P. Sloan foundation and National Science Foun- 
dation under Grant No. DMR-1151440. 



Appendix A: Projective symmetry group in parton 
construction 

Parton construction is a convenient way to obtain 
quantum states with topological order. The basic idea is 
to write down a topologically ordered state directly us- 
ing the anyonic quasiparticle degrees of freedom. For ex- 
ample, in the Schwinger-fermion representation of quan- 
tum spin liquids with Z2 topological order, fermionic 
quasiparticles are used to represent a physical spin- 1/2: 

Si — l/2j}^cTaj3fip, where i labels sites and a,/? label 
spins. Note that this construction enlarges the Hilbert 
space from 2 to 4 per site, and one eventually needs to re- 
move the unphysical states (empty and doubly-occupied 
sites) to obtain a physical spin-1/2 wavefunction. This 
removal of unphysical states can be accomplished by the 
so-called Gutzwiller projection: Pq = ni"i(2 ~ "■i)j 
where = fj^fia is the fermion number on site-i. 

In this approach, on the mean-field level, a Z2 quantum 
spin liquid can be represented as a free fermion state 
\^mf) of fia fermions, which is the ground state of a 
spin-singlet mean-field Hamiltonian: 

HmF = X^3fLfJa + A,,e„^/^/j^ + h.c. (Al) 

Hmf has both hopping and pairing terms on a lattice 
{ea/3 term is the spin singlet pairing). The physical spin- 
1/2 wavefunction of the Z2 QSL can be obtained by 
Gutzwiller projection: \iPqsl) = Pg\'4^mf)- Here the Z2 
gauge fluctuations emerge exactly because of this projec- 
tion: two mean-field states differing by a gauge trans- 
formation Xij ^iXij^j^^ij ^i^ij'^j i^i = ±1) give 
exactly the same spin wavefunction. Therefore such local 
Z2 fluctuations correspond to redundancies in the formu- 
lation and are gauge fluctuations. The fia fermions are 
the quasiparticles carrying the Z2 gauge charge. The low 
energy effective theory of the state IV'qsl) is described 
by Z2 gauge charges fia coupled with a dynamical Z2 
gauge field. 

How can we make sure that the QSL wavefunction 
iV'QSi) is symmetric under a symmetry group SG, such 
as lattice translations? Naively one would require the 
mean-field Hamiltonian Hmf to be invariant under SG 
transformations. In fact, this is not required. Because 
two mean-field states differing by a Z2 gauge transfor- 
mation label exactly the same physical state, one only 
requires Hmf to be "projectively symmetric". Namely, 
Hm f before and after an SG transformation can differ by 
a Z2 gauge transformation. This is the key observation 
underlying the PSG. 

For any element g in SG, there will be a certain Z2 
gauge transformation Gg associated with the g such that 
the combination Gg ■ g leaves the mean-field Hamiltonian 
Hmf or the state IV'mf) invariant. The collection of all 
such combinations form a group, which is defined to be 
PSG: PSG = {Gg ■ g : Gg ■ g leaves Hmf invariant, V5 G 
SG}. 
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Let's look at the mean- field Hamiltonian Eq. (jAl[) 
again. Clearly, one can do a global gauge transforma- 
tion: Si = — l,Vi and Hmf is invariant, which is also the 
only non-trivial gauge transformation which leaves Hm f 
invariant. This means that in PSG, there will be two 
elements corresponding to the identity element in SG: ei- 
ther a = 1, Vi, or ei = — 1, V«. In general, for any element 
g G SG there will be two elements, Gg ■ g and Gg ■ g in 
PSG corresponding to it. And the gauge transformations 
Gg and Gg differ by the global Z2 gauge transformation. 
Mathematically, the algebraic relation between PSG and 
SG is given by^: 

PSG/IGG = SG. (A2) 

Here IGG~Z2, the group of global gauge transforma- 
tions. 

Eq. (|A2[) is the key mathematical structure underlying 
PSG. It indicates that PSG is a group extension of the 
group SG by IGG. When IGG is abelian, which is true in 
the Z2 case, we know that IGG is in the center of PSG, 
because global gauge transformations obviously commute 
with any PSG element. In this PSG is a central 

extension of SG by IGG. Further, the classification of all 
different PSGs becomes the classification of all possible 
central extensions. There is a nice mathematical theo- 
rem on central extensions of groups stating that all such 
central extensions are classified by H^{SG, IGG){see, for 
example, Ref^S'). 

At this moment, it appears that PSG is a feature 
of the parton mean-field states only. Whether PSG is 
physical or not beyond the mean-field formulation is not 
completely clear. To see the physical meaning of PSG 
and H^{SG, IGG) beyond the mean- field formulation, 
we need to consider the low energy effective theory, which 
is discussed in the main text. 



Appendix B: The operator realization of twisted 
extended ribbon algebra and the quasi-quantum 
double 

In this section we write down the explicit forms of the 
ribbon operators i^('''9)(r) and the operators acting on 
ends of ribbons D(^f^ g^{A),D(^f^^g^{B), and demonstrate 
the algebra satisfied by these operators. As in the main 
text, we also only consider an Abelian g roup SG GG, 
These operators are defined for h 6 GG, g G SG x GG. 
In addition, we study the braiding and fusion properties 
of the quasiparticles created by these ribbon operators, 
and show that they are mathematically described by the 
quasi-quantum doublcj^i^ 

1. The operator realization of the twisted extended 
ribbon algebra 

Let us start by noting that the definition and 
some properties of ribbon operators were presented in 



Sec. lIVB] Here we start by recalling the form of the non- 
zero matrix element of the ribbon operator F^'^'^'^ (T): 

{&n\F^^^3)^r)\i)^fA-fB-fAB-wl{g), (Bl) 

where wf^{g) is defined in Eq. (j72p . and the relation be- 
tween the initial and final states, |i) , |fin), is explained in 
Eq. ([7T|) : the /aj/bj/ab are rather complicated phase 
factors depending only on the degrees of freedom living 
on ends of F, which we here explicitly define: 

^ _ uj{hAh~'^ ,h,bM)uj{hN,h,hAh~'^)ijj{e'^\^,hNC^^^,h) 

(jj{cN+i,bNC^\i,h)u(hAh-'^ ,h, h^'^)uj{h^'^ ,h, hAh-^) ' 
j,^ _ ChBih,hB)chhB{ci,h) 
ChhB{hB,h)cc-,{h, Hb) 
U}{h, h~^c^^bo, e^^)w(ci, h, h~^c^^bo) 
uj{hB,h,hg^h^'^) 
fAB^i^-\h,h^\hAhB). (B2) 

Here the fiux in Ia is Ha — bJj^CM+ieN+i and the flux in 
t_B is liB = 6q ^ciei in the initial state {gij})- (see 

FigCSfor the definitions of bo,bN--- degrees of freedom 
living of the ends of F) . Although here we use the specific 
geometric configuration of the ribbon F in Fig |15[ the 
definitions of i^('*'f)(F), Siud D(^h,g){ A) {D,^h,g){B)) below, 
can be easily generalized to any geometric configuration 
of F. 

Next, we define the operator ZJj/j g)(v4) (_D(^ ex- 
plicitly on end-A (end-B) (also mentioned in the main 
text), where h G GG, g € SG x GG. They are defined as 
an operator in the whole Hilbert space H (In fact they are 
also well-defined in/C(F), because D(h^g){A),D(h^g)[B) do 
not change the flux of a 2-simplex inside F), via the ma- 
trix elements: 

(f| D^y,,g){A) |i) = 5h^,h ■ CH{g,bN)W(,{i) 
{l\D(,,^g){B)\i)^5hB,h-CH{g,ci)We{i), (B3) 

where the phase factor We{i) is defined in the main text 
in Eq.®. 

With these definitions, after straightforward but com- 
plicated algebra, one can show that, in both the Hilbert 
space % and its subspace /C(F): 

D(h2,g2){A) ■ Di^h^,g^)(A) = Sh,M ■ ChA92,9i)D(h,.,g2g,){A) 

-^(/i2,S2)(-^) • ^{hi,gi){B) = ^h^M ■ ChA92, gi)D(ht,g2gi){B)- 

(B4) 

In the sub-Hilbert space /C(F), we have more identities: 

-c,, {hMch, (32, 9i)D(h..g,) (A)^^("-3^) (F), (B5) 

i?(,,,,,,,)(i?)F(''-^'^^^)(F) 
=c,, {h,,h2)ch, igu92)F^^'''^\T)D^^^^g^^{B), (B6) 



42 



and, 



(B7) 



Eas. (|B4IB5IB6IB7t are summarized in 

Eas. (|81al81cl81dl81bp in the main text. 

In addition, it can be easily shown that D^iig-j{A), 
D(h.g){B) and F^^'^^iV) operators all commute with the 
global symmetry transformations in SG. 



2. Braiding 

a. Braiding between the quasiparticles at end-A 's 

Next, we describe the braiding and fusion proper- 
ties of the quasiparticles created by the F^'^iT) opera- 
tors when applied on the ground state \gs). For these 
purposes, we must consider multiple ribbons. Start- 
ing from multiple ribbons Fi, F2, ...F^r, we define a sub- 
Hilbert space /C(Fi, F2, ...F^v) C H to be the one spanned 
by those states satisfying zero-flux rule everywhere in- 
side Fi,F2,...FAr. Now ah the operators F('''9)(r„), 
D{h,g){An)-, D(^h,g){Bn), = l,2,...iV are ah well-defined 
in /C(Fi, F2, ...F^v)- {-An and Bn are the two ends of the 
ribbon-F„). 

In order to make the braiding and fusion algebra of the 
quasiparticles at both the An and Bn ends on the same 
footing, it is convenient to introduce a slightly modified 
ribbon operator: 



(B8) 



where the phase factor Jab is given in 
also commutes with the global symmetry, and satisfies 
the same operator algebra as F^'^^^^T) with D{A), D{B) 
operators (see Ea. (|B5IB6|) ). (But the algebra in Ea. (jB7|) 
IS no longer satisified by F'-'''9\T).) We can then con- 
struct an excited state using multiple ribbons Fi,...Fjv 
that all share the same end-B, but end-A„ are all dif- 
ferent from each other. Here by "sharing the same end- 
s'' , we require that the vertices i b„ , the triangles t b„ , 
and the edge-i3„'s are all the same: Bi = B2 = ■■■ = 
Bn = B, edge-Bi — edge-i?2 — ■■■ — edge-B^ = edgc-B, 
Let's only consider the case that Fi,...Fjv do not over- 
lap/intersect with each other except over a finite length 
starting from end-i?, which is enough for our purposes. 
We can choose an ordering of the ribbons: F„ is on the 
counter-clockwise side of F„_|_i, Vn. (See FiglTTk for a 
geometric illustration.) An excited state is given by : 

(B9) 

where we have used fci = (hi, fci)...fcjv — (/ijv, ffjv) to save 
notation. (B) is a projection operator that enforces 
= ub & SG at end-B, which commutes with all the 



F-operators. IV'us"'^") only contains quasiparticles at 
end-An and end-B, and will be very useful to understand 
the braiding and fusion algebra. We denote the sub- 
Hilbert space spanned by all the states IV'mb"''") with a 
fixed ub as Cub (^i: ^2, ■■■An, B), because one can show 
that this sub-Hilbert space only depends on the ribbons' 
ends, but does not depend on the paths of the ribbons. 

We will soon study the braiding and fusion op- 
erations of the quasiparticles at the end-A„'s in 
£„^(Ai, A2, ...An,B). These operations only act within 
the sub-Hilbert space £ub{^i,^2, ■■■An,B) for a fixed 
ub, because mb is always unchanged in them. In fact, 
we will show that the braiding and fusion algebra (or, 
the superselection sectors of end-An quasiparticles) only 
involve D(^,i^ fi^-){B), with /ii,/i2 G GG. 

However, when a general D(^g)(_B) operator 
(g = hg ■ g, hg £ GG and g e SG) acts on a state in 
£„g(Ai, A2, ...An,B), it will send the state to a different 
sub-Hilbert space Cg. u^iAi, A2, ■■■An , B). In addition, 
the global symmetry transformation of g also sends 



Cub{Ai,A2,...An,B) to £5 



XAi,A2,.-An,B). 



Therefore if one wants to study the general D^^f^g-^^B) 
operators and the global symmetry transforma- 
tions, one should consider a larger sub-Hilbert space: 
C{Ai,A2,-An,B) = ®UBeSGCuB[Ai,A2,-AN,B). 
Although we do not have a proof, we believe that 
£(^1,^2, ...AntB) contains all possible excited states at 
end-A„'s and the end-i3. (For single ribbon states IV'mb' 
we do have a proof in the main text that they span all 
possible excited states at the end-A and the end-B.) 

Note that because /ia„ = 1 in \gs), Jab — 1 due to the 
canonical form of 3-cocycle. Therefore [t/j'^i - '^n^ ^.g^jj |-,g 
equally created by the original i''('''f)(F) ribbon operator: 

lA'^-'") = PuB{B)F'''{Ti)F''nT2)-F'''^{TN)\gs), 

(BIO) 

which is given in Ea. (j89p in the main text. In fact, to 
study the braiding and fusion operations of the end-A„ 
quasiparticles, it does not matter whether _F('*'ff)(r) or 
p{h,g)^^^ is used, because they give the same algebra. 
This is why we only use _f(''^9)(F) operators in the main 
text for simplicity. 

The geometric illustration of the braiding process be- 
tween an end-An and an end-A„_|-i has been discussed in 
the main text, but is also shown in Fig. [T7b. The follow- 
ing operator identity is crucial to compute the braiding 
algebra: when Fi and F2 share the same end-B but have 
different end-A's, and Fi is on the counter-clockwise side 
of F2: 

p(hi,gi)(j-^-^p(h2,g2) (-p^-j 

=c..(g2fer\fei)^"'-^-''-'''(ri)F(''^-^-Hr2) "i^'^!;^i . 

w(/l2,/ll,/ls) 

(Bll) 

Here hB should be understood as an operator that mea- 
sures the gauge flux in the 2-simplex tB^ The order be- 
tween the the _F-term and the w-term on the right-hand 
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a) 





FIG. 17: Multiple ribbons sharing end-B, describing quasiparticles at their end-A's. (a) Example of ordering multiple (five) 
ribbons, with strings Fi . . . Fs, in a "counter-clockwise sense". The ribbons are allowed to have overlaps only over a finite length 
starting fi"om end-B. The figure shows a realization possible on a triangular lattice, (b) The counter-clockwise 180° braiding 
of particles A„ and j4„+i in an A'^-particle state, a generalization of Fig. 1161 The (blue) strings F' apply to the braided state. 



side is therefore important. Eq. (|Blip already describes 
the physical counter-clockwise braiding (180°) operations 



^cc^^ between quasiparticles at end-yl„ and end-A„+i 
in £.{Ai, An, B) completely: 



CC \VUB 



■ihn,gn).{h„ + i,gn + i). 



=p„,(i?)F'=nri)...F(''"'^"Hr„+i)i?(""+^'^"+^Hr„)...#'=~(r^)|gs) 

W(Al„-|_i, tin, 'ln+2'in+3---"-iV) 

w(/l„, /l„+l, /l„+2^n+3---'lAf) / , -1 , M/fei. 

= -77 T T 7 T-^Ch„+A9n+lhn ,/ln)|l/'«B 

W(ft„+1, tin, tln+2nn+3---'lN ) 



.(h„ + i,3„ + i/i„ ^),(/i„,g„)...,fejv 



(B12) 



Here the Hb in Eq. ljBlip picks up the accumulated flux 

in tB- hn+2hn+3---hN ■ 

Although everything about braiding can be understood 
from Ea. (IB12[ ). the w-factor in it is not convenient. Can 
we get rid of this phase factor and find the underlying 
algebraic structure satisfied by Uq^^^I 

When acting on the ground state, /is = 1, and the 
w-tcrm in Eq. (jBll[ ) is 1 due to the canonical form of 
a 3-cocycle. In this particular case Eq. ljBllI) becomes 
Eg . (j99| . The braiding of an excited state with only two 
ribbons ri,r2 is indeed described by Ea. (P7|) . Let's de- 
fine a formal braiding operator i?^^ which implements 
the physical braiding R\f(j in 

R]fc ^R'^'-Di^^B) ® D(^\B) ■ a 
=a-R'-iDi.^\B)(g,Df\B) 



hi,h2&GG 



the tensor Rihugi),{h2,g2) ^ Sh.i,g2Sg^ t is also defined in 
the main text, and ^ ^{B) with hn,Pn G GG is de- 
fined to be a formal operator that only transforms the r„ 
operator F('''9)(r„) as if we are in a single ribbon state 
(see Ea. (|B6|) ). More precisely, because 

(B14) 

we have 



ihi,giPi ),(/i2,g2P2 ) 



(B15) 



where a is the permutation operator: o'IV'ub ) = IV'ub )> 



The key property of these formal tensor product oper- 
ators is that they commute with any local operator at 
end-^i, Vi. That is why they are topological operators 
for the quasiparticles at end-A^'s. 
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Because of Ea. (IB4[) . the multiplication of the formal 
operators D^^^^^{B) {yh,p £ GG) satisfy the following 
algebra: 

D^-\b)d'^;\b) = ^1%d'i:\b), (B16) 

where we used the tensor iTjl = 
hi,h, Sh^M ^gk,gigj Chk{9i,9j) defined in Eq.®. 

Ea. (jB16l) tells us that the quasiparticles created by 
ribbon r„ at end-A„ form representation of this al- 
gebra. Mathematically, the algebra in Eq. (jB16P is 
called the multiplication in the quasi-quantum double 
D"(G'G)55'6l, where Q e H^{GG,U{1)) is the 3-cocyle 
induced on GG by the cocycle w e H^{SG x GG, U{1)) 
in our model by restricting the elements x,y,z G GG 
in uj{x,y,z). Multiplications in D'^(GG) is associative. 
A representation of the multiplication algebra Eq. (jB16[ ) 
is called a representation of the quasi-quantum double 
D'^(GG). 

Because we will show that the braiding algebra is com- 
pletely determined by D^^'^^^{B) operators with h,p E 
GG, one knows the braiding properties of a quasiparticle 
at end-A„ if we know which representation of D"(GG) 
that this particle is in. In fact, the quasiparticle species 
(or more precisely, its superselection sector) is labeled 
by an irreducible representation of D"(GG). And differ- 
ent irreducible representations of D'^(GG) correspond to 
different quasiparticle species. 

The physical counter-clockwise 360° braiding for the 
two particle states is {Rqq~^^)^. Its action in the basis 
IV'u^ is actually very simple, because the oj-terms in 
Eq.dHTH cancel out for {R'^'^'^'^)^: 

(B17) 

It is also tempting to formally define the general oper- 
ator 

=a ■ R'''^dIJ'\B) ® D^^+^\B) 

h,^,h„ + i£GG 

(B18) 

Note that the formal operator R^;^^^ does not has a hat, 
which distinguish it from the physical braiding Fi^'^^^ in 
Ea. (|BT2l) . The w-factor in Eq. ((BT2l) tells us that, for 
multiple ribbon states, the physical braiding R^'^^^ is 

not implemented by in the basis {IV'uL' 

But we will show that if we change into certain different 

basis, R^'q~^^ is still implemented by R^'^^^l 

Let's firstly consider three-ribbon states: \ipu]^^'''')- 
The physical Rqq between the end-A2 and the end-Aa 



are still implemented by R^q. But in this basis, Rqc 
and the physical braiding Rqq differ by a phase factor 
Zihl'.Mhf) ' w-term in Eq. (jBll[ ). However, this 

braiding can still be implemented by Rqq if we choose a 
different basis. 

Let's define the state: 

l^i^'^^'^^^'^^) ^^-'(/ii,/i2,/^3)|V'^L'''^) (B19) 

For reasons that will become clear in a moment, we also 
define: 

All states {IV'ls'''''^''''-')} form another basis of 
£„g (Ai, ...Ajv, S), with differ from the original ba- 
sis { I "013^ '^'^^''^^ ■*■')} only by phase factors. 

It is then clear now that the physical braiding Rqq is 
implemented by R^c basis {IV'S ''^^ '^^'''')} (but not in 
the basis {IV'mb^ '^^^''^^ and the physical braiding R^^ 
is implemented by R^j^ in the basis (but 

not in the basis {IV'Ib^''"'^^''^^'''')})- Such a basis change 
is necessary to maintain the same algebraic form of the 
formal braiding operator i?^'^^^ in Eq. ljBISp . 

The two different parentheses configurations: 
(^1,(^2,^3)) and ((^1,^2)1^3) can be viewed as two 
different ways to "multiply" quasiparticles. For reasons 
that will become clear later, it is better to call these 
operations as "co-multiplications". And different orders 
of "co-multiplications" do not give the same results, 
which differ by a basis change: mathematically, the 
"co-multiplications" are not associative, but are quasi- 
associative. It turns out that the "co-multiplications" of 
quasiparticles here have a clear physical meaning: the 
fusions, which we will discuss shortly. 

One can generalize the above observation for 3-ribbon 
states to multiple-ribbon states, which turns out to be 
also very useful to represent the fusion algebra. We de- 
fine: 

|^(fcl,(fc2,...(fejV-2,(fcjV-l,fcjv)))---)^ — |^fclfc2...fcjV~2fe]V-lfcjV^ 

(B21) 

We now consider an arbitrary parentheses configuration 
between fci, A:2, ^at. For convenience, we use a tree 
diagram to represent it. For example, (^i, (fe, fcs)) and 
((/ci, fe), fcs) can be represented as figure- (a) and figure- 
(b) below, while {ki, {k2, ...{kN-2, {kN-i, kN)))---) is rep- 
resented as figure- (c). 




/ii/j2/i3 h\h2hz h\h2.--h]si 



(a) (b) (c) 
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Note that we define a tree diagram not only as a ge- 
ometric object: First, it contains fci = (hi, gi), ...kj^ — 
(/iATjffAf) assigned for the top end points. Second, every 
edge (Hne segment) in a tree diagram is also assigned a 
group element € GG, which is specified as follows. The 
top edges are assigned as /ii, /i2, .../iat, and every lower 
edge coming out of merging two upper edges is assigned 
by the product of the group elements in the upper two 
edges. 

We can use a tree diagram to represent any parentheses 
configuration: 



ki k2 ....kN-2 kN-i fcjv 



A general tree diagram: 




We now define a basis for a fixed tree diagram Treea, 
where a labels the tree configuration: 

|V„«(Tree„)) = u;(Tree„)|V;^;^'=-'=«-^'="-^*'-"), (B22) 

where w(Treea) is a phase factor which we define below. 

Any two tree diagrams with the same set of assigned 
ki = {hi,gi), ...Un = {hN,gN) on the top end points can 
be deformed into each other by a finite number of so- 
called F-moves. An F-move is a local deformation of a 
tree diagram, namely Treea and Tree/3 are only different 
locally as shown below: 

Treea Tree^ 



chooses. This is a consequence of the 3-cocycle condition. 
3-cocycle condition dictates that the F-moves satisfy a 
crucial self-consistent condition: the pentagon equation, 
which in turn indicates that li'(TreeQ) is well-defined — a 
consequence of the Mac Lane's coherence theorem^ We 
refer interested readers to Refl67l by Kitaev for detailed 
discussions. 

With this definition of [i/'mb (Treea)), one can show 
that the physical braiding operation R^^'^^ between end- 
An and end-^„+i is implemented as Rqq^^ in any ba- 
sis in which /c„ and fc„_|_i are parenthesized together: 
...(fc„, fc„ -|- 1).... 

Mathematically, the formal operators Rqq~^^ does not 
satisfy the Yang-Baxter equation, a self-consistent equa- 
tion for braiding algebra. This dictates that an appro- 
priate changing of basis is required, which is discussed in 
detail above. With this changing of basis, the formal op- 
erators satisfy the so-called quasi- Yang-Baxter equation. 
These mathematical structures are exactly those in the 
quasi-quantum double D"(GG). 



b. Braiding between the quasiparticles at end-B 's 

Similarly we can construct an excited state using mul- 
tiple ribbons Fi, ...r„ that all share the sa end-A: Ai = 
A2 = ... ^ An = A, edge-^i = edge-A2 = ... = 
edge- An = edge- A, and that do not overlap/intersect 
with each other except for a finite length starting from 
end-A, while end-i?ri are all different from each other. 
We also choose an ordering: F^ is on counter-clockwise 
side of fi_|_i. (See Fig|T5k for a geometric illustration.) 
We can create an excited state: 




hah),hc hahhhc 



U3 ^{ha,hi,hc)\ipnB(T^eea)) = I^Iub (Tree^)) 

FIG. 18: A general F-move. 

When a local F-move occurs, we define the state 
a;"i(/ia,/ifc,/ic)|V'«B(Treea)) = |V^„b (Tree^j)). In this 
fashion, starting from the tree diagram of the original 
basis (fci, (^2, ...(fcjv_2, (fcjv_i, fcjv)))...), we can find the 
phase factor w(TreeQ,) accumulated during a sequence of 
F-moves for any tree diagram a. 

Apparently there are many different possible F-move 
paths that can connect a given tree diagram with that 
of (fci, (A;2, ...(fcjv_2, (A;Ar_i, /cat)))...). One may wonder 
whether the accumulated phase factor ^(TreeQ.) is the 
same or not for different F-move paths. It turns out 
that w(Treea) is independent of which path that one 



l^t-"-) = P„.(A)F'=^(fi)F'=^(f2)...F'="(f„)|5s), 

(B23) 

which only hosts quasiparticles at end-Si and end-A. 
Here Pu^ is a projector which enforces Ui^ — ua for 
a fixed element ua G SG. (Note that it is important to 
use F('''f)(r) operators, not the ^('''"^(f) operators to 
construct IV'uii"'''")-) We will use IV'ui'"''™) to study the 
braiding properties of the quasiparticles at the end-i3„'s. 
We denote the sub-Hilbert space spanned by {IV'^^'"'^™)} 
as Cua{Bi,B2, ...,Bn,A). 

The geometric illustration of the counter-clockwise 
180° braiding of the end-i?„ and the end-i3„+i is shown 
in Fig |19b . The following operator identity is crucial to 
understand its underlying algebraic structure. When Fi 
and r2 share the same end-A but have different end-B's, 
and Fi is on the counter-clockwise side of F2: 

=C.,(/.2,3l/^2-^)^<''-^^nfl)F('^-^^'^^"^nf 2)^4^4^ 

(B24) 
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FIG. 19: Multiple ribbons sharing end-^, describing quasiparticles at their end-_B's. (a) Example of ordering multiple (five) 
ribbons, with strings f i . . . Fs, in a "counter-clockwise sense". The ribbons are allowed to have overlaps only over a finite length 
starting from end-yl. The figure shows a realization possible on a triangular lattice, (b) The counter-clockwise 180° braiding 
of particles Bn and Bn+i in an A''-particle state. The (blue) strings F' apply to the braided state. 



Here Ha should be interpreted as an operator measur- 
ing the gauge flux in tA- /ia = 1 in \gs) and Ha — 
h^'h^\..h],' in IVi^--'^"). 

Eq. (|B24p already describes the physical counter- 



clockwise braiding (180°) operations R^;^^^ be- 
tween quasiparticles at end--B„ and end--B„+i in 
C{Bi, i?7v, A) completely: 



j^^^'^+^\^^ki,...{hr,,g„),(hr, + i,g„ + i)---,kN^ 



-R 



^cc 

n,n+l 



CC 



P„^(A)F'=i(fi)...F(''-f")(f„)F(''"+i'^'"+i)(f„+i)...F'^"(f^)|5s) 



=P„,(v4)F'=nfi)...F(''-S")(f„+i)i^(""+^'^"+^'(f„)...F'="(fAr)|gs) 

(B25) 



w(/l„+l, hn, hn+2hn+3---h]s) 



Uj{hn,hn+l,h 



Mn) 



/, , _i SI 7/ci,...(/i„+i,g„ + i),(h„,g„/i J )...,fejv. 



Here the h']^ in Ea. (IB24[ ) picks up the inverse of the 
accumulated flux in Ia'- hn+2hn+3---hN. 

We can also define the formal braiding operator: 

=cr-i?"-D(")(yl)«) (v4) 

h^.h,^ + ieGG 

where a is the permutation operator, and we define the 
tensor i?(''i.9i).('»2,S2) as: 

(B27) 



r 



In fact one can show that for an abelian group, 
Chi{h2,h2^) = c;ii(/i^\ /12), Vfti,/i2. dI^^^^^^^{A) with 
hn,Pn G GG is a formal operator that transforms the r„ 
operator F'-'''3\r„) as if we are in a single ribbon state 
(see Eq. ()B5P ). More precisely, for instance, because: 



='5^-Ju-S/('*r\/^i)%'(Pi,3i)lC"'^^^^) (B28) 
we have: 

. |^(^'ii,pisi),('i2,P2S2)^ (B29) 



D 



(n) 



{A) also satify the multiplication algebra 



47 



Ea. (IB16| ) in the quasi-quantum double D'^{GG). 

One can show that for two-ribbon states, the w-factor 
in Eq. (|B25P is unimportant because it equals one due to 
the canonical form of a 3-cocyle. In this case one can 

show that the physical braiding i?^^ is indeed imple- 
mented by the formal operator R]j^. 

For multiple- ribbon states, the cj-factor in Ea. (jB25P 
becomes important. Similar to the braiding of end-A's, 
one way to get rid of the cj-factor in Eq. (|B25p in the 
braiding algebra of end-B's is to introduce appropriate 
basis changes, which we do not dicuss here. After the 

appropriate basis change, physical braiding ff^^^^ can 
still be implemented by the formal operators R^^^^ 

Finally, the physical counter-clockwise 360° braiding 
{R^^^'^^Y c'^lso has a simple algebraic form: 

(B30) 

3. Fusion 

We study the fusion of the quasiparticles at the end- A's 
only. Let consider multiple-ribbon states with ribbons all 
sharing the same end-i?, but with different end-A's. 

Firstly let's consider two-ribbon states IV'ub'^^)- 
cording to the twisted extended ribbon algebra Eq. (jB6l) . 
the p)(i3) operators with both h,p G GG transform 
two-ribbon states as: 

XI [Cg{hl,h2) 
hi-h-2 = h,hiJi2eGG 

=A;i^p)i^f)(B) )(B)|V4i-^^)'(^-^^)), (B31) 

where we have used the formal operators d\'^\b) 
introduced in Eq. (|B13l) . and the tensor A^^ — 
^g,,gj ^9k,gi Shk.hihj Cg^{hi, hj) defined in Eq.llSS]). 

This motivate us to generally define formal operators 
for multi-ribbon states: 

= A;^.Df ^(B) ® (B32) 

Basically, when acting on a multi-ribbon state, 
^(n)(n-i-i)^^^ only transforms the r„ and r„+i ribbons 
as if we are in a two-ribbon states. Mathematically, the 
operation in Ea. (IB32p is called co-multiplication in the 
quasi-quantum double D'^(GG'), and is often denoted in 
mathematical literature by: 

A{Dr)=AlW^(g)Dj. (B33) 




FIG. 20: Braiding quasiparticle Ai with the fused quasiparti- 
cle A2A3 can be calculated using the shown formula for a two- 
step process. The resulting braiding operator is determined 
by the action of two topological operators (Eq. (|B35|l ). analo- 
gously to the simple case of two-particle braiding; however, in 
this case one topological operator relates to particle Ai, but 
the other is a "comultiple" (Eqs. (|B32|) . (|B33|l ') of topological 
operators acting on particles A2 and A3. 

^(ri)(n+i) j-^-j |-^j.j^g very useful, we will show 

soon that braiding properties of the fused quasiparticle 
of r„ and r„+i ribbons are completely determined by 

Because dI^^^^^^^B) has a physical interpretation of 
acting Dr{B) on the two-ribbon states, dI^^^^^^^B) 
clearly also satify the multiplication algebra Eq. (jB16|) in 
the quasi-quantum double D"(G'G). Therefore, if we 
know the irreducible representations of D'^(GG') (i.e., 
superselection sectors) for the quasiparticles at end-A„ 
and end-A„+i, the co-multiphcation Ea. (|B32p induces 
another representaion of D"(GG), which is generally re- 
ducible. One can decompose this induced representation 
into its irreducible components. Every irreducible com- 
ponent corresponds to one fusion channel. This procedure 
defines the fusion rule. 

Now let's consider three-ribbon states \'4'u]^^'^^)j which 
allows us to study the braiding of a fused quasiparti- 
cle with another quasiparticle. Let's imagine that we 
fuse the end-A2 and the end- A3 quasiparticles first, and 
braid the fused particle with the end-j4i quasiparticle. 
We should physically braid both end-A2 and end- A3 with 
end-Ai, which we define as R^^'^^ = R^q • Rcc (^^^ 
Fig. d^ ). To understand the algebraic structure of this 
procedure, the basis change that we introduced earlier 
becomes useful now. One can straighforwardly show 
that: 

(B34) 

=CT • R'-^Di^^B) ® D(2)(3)(^)|^(^fci,(fc2,fc3)))^ (B35) 
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where the permutation 



operator a is defined as 
Therefore the braiding 



algebra of the fused particle satisfy the same algebra as 
in Eq. pTS)) after using Df^^^\B) operator. 

Similarly, one can imagine to fuse the end-Ai and the 
end-A2 quasiparticles first, and braid the fused particle 
with the end- A3 quasiparticle. This physical process is 



n(l,2),3 ^ Al,2 

^cc — ^cc 



2,3 

cc- 



One can also show that: 



••cc 



), (B36) 

where the permutation operator a is defined as 

This discussion can be easily generalized to multiple 
ribbon states. One can show that the braiding alge- 
bra of the fused particle is always represented using the 
in a basis where (fc„,fc„4-i) are parenthe- 
sized together. 

Finally, the fusion algebra is formally represented by 
the CO- multiplication in D"(G'G) in Ea. (IB32p only in a 
basis where {kn,kn+i) are parenthesized together. Co- 
multiplication is not associative but is quasi-associative; 
namely it becomes associative after the changing of ba- 
sis: the F-move introduced earlier. Because F-moves sat- 
isfy the pentagon equation, one can show that the fusion 
algebra also satisfy the pentagon equation^iiS, the self- 
consistent equation for fusion algebra. 

4. Summary 

In this section we find the operator realizations of 
the twisted extended ribbon algebra, and also study the 
braiding and fusion properties of the topological quasi- 
particles created by ribbon operators. We find that 
the topological order in our model is described by the 
quasi-quantum double D'^(GG'), where the cocycle u £ 
H^{GG, U{1)) is the one naturally induced by the cocycle 
uj e H'^{SG X GG, [/(!)) in our model. 

The core mathematical structures of the quasi- 
quantum double D"(G'G) include the multiplications 
in Ea lB16[ the co- multiplications in Eq. (jB32|) and the 
changing of basis described by Fig. ([T^ (mathematically 
called associator). The superselection sector of a quasi- 
particle is determined by the irreducible representation of 
the multiplication algebra Eq. (jB16p . The braiding alge- 
bra of quasiparticles is determined by the formal braiding 
operator in Eq. ljBlSp . together with the changing of basis 
(associator), which satisfy the quasi- Yang Baxter equa- 
tion. The fusion algebra of quasiparticles is determined 
by the co- multiplication algebra Eq. (jB32p . together with 
the changing of basis (associator), which satisfy the pen- 
tagon equation. One can further show that the braiding 
algebra the fusion algebra are compatible: they satify the 
hexagon equation^^. 

We have not studied the interplay between the global 
symmetry SG and the topological order D'^(GG) here. 



For instance, we have not used D(/j_g)(A), Di^h^g-^{B) 
operators when g ^ GG except for stating their ba- 
sic properties in the twisted extended ribbon algebra 
Ea. (|B4IB5IB6IB7p . However, in SecEl we carefully study 
the interplay between the global symmetries and the 
topological orders in some examples. We believe that 
those studies can be generalized to any phase in our clas- 
sification. 



Appendix C: Particle statistics directly from 
crossing strings 

Here we provide an alternative, direct approach to 
braiding statistics, and compare it to that of Sec- 
tion!^ 

The braiding statistics of two quasiparticles is a topo- 
logical property that therefore cannot depend on the de- 
tails of the ribbon operator at its ends. When the system 
is a torus, the statistical phase of quasiparticles (mean- 
ingful for Abelian quasiparticles only) follows from the 
commutation of operators Ta which describe the follow- 
ing process: creating a particle — anti-particle pair, tun- 
neling the particle across the system in direction a = x,y, 
and finally annihilating the pair, as described by the for- 
mula^ 



i29 



(CI) 



with 6 the (exchange) statistical angle of the quasiparti- 
cles, see Fig. [2TJa). The ribbon operator i^('''9)(r) rep- 
resents exactly the operation T along its F. Since the 
string ends are not involved, we can actually also use 
7«r(g) (seeFig.IH]). 

The angle in Eq. (jCll) depends only on the commuta- 
tion relation at the intersection point of the two strings, 
so even if the system is not a torus, we can consider a 
braiding operation performed locally^! and realized by 
an open and closed string that intersect at a single point, 
Fig.inia). 

Fig.HHshows the ingredients for calculating the phases 
Wmn = (f| W/i" {gm)Wf^" ((;„) |i), which reveal the statisti- 
cal angle: 



i2e 



W12 
W21 



(C2) 



The phases Wmn differ due to the 3-simplices positioned 
on top of two lattice triangles which are shared by the 
two ribbon operators, as well as due to string phases Op" 
on two edges. Fig. W\\h) . 

Using the definition of the 2-cocycle c^, Eq. we 
get for the phase factor ratio due to the 3-simplices 

-r^ = ChAh2,a.i)ch2{hi,b^^^), (C3) 
^21 

with i the lattice site of intersection. 

Since the definition of statistical angle only makes 
sense for Abelian quasiparticles, we consider only trivial 
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a) 






FIG. 21: Quasiparticle statistics from ribbon commutation, (a) Commutation of two strings (ribbon operators, black and red) at 
the intersection point (blue square) gives the braiding statistics. (Top) Strings representing tunneling of particle — anti-particle 
pairs across the system (the surface of torus) can give the self-statistics. (Bottom) Braiding setup independent of system shape, 
(b) Near the intersection of two ribbon operators w^^i^gi), Wf^^{g2), the phase difference occurs due to 3-simplices on top of 

two blue shaded triangles, detailed in (c) with phases Wmn = (f| i^h^ (Sm)w'Jj" (pn) |i) shown in top {W12) and bottom (VK21) 
rows. Note that elements bi are defined to be oriented oppositely to our lattice edge orientations, but consistent with definition 
of ribbon operator phase, Fig. 1131 hence the appearance of in Eq. (jCSj. 



cocycles Ch- This means that we can rewrite the 2-cocycle 
using a 1-cochain, as in Eq. ([57]) . 

In that case, the phase ratio 612/621 is given by 
Eq. (j74p . and can easily be simplified, giving 



612 
621 



£h2 {bi 



(C4) 



see labels in Fig. [2]Ta). 

Using the property £h{gi^^ 
^h^ {91) / ^h^ {92) , which holds for any 1-cochain de- 
scribing a trivial canonical 2-cocycle, the total phase 
ratio ^ = l^fl^ finally becomes 



-i2e 



^ £h^{9l)£h2 {92) 

£hd9'i)£ 112(92)' 



(C5) 



The group elements appearing here are: gi = Yij^j^ 
92 — Yik^k, i-e. the values for isolated strings, and 
91 ~ '^sJlj-'^i' 92 — ^lYik^k, i-e. the values after the 
other string is already applied. 

After using the 1-cochain, Eq. dHT]), on Eq. PHO)) . the 
result obtained by calculating the braiding matrix by 
using the algebra of local operators at the string ends 
almost exactly matches the one here (up to replacing 
9i — ^ ^2 '51 here). The difference involves a deeper subtle 
issue we will address now. The obtained angle 26 (either 
following from Eq. (|100l) or Eq. (|C5p ) should represent 27r 
braiding of Abelian quasiparticles in systems with trivial 
cocycle Ch- However, for the statistical angle to be phys- 
ical, one expects that it is invariant under changes of the 
3-cocycle by a 3-coboundary, as introduced in Eq. (ITUl) . 
because this change does not alter the physical content of 
the theory. Such transformations of by a coboundary 



lead to the following transformation of the 1-cochain e 
by an arbitrary function u{x, y): 



Shig) -> £h{g) 



u{h,g) 
u{g,h)' 



(C6) 



One can see that the factor e/n (ft.2) e/12 e^i) in the result 
Eq. (|C5I) is invariant, but the rest of the expression is not! 
(The same is true for result following from Eq. (|100p .) 

The resolution of this puzzle is instructive. Namely, 
the ribbon operator F*^'''f'(r) we used here and in Sec- 
tion |TVDl when acting on the ground state, does not nec- 
essarily create a quasiparticle pair state with well defined 
flux h and gauge charge hg. As explicitly shown in the 
examples of Section |Vl the gauge charge operator may 
act non-trivially within the quasiparticle Hilbert space. 
(This Hilbert space is created by the action of _F(^'S)(r) 
on the ground state, but has to be specified further by 
discriminating different values of SG elements uc at the 
ends C = ^, -B of string F, Section IV Al ) This means 
that physical states, having well-defined gauge charge, 
can actually be non-trivial linear combinations of states 
spanned by i^('*'9)(F). 

Of course, the braiding matrix Eq. (|100p contains all 
information necessary to specify braiding properties of 
quasiparticles, one only has to pose the right question, 
which would involve explicit use of physical flux and 
charge states constructed using the formalism in Sec- 
tion El 

In the present discussion of Abelian quasiparticles, it 
becomes obvious that the physical quasiparticle states 
with fixed flux and charge are simply obtained by absorb- 
ing the non-invariant factors in Eq. (jC5p . This is achieved 
easily by using ^('''^'(F) = FC^'^) (F)e,j(5), which leads 
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to the statistical angle independent on the cocycle within 
a fixed cocycle equivalence class: 

e-'^' ^ ehdh2)ehM. (C7) 

(The same is obtained starting from expression given by 
Eq. mJh .) 

Let us now use Eq. (jC7p on two Abelian topological 
theories with Z2 order, namely the toric code (TC) and 
double semion (DS) models (see Appendix ID|) . These two 
orders are physically distinguished exactly by the statis- 
tics of their quasiparticles. 

As demonstrated explicitly in Appendix [Dl the TC is 
recovered by choosing the trivial cocycle w (51, 52, 53) = 
l,V(7i,52,53 S G. From these constraints we need to 
obtain the values of the 1-cochain eh{g)- The constraint 
^{91,92,93) = 1 leads to: 

c„(5i,52) = l, (TC), (C8) 

which implies that 

eh{gi)sh{92) = eh{9i-92), (TC), (C9) 

i.e. the 1-cochain e actually becomes a 1-cocycle. There 
are two representations of the Z2 group, i.e. in total four 
solutions for e given by eo{g) = 1 or eo(g) = (—1)^, and 
ei(5) = lor£i(g) = (-l)9. 

In the present formalism, a particular ribbon is deter- 
mined by its flux hy and one of the particular solutions 
for the 1-cochain e. 

The TC has four quasiparticles, {U, e, m, em}, the 
trivial, charge, flux, and charge-flux bound state, re- 
spectively. Assigning hy = to Il,e, and = 1 to 
TO, em, we make a choice (not unique) of assigning each 
of the four solutions to one particular quasiparticle: 1) 
To 1: eh{9) = 1; 2) To e: Shig) = (-1)^; 3) To to: 
ei(.g) = (-1)9; and 4) To cto: eo{g) = (-1)^ . 

This solution recovers the well-known TC result that 
U, e, TO are bosons, em a fermion, and non-trivial mutual 
statistics is given hy 29 = tt amongst the {e, to, em}. 

According to the explicit demonstration in Ap- 
pendix ini the DS theory is obtained by choosing the 
second representative 3-cocycle for G — Z2 given by 
^{91,92,93) = -1 for 51 = 32 = ffa = 1, and 1 other- 
wise. Repeating the analysis from the TC case above, 
the difference lies in the constraint: 

ci(l,l) = -l, (DS), (CIO) 

which leads to solutions given by £0(5) = 1 or £0(5) = 
(-l)f, and £1(1) = i or £i(l) = -i, while £i(0) = 1 
always. 

The non-trivial excitations in the Abelian topological 
phase described by DS theory are two scmions si, S2, 
and their bound state S12. Assigning — to l,si2, 
and hy = 1 to si,S2, we make a choice (not unique) 
of assigning each of the four solutions to one particular 
quasiparticle: 1) To 1: £i(l) = i; 2) To S12: £o(ff) = 




FIG. 22: Recovering the two Z2 topological phases from the 
exactly solvable models: toric code (TC) and double semion 
theory (DS). The dual lattice is honeycomb (dashed). Red 
dots mark the positions of group elements Qij assigned to 
edges ij, becoming "spin-1/2" states in TC and DS. A honey- 
comb lattice plaquette ph is shaded. The plaquette operator 
Bp (Fig. [9)1 acts on the six "spins" of ph. The DS model dif- 
fers from TC by having an additional phase in its plaquette 
operator. This phase depends on spins on the six outer legs 
of ph, one of which is marked by green line. These six "outer" 
spins belong to the six tetrahedrons in our Bp, Fig.|9] The Qt 
operator in both models just acts on the three spins nearest 
to a site (spins on blue lines). 



(-l)f, £i(l) = -i; 3) To Si: £o(g) - (-1)^, £i(l) = i; 
and 4) To S2: £i(l) = -i. 

This solution indeed recovers the well-known quasipar- 
ticle properties: l,si2 are bosons; si, S2 are semions 
{29 — tt); and the non-trivial mutual statistics is 29 = n 
between S12 and either si or S2- 

Appendix D: Explicit form of models for Z2 
topologically ordered phases 

Here we consider the well-understood case of Z2 topo- 
logical order, to demonstrate that our general model ex- 
plicitly yields models for the two inequivalent phases hav- 
ing such order. Two well-known models for the two dis- 
tinct Z2 topological phases are Kitaev's "toric code"— 
(TC) and the "double semion" theory^" (DS). We will 
make a direct comparison to the TC and DS model vari- 
ants presented in Refl27l 

For that purpose, we set the symmetry group to be 
trivial, SG — Z\, and the gauge group to GG = Z2. 
There are only two group elements, which we can label 
by G = {1 = -1-1, a = —1}, with a = a~^. The group ele- 
ment gij = ±1 assigned to the ij edge (of the triangular 
lattice) we will call "spin" (±l="up/down") positioned 
at the mid-point of the edge ij. We will now switch to 
the dual lattice, which is honeycomb, as shown in Fig. [521 
The spins are still on the edges of the honeycomb lattice, 
and the honeycomb plaquette ph contains the same six 
spins as our original plaquette p. 

Let us start with the operator Qt, which will be the 
same in both topologically ordered theories. This oper- 
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ator forces the product of three spins adjacent to a hon- 
eycomb lattice site to +1, see Fig. (25] It therefore acts 
as the projector 

on three spins i neighboring the vertex v. Comparing 
directlyto the variant of TC and DS models presented 
in RefOT, the first term in H, Eq. ([57]) . becomes up to 
an overall constant just the standard vertex contribu- 
tion Y.V Yliev to H'^^ or in ReflU Further, 
the product OiGp appearing in the plaquette terms of 
H exactly becomes the Il-uep^i term in H'^'^ or H^^ 
which projects the flux in the honeycomb plaquette ph 
to zero (this is the Pp factor in Refl27l). 

At this point, we can discuss the physical interpre- 
tation of Qu and the outcome of this subsection more 
precisely. In our model, Qy = l,Vi; gives the zero-flux 
rule on the lattice. One should recall that if we had been 
considering the dual of the globally symmetric SG — Z2, 
GG — Zi, the zero-flux rule would have been automat- 
ically satisfied (see Eq. (H^l) ). This fact is equivalent to 
saying that in the undualized theory the domain walls 
separating regions of up and down spins have to form 
closed loops, i.e. "closed strings" in the language of 
string-net models^. No matter if we start from the du- 
alized SG = Z2 model or the GG = Z2 model, the re- 
striction to states satisfying = 1 (which is automatic 
or chosen, respectively) leads to Z2 gauge theory mod- 
els, and this will be clearly shown below. In fact, Ref. [13 
clarifies that there are still two physically distinguishable 
Z2 gauge theories, which represent two different topologi- 
cal orders; both of these were obtained there by dualizing 
Z2 spin models, automatically enforcing the closed string 
rule. We will explicitly show that our G = GG = Z2 
model indeed reduces to the two inequivalent Z2 gauge 
theories when = 1. 

One can now ask how do the Kitaev's "toric code"— 
(TC) and the "double semion" theory^o (DS) fit in this? 
Simply, the two different Z2 topological orders are also 
described by these two different models. (TC and DS can 
be seen as Ising matter coupled to the two different Z2 
gauge fields. ^^) The TC and DS models are expected to 
arise when the restriction to = 1 is lifted, i.e. they 
are "open string" models. We will explicitly show below 
that our G — GG = Z2 model indeed gives rise to TC in 
the absence of that restriction. 

Let us move on to the plaquette operator and the 
Hamiltonian of our model. Under the six group 

elements g on edges (of triangular lattice) which share 
the lattice site i are multiplied by a. Therefore the 
six spins on the honeycomb ph plaquette are "flipped" , 
(7 = ±1 ^ =pl, i.e. acted on by Pauli matrix ax- Obvi- 
ously, B^=i does not change the state (no flipped spins). 

Next we consider the phase factor of 13. There are 
exactly two inequivalent classes of 3-cocycles w for the 
group G = Z2, as dictated by H^{Z2, U{1)) = Z2. They 



are represented by the two distinct choicesas: 

L^(-l,-l,-l) = ±l. (D2) 

In both cases, a;(5i, ^2, .gs) = 1 for any 31,32,53 e G 
when they are not all equal to —1. 

To obtain a model for TC phase, we choose a trivial 
cocycle, i.e. ^(31,32,53) = 1, Vgi,g2,53 e G. The ampli- 
tude from Eq. ([5^ then becomes trivial, -Bj*^ = 1. The 
action of the plaquette operator is 

with six spins I in the honeycomb plaquette ph. This 
is equivalent to the standard plaquette operator of TC 
model in Ref. [13, although in the present case there is 
a constant term added to the a^^ term. This constant 
shift can be removed by smoothly deforming the B^'^ 
operator without ever closing the gap in the Hamiltonian. 
(The constant shift is there to make our Sj'-^ a projector, 
which is its general property.) Our model with cj = 1 is 
consequently equivalent to the TC model. Under the 
Qu = 1 constraint, obviously we obtain the standard Z2 
gauge theory on the honeycomb lattice. 

We now turn to the DS model^°. Compared to the 
TC model, the only difference is that the honeycomb 
plaquette operator now has an additional phase factor 
determined by the value of six spins on "outer legs" of 
the plaquette^''; one of six such legs of the plaquette ph 
is marked by green line in Fig. [2H This phase in the DS 
model is assigned by the rule: if exactly two, or exactly 
six, of the "outer leg" spins are in state —1, there is a 
phase factor —1 in the plaquette operator. If we again 
restrict to a "closed string" theory, i.e. Qu — \, the lines 
of —1 spins on the honeycomb lattice physically repre- 
sent the closed domain walls of the dual Z2 spin theory, 
and therefore there is always an even number of —1 legs 
entering a plaquette ph. As mentioned, this closed string 
restriction of DS model is the second type of Z2 lattice 
gauge theory. 

Let us show that under = 1 our model indeed re- 
duces to the Qv = 1 restricted DS model from Ref|13. 
The six "outer leg" spins of the honeycomb lattice are 
actually lying on the six outer edges of the plaquette p 
in our model on the triangular lattice (compare Figs. [22] 
and [21). Although our Bp does not change them, they are 
on the edges of the six tetrahedrons which determine the 
phase of the operator. 

Given the restriction Qy = 1, i.e. the zero- flux rule on 
the triangular lattice, and using uj{—1, —1, —1) = —1, it 
is straightforward to show that the phase coming from six 
tetrahedrons (see B^"^^ in Eq. ([5T|l ) is equal to —1 only 
in two cases: 1) When the only outer leg spins in state 
— 1 arc on two legs separated by 120°; 2) When beside 
the two legs as in case (1), there are two more legs with 
spin —1 separated by 60°. 

These rules seem strange, however, this is remedied 
by a simple change of basis. As explained, for an 
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arbitrary basis state there is an even number of —1 
outer legs, i.e. segments of dual domain walls, en- 
tering the honeycomb plaquette ph. These segments 
of domain walls have to be connected (they do not 
have ends!) in some way along the six edges of the 
hexagonal plaquette ph. Given the configuration of 
— 1 legs entering the ph, no matter how we choose to 
connect them, the number of ph edges we will use for 
that will have fixed parity. It is therefore well-defined 
to use this number, #int.legs(p/i), for any basis state. 
After redefining the basis by [state of plaquette p) — > 
(_l)#int.iogs(ph) ig^g^^g plaquette p), our "strange" 

phase factors become exactly the DS model's phase 
factors described above. Therefore, within the zero- 
flux manifold of states, our model with G = Z2 and 
1,— 1,— 1) = —1 exphcitly reduces to the "closed 
string" version of the DS model, which was firstly 
understood in that form in Refl27l. 



Appendix E: All two-particle states are given by 
action of ribbon operator, and local operators act 
projectively on these states 

Let us first show that the local operators form a pro- 
jective representation of the group G in the Hilbert space 
C{A, B) of two excitations at the ends A, B of the rib- 
bon r. Let us for concreteness focus on the projected 
subspace Cu^ {A, B) defined in Eq. (j84| , and the exact 
same results follow for Cub{A^B). For brevity, let us 
also use the double index notation, i.e. i = {hi,gi),j = 
(hj, gj), . . ., with hi,hj, . . . G GG and gi,gj, . . . G G. Di- 
rectly combining the definition of excited states Eq. (|84|) 
and the operator algebra from Eq. (|8ip . we get: 



DM)\^uJ=k,,h-^%Hhk,h-')c-^{g„g,) 



DjiB) IV'^^) = Sh^j^^ chk{9j ^gk,gj) 



WuA 



(El) 



(E2) 



where we used the usual factorization gj = h'j ■ gj with 
gj G SG, h'j G GG. Applying the operators twice and 
using the 2-cocycle identities Eq. (j87bp . (|87ap . we get: 

D,{A)D,{A) I^^J = c,-^{g,,g,) D,,{A) (E3a) 

D,{B)D,{B) I^^J = CHA9^,9,) D^,{B) , (E3b) 

where we used the obvious shorthand notation ij = 
{hihj, gigj). The 2-cocycle c/j, with h G GG, therefore 
determines the projective representation of G. 

The projective representation of local operators is ac- 
tually unitary, as we will next show. In the basis ["iAM^) 
labeled by k (ua is fixed) , the matrix elements of opera- 
tor £'(h.g) (A) from Eq. (|ET|) are 



Using the multiplication law, 
-1 



Eq. (|E3al) . it further 
follows that M{h,g)M{h-\g-'^) = c?,(g, g-i)M(l, 1), 
where the identity operator M(l, 1) in the basis 
is represented by the unit matrix. This determines the 
inverse of the matrix M{h, g). 

On the other hand, in the orthonormal basis the matrix 
of adjoint operator is given by the conjugated transpose 



t 

k'k 



Mih,g)*'- 



The matrices M''' and M ^ differ only by a phase factor 
containing five 2-cocycle factors. Using the cocycle rules 
Eqs. (j87b|) . ((57a|) and (157c)) it is easy to show that the 
phase factor is equal to one, and the representation M is 
unitary. The same can be shown for the representation 
ofD{B). 

We note that the projective representation of G formed 
by local operators D{A) or D{B), as discussed here, is not 
the algebra defining the quasi-quantum double. The op- 
erator algebra within the quasi-quantum double is formed 
by the "topological operators", and depends on cocycle 
of elements in GG only. This is discussed in detail in 
Appendix |B] 

Let us next consider the problem of determining the 
two-quasiparticle Hilbert space C{A, B), where we fix the 
two excitations at positions A and B. We will prove that 
the Hilbert space C{A,B) — ®uc^uc{AtB) (C is either 
Aor B), as introduced after Eq. ([5^ . indeed contains all 
the states of the two excitations, i.e. it contains C{A, B). 
Our proof follows the one given in Ref. [5^ for the model 
generalizing the Z2 toric code to arbitrary finite groups 
(the "generalized toric code"); however, our proof has 
significant changes, rooted in the fact that our model 
differs significantly from the generalized toric code. (For 
the present proof, the presence of global symmetries in 
our model makes the most important difference.) 

We note that since clearly the two-quasiparticle Hilbert 
space C{A, B) is a subspace of /C(r), the action of -F'(r) 
and D operators appearing below is always well-defined. 

To establish our result, we need to define additional 
useful tensors related to the ones in Sec. IIV CI We 
recall the double index notation using Latin indices: 



i = {hi, 
and gi,gj 



with hi , h 



11 ■ 



G GG 



G G, while the Kronecker delta function 
g. . The antipode tensor S is defined such 
Dk is the inverse of the element 
Di in the algebra F (g) D. This definition results in 
the following identity and solution for S: 



that the element S^F'- 



St'A 



I "fcn"m — ^P"^ 



(E5) 
(E6) 



One can check that the defining identity Eq. (|E5P holds 
for the given form of S by using the 2-cocycle identities 
Eq. ((87b)) and Eq. (|87d)) .^ 

Recall from Sec. IIV O the functions = (5^ . and e' = 

M{h,g)k',k ^ Cg^{h^^ ■,h)c^li{g,gk)5g^, ,ggj^j^-i5hk.hf^, ■ '^g,,!: which define the unit and counit of the algebras 



(E4) 



Using these functions we now define 
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and 



C{A) 



GG 



(E7) 



where Nqg = \GG\ is the order of the gauge group and 
No = \G\ is the order of the group G. It is impor- 
tant to notice that by its definition, the constrains 
only the gauge group element hi = t oi the double in- 
dex i = (hi,gi). This means that the operator C{A) in 
Eq. (jE7l) projects out any non-zero flux in triangle Ia, 
i.e. it ensures there are no flux excitations at that lat- 
tice plaquette. However, the action of operator D{A) in 
C{A) is still non-trivial, e.g. since the element £ G in 
principle modifies the lattice site element u^^. 
We now arrive at the central identity: 



NggNg 
1 



(E8) 
(E9) 



-1 



NggNg 

where we again used the double summation convention, 
i-^- = J2h eGG g eG' ^^^^ equality can be 

proved by using the explicit definitions of all tensors, the 
operator algebra, and the 2-cocycle properties. The last 
line of Eq. (jESp is slightly more subtle. Consider an ar- 
bitrary matrix element of the operator in Eq. (jE9p . The 
delta functions ensure that the ribbon and local operators 
do not change any edge or site elements between the ini- 
tial and final state. Further, their phase factors contain- 
ing cocycles are also trivial due to our choice of standard 
cocycle, see Eq. (|87cl) . Consider then the matrix element 
contribution from the ribbon operator: S{uij^uY^ ,gr), 
with the obvious factorization r = {hr,gr) — {hr, h'^. ■ gr), 
hr,h'^ € GG, cjr G SG. But there is a summation over 
r, and therefore over cjr £ SG, so the delta function is 
harmless. The same is true for the delta function com- 
ing from the local operator. We therefore see that the 
quantum amplitude is 1, and indeed the action equals 
the identity operator 1 in the algebra Da (g) F (g) KDb- 

Let us now consider an arbitrary state |^) in the 
Hilbert space C{A,B). We can define: 



\7f) ^ C{A)F^{T) \^) 
G, = NGGNGTsnipSPF^{T) 



(ElO) 



so that 



|V) = GJ77') (Ell) 

holds due to the identity Eq. (|E8|) . The states \7f} ac- 
tually have no flux excitations at A due to the action of 
C{A). Elementary gauge and charge excitations must be 
in pairs and by construction are symmetric, so \r]'') can 
only describe the ground state \gs) or a one-particle state 
that breaks the global symmetry. Since the site elements 
Uij^,Uig are the local (at A,B) degrees of freedom acted 
on by the global symmetry, we conclude that the \rf) 
states can be written as linear combinations of the form: 



Ui,U2GSG 



(E12) 



with some coefficients fc^i.ua 



(These states include the 

ground state itself." 

Eqs. (|E11[) and (|E10[) show that an arbitrary two- 
particle state is a combination of the form 



K„ 



,F"^{T)PuM)PuAB)\9s) . 



ui ,U2GSG 



(E13) 

By the definition of the ribbon and projector operators, 
the \tp) is the zero vector unless ui ■ — g„i, where 

m = {hm, gm) ~ {hm, ' ^mj /^m € GG , (jm G SG. 

One of the three sums over SG in Eq. (jE13p is therefore 
superfluous for physical states in £{A, B). 

We can use ui ■U2^ — g„i to eliminate either of the ui^2 
in Eq. (jElSp , which immediately shows that the arbitrary 
state G C{A, B) is indeed a linear combination of 
projected ribbon states from Eq. (|84|. 



Appendix F: Symmetry fractionalization for 
multiple quasiparticles 



In the examples of Sec. lVBl we describe the scheme to 
find the fractionalized symmetry transformations U„{G) 
[G = A, B) for a quasiparticle at location C in the single- 
ribbon states, see Ea. (|108p where we assumed Eq. (|110p . 
In this section, we show that at least when the non-trivial 
symmetry fractionalization only occurs in a Z2 gauge sec- 
tor (namely for visons in a gauge sector Z2 G GG), the 
definition of local operators Ua{G) in Eq. (|108l) imple- 
ments the fractionalized symmetry transformations for 
multi-particle states created by ribbon operators defined 
for ribbons ri,r2, ...,rjv that all share the same end-S: 



uMi)uM2)-UMn)u„{b)\^':^'^- 



') 



(Fl) 



In our simple examples in the main text, it is al- 
ways true that when non-trivial symmetry fraction- 
alization occurs, it always only occur in one Z2 
gauge sector. In addition, we belive that all states 
{|V'mb'°^ '^")}' Vms, fci, A;2, fcjv span a Hilbert space 
£(^1,^2, ...,Am,B) that contains all possible excitations 
at those locations. Therefore Eq. (jFl[) . which we will 
prove below, indicates that the Ua{C) operators are the 
general fractionalized symmetry transformations in our 
examples. 

It is straightforward to show that (both GG and SG 
are Abelian): 



fclfc2 



(F2) 



while 
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D 



{h 



-h2,a-K,UA.u-^--h„,cr-\uA^u-^' ■ ({Q, hi, h2, ...hN) ■ H^^'-?;^-^ ") 
^ihu<y-\uAi<^h2.a-\uA2-^hr,,<7~\uAr,^hih2...hM.,<T-\ug^ ■ £,{g,hi,h2,--hN) ■ 

fcl ^2 . . . 



h-. .a ^.UA^ /in 



'"/ll/l2---^iV,<5" ■'■ 



_^^.^(5,/ii,/.2,...M-l€-t;") 



r 



(F3) 



where we used the assumption: ^h.a-^,h' = li ^h^h' £ ^x,y,z- The extra phase ^(g, /i2, .../iat) is defined as: 
GG and a G SG, and the basic properties of the tensor 

I 



$(5, /ii, /i2, .../ijv) = c^-i(/ii \ ^i)c^-i (/i2 \ 'i2)---c^ii ihj/,hN) ■ c„-i {hi,h2h3...hN)c„-i {h2, h3...hN) — .c^-i (ft.jv_i, /iat) 

(F4) 



r 



In our examples, non-trivial symmetry fractionalization 
occurs for a single Z2 subgroup ■'^ e GG and GG — 
X GG' . Let us denote a group element h G GG 
as /i = {hsF,h'), with /isF € and h' e GG' . This 
indicates that C(j-i (/iq, ft-b) 7^ 1 can occur only when both 
ha,SF and hi, sF have non-trivial components in Z^^, 
i.e., ha,SF = ^h.SF = 1, and generally c„-i{ha,hb) = 
Ca-^ {K.SF, hb,sF)- Let's denote c^-i (ft-a.sF = 1, ^b,SF = 
I) = rj (we use the {0, 1} notation for Z^^.). Under these 



conditions, it is easy to show that ^{g,hi,h2, ■■■hN) = 
^-Nsp/'i^ where Nsf is the total number of visons (at 
the end-Ai and end-B) in the ^ sector for the state 



h\ k2 . . . fcjv 



). Therefore, using Ua{C) defined in Eq. (|TII5)) 



for an excitation with a gauge flux h, 



U„{C) = ^^^MhJF^) ■ ■ D(^u,a-^){C), (F5) 

it follows that Eq. (|Fll) is indeed satisfied. 
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does not modify the argument below. 
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One may wonder how we fix our convention — because for 
a usual finite abelian gauge theory GG in 2-1-1 dimensions, 
the gauge charges and gauge fluxes are self-dual. In fact, 
for a general gauge theory with finite gauge group GG, the 
gauge fluxes and gauge charges are physically distinct. For 
example, gauge fluxes are labeled by conjugacy classes of 
GG, while gauge charges are labeled by the irreducible rep- 
resentations of the centralizer group of a conjugacy class. 
We will show that in our exactly solvable models, it is 
clear that the H^{SG,GG) classiflcation corresponds to 
the symmetry fractionalization classes of the gauge fluxes. 
In Refl27l. string-net models are used to construct the ex- 
actly solvable models for Dijkgraaf- Witten theories. It ap- 
pears to us that this construction may not be general; 
namely not all Dijkgraaf- Witten theories can be described 
by string-net models using the construction in Ref. [60l . In 
the present work, we use a different way to construct ex- 
actly solvable models for Dijkgraaf- Witten discrete gauge 
theories, which is general. 
^® see, for instance, the webpage by John Baez at: 



http://math.ucr.edu/home/baez/qg-winter2005/w05week06.pdf 
For simplicity, let's consider G — Z2 and use the trivial 2- 
cocycle oj{x,y) = 1, Va;,?/ £ Z2. We can construct the pro- 
jector P for a circle following the discussion in Sec lIII A2\ 
after triangulating it by A'^ line segments. Let's use gi — ±1 
as the group element for the colored line segment /. The 
2-fold ground state sector Img{P) of the induced TQFT 
can be easily found: = 1/2(^-1)/^ ^^^^^^^^ \{gj}), 

and IV2) = 1/2(^-1)/^ I^n.S7 = ~i Naively the sec- 
ond ground state corresponds to a trapped Z2 gauge flux 
inside the circle. However, the degeneracy between 
and |^/)2) is not protected. One can consider a local pertur- 
bation SH = e\gi = —l){gi = +1\ + h.c. for a certain line 
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segment I. Straightforward perturbative calculation shows 
that this perturbation lifts the ground state degeneracy by 
a finite energy gap. 

As usual with duality, the flux, i.e. Y\^- Uij, automatically 
vanishes through a closed loop ij on the lattice. The dual 
SG therefore has no flux particles. Also, note that the du- 
ality transformation is not a 1 — to — 1 mapping, since con- 



figurations {ui} and {ui ■ s ^} are mapped to the same 
configuration Uij in SG, with SG understood as the group 
SG acting in the dual theory. 

Although we do not have a general proof, we believe that 
this is true, which is also confirmed in the examples in 
Sec|V]that we solved. 



